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ABSTRACT 


This thesis is concerned with derivation of a macroscopic plasma 
Lagrangian, and its application to the description of nonlinear three- 
wave interaction in a homogeneous plasma and linear resonance oscilla- 
tions in a inhomogeneous plasma. 

One approach to obtain -the— Lagrangian is via the inverse problem 
of the calculus of variations for arbitrary first and second order 
quasilinear partial differential systems. Necessary and sufficient 
conditions for the given equations to be Euler- Lagrange equations of a 


Lagrangian are obtained. These conditions are then used to determine 


the transformations that convert some classes of non- Euler- Lagrange 


equations to Euler-Lagrange equation form. The Lagrangians for a linear 
resistive transmission line and a linear warm collisional plasma are 
derived as examples. 

Using energy considerations; the correct macroscopic plasma 
Lagrangian is shown to differ from the velocity- integrated Low Lagrangian, 
by a macroscopic potential energy that equals twice the particle thermal 
kinetic energy plus the energy lost by heat conduction. The generalized 
variables are the macroscopic plasma cell position (Eulerian coordinates) 
defined in Lagrangian coordinates, and the vector and scalar potentials 
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defined in Eulerian coordinates. With the continuity and heat flow 
equations treated as constraints, the Euler-Lagrange equations are shown 
to be the force law and Maxwell’s equations. The effects of viscosity, 
heat conduction, and elastic collisions are included in this variational 
principle. The corresponding macroscopic Hamiltonian, and the micro- 
scopic Hamiltonian corresponding to the Low Lagrangian, are also derived. 

Under the assumptions of scalar pressure and adiabatic processes, 
the macroscopic Lagrangian is approximated by expansions in weak pertur- 
bations of the generalized variables. The averaged Lagrangian method 
is then used to derive nonlinear three-wave coupling coefficients in a 
warm homogeneous two-component plasma. The effects of wave damping 
are included phenomenologically in the coupled mode equations. The 
general results are then specialized to make detailed quantitative 
comparisons between theory and available experimental results on parame- 
trically excited ion-acoustic waves. 

The approximate quadratic Lagrangian is also used to estimate the 
electrostatic ( Tonks-Dattner) resonance properties of an inhomogeneous 
plasma. The Rayleigh-Ri tz procedure is applied directly to the Lagrangian 
corresponding to a system of Euler-Lagrange equations. Use of an appro- 
priate set of trial functions then leads to frequency and eigenfunction 
estimates in excellent agreement with the existing theoretical and experi- 
mental results for a low pressure positive column. Since this method 
mainly involves evaluating finite integrals, and solving algebraic eigen- 
value equations, it is found to be more efficient than numerically solving 
differential equations, and more accurate than inner-outer expansions. 
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1 . INTRODUCTION 


In theoretical descriptions of plasmas, three approximate models 
are commonly employed. These are the cold, .the microscopic, and the 
macroscopic plasma models, all of which use Maxwells equations. To 
complete the system of equations, the cold plasma model uses Newton’s 
force law, 


m ~ = q(E + v X B) , 
at ~ ~ ~ 


(1.1) 


where E and B are the electric and magnetic fields, and v, m, and 
q are the particle velocity, mass, and charge of a species. The 
plasma is then regarded as consisting of interpenetrating cold fluids 
with charge and current densities, 2 qn and 2 qnv, respectively, where 2 
sums over all particle species, and n is the particle number density. 

The microscopic plasma model uses Newton's force law in (l.l) for each 
particle, and the Boltzmann -Vlasov equation (Clemmow and Dougherty, 1969) 


dv 

+ v • Vf + Vf = 0 

St ~ dt V 


where f(x, v,t) is the Boltzmann distribution function for each particle 
species. The expressions for charge and current densities now become 
£ qjfdjv and 2 qjfv^dv^ respectively. The plasma is regarded as a system 
of charged particles evolving under the influence of their own electro- 
magnetic fields, and externally applied fields (if any). In principle, 
complete solutions to the particle force law for all particles will auto- 
matically generate the solution of f(x,v^, t) because (l.2) is a statement 
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of particle conservation along the particle trajectories in a six- 
dimensional phase space (x, v) . When it is not necessary to obtain the 
particle trajectories, solutions of f(x,v,t) are obtained by use of 
(l. 2 ) together with Maxwell's equations* 

The macroscopic plasma model uses velocity moments of the Boltzmann- 
Vlasov equation (l t 2 ). Examples of these equations are expressed by 
(3.29) and (3A6)* The charge and current densities are now written as 
£ qn and £ qnv^, where is the drift velocity, i.e. the averaged 

local velocity of a particle species [see ( 3 « 13 )]« Here, the plasma is 
approximated in terms of localized variables such as density, drift 
velocity, pressure, and heat flux, In terms of degree of approximation, 
the macroscopic model falls between the other two models. 

The appropriate Lagrangians for the cold plasma model (Galloway and 
Crawford, I9TO) and the microscopic plasma model (Low, 1958) are already 
well-known 0 However, the Lagrangian for the macroscopic plasma model, 
that corresponds to Maxwell's equations and the moments of the Boltzmann- 
Vlasov equation, has not been established. The gap will be filled in 
this thesis. 

The interests in developing this variational principle stem from 
the fact that current theoretical investigations in nonlinear wave-wave 
and wave-particle interaction properties of homogeneous plasmas, and in 
linear properties of inhomogeneous plasmas, are at the limit of analytic 
tractabili ty . While a suitable variational principle does not provide 
new fundamental laws, it leads to a relatively concise formulation and 
easy manipulation for these otherwise difficult problems. 

To obtain the suitable macroscopic Lagrangian, the first problem that 
arises is the inverse problem of the calculus of variations, i.e, the 
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derivation of Lagrangians from arbitrary equations. This mathematical 
approach will be examined in Section 2 in contrast to the approach via 
energy considerations commonly used for physical problems. From this 
general point of view, it will be shown that energy dissipation effects 
can be included in variational (minimal) principles, in general, and the 
results will be demonstrated with examples. 

The approach of the inverse problem treats all dependent variables 
as generalized variables. The plasma variational principle to be presented 
in Section 3 will treat only the macroscopic plasma cell position and the 
electromagnetic potentials as generalized variables. This type of formula- 
tion for a system of discrete charged particles, was described in a 
relativistically covariant form by Landau and Lifschitz (1965), and 
in the non-relativistic form by Goldstein ( 1950 ). Extensions of these 
Lagrangian densities to the microscopic plasma model, to include a velocity- 
distributed system of particles,, have been proposed by Sturrock (1958a) 
and Low (1958). Based on Low's Lagrangian, and using energy considerations, 
we shall obtain the corresponding Lagrangian and Hamiltonian for the 
macroscopic plasma model, with the effects of viscosity, heat conduction, 
and elastic collisions taken into account. 

In applications of Lagrangians to problems involving homogeneous 
plasmas, there has been progress in the areas of linear waves (Kim, 1972 ), 
nonlinear three-wave interactions (Galloway and Crawford, 1970), wave- 
background interactions (Dewar, 1970), wave kinetic equations ( Suramlishvili, 
196U and 1965; Galloway, 1972 ), higher order nonlinear wave processes 
(Dewar, 1972 ; Dysthe, 197*0 7 and statistical analysis of plasma turbulence 
(Kim and Wilhelm, 1972 ). For problems in inhomogeneous plasmas, results 
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have been presented in the area of energy principles (Newcomb, 1962), and 
in the use of the Rayleigh-Ri tz procedure to obtain approximate solutions 
(Dorman, 1969) . As compared with other branches of physics, these cases 
comprise a disproportionately small fraction of the theoretical effort 
in plasma physics. 

In this work, the new macroscopic Lagrangian will be used in two 
ways. In Section 4 we shall be concerned with the general description 
of nonlinear three-wave interactions in a homogeneous plasma by use of 
the averaged Lagrangian technique (Whitham, 1965) . These results will 
be specialized in Section 5 to parametric amplification of ion-acoustic 
waves to make quantitative comparisons with available experimental data. 

The second application of the Lagrangian will be pz’esented in Appendix B, 
where the Rayleigh-Ri tz procedure is applied to obtain approximate solu- 
tions for electrostatic resonances in a low pressure positive column. It 
will be shown that the results compare favorably with available experimen- 
tal data and conventional numerical calculations by others (Parker, Nickel, 
and Gould, 1964 ) . 

Some conclusions are drawn in Section 6 , where new contributions and 
future extensions of this research are briefly discussed. 



2 . INVERSE PROBLEM OF THE CALCULUS OF VARIATIONS 


2.1 Introduction 

In the calculus of variations Hamilton’s principle is applied to 
the integral of a given Lagrangian density, often referred to. as the 
action integral, to obtain the Euler-Lagrange equations extremizing the 
action integral (see for example, Courant and Hilbert, 19 66). The 
inverse problem of the calculus of variations is to find the conditions 
that arbitrary differential equations must satisfy to be the Euler-' 
Lagrange equations of a certain Lagrangian density, and to determine 
an appropriate Lagrangian density from the given differential equations. 
The term ’Lagrangian density’ as used here will include cases with 
multiple independent variables. Recent developments of a method of 
studying weakly nonlinear, wave propagation in distributed systems makes 
the inverse problem for a system of partial differential equations of 
particular interest (Whitham, 1965; Galloway and Crawford, I97O; see 
Chapter . A general approach is required to obtain appropriate Lagran- 
gian densities for systems of equations including effects of energy loss 
due to heat flow, viscosity, and collisions. 

The inverse problem of the calculus of variations attracted 
attention over a century ago, when Jacobi ( 1837 ) examined the character- 
istic properties of an ordinary Euler-Lagrange differential equation of 
second order (Kiirsch^k, 1906; Akhiezer, 1962) . More involved problems 
have since been studied. Such work includes that by LaPaz (I93O), who 
treated the case of one dependent variable in many independent variables 
using the necessary property of self-adjointness of the Euler-Lagrange 
equation. Douglas ( 19 ^ 1 ) has studied the case of many dependent variables 



in one independent variable, by considering the characteristic forms of 
the coefficients in a system of Euler-Lagrange equations. Van der Vaart 
(1969) has applied the results of Douglas ( 1941 ) to a system of ordinary 
linear differential equations of second order. A survey of the liters- 
ture on this approach to the inverse problem may be found in the book 

by Funk ( 1970 ). 

The general case of many dependent and independent variables has 
been considered by Vainberg ( 1964 ). By treating the dependent variables 
as points in a coordinate system of functions, the invariance of an action 
integral under the variation of the functions is shown to be analogous 
to the invariance of a potential under the variation of the path of inte- 
gration (Tonti, 1969 a) * This analogy has led to definition of poten- 
tiality conditions for the operators of the system of equations. For 
differential equations, these potentiality conditions are then the necessary 
and sufficient conditions for solutions of the inverse problem of the 
calculus of variations (Tonti, 1969b). 

The treatment to be presented in Sections 2.2 and 2.3 will be con- 
fined to quasilinear differential systems of first and second order, 
respectively. We shall emphasize special forms of the Euler-Lagrange 
differential system, following the approach used by Douglas (l94l), while 
generalizing to the case of partial differential equations. Since we 
are looking for some scheme that generates a Lagrangian density, the 
conditions for the given equations to be Euler-Lagrange equations will 
be established in such a way that, once satisfied by the differential 
equations, an explicit Lagrangian will be derivable. The cases of linear 
and weakly nonlinear differential equations are treated as examples in 
Sections 2.2.3 am* 2.3*3; respectively. 
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In supplement to these conditions, the nonuniqueness in the form of 
the given differential system requires discussion: there are operations, 
such as multiplication of the system by some matrix expressions (Davis, 1929 ) 
and changes of the dependent variables, that can convert apparently non- 
Euler ^Lagrange equations to equivalent Euler-Lagrange equations. In 
Section 2.4, we shall discuss one of these techniques, differential 
transformation of the dependent variables. 4s examples, we derive 
appropriate Lagrangian densities for a resistive transmission line, and 
for a warm collisional plasma in Sections 2.4.2 and 2.4.3, respectively. 

2.2 First Order Differential Equations 

A system of first order quasilinear differential equations has the 
general form, 

cW+D^O (i-1, N; L, M) , (2.1) 

where the coefficients C?^ and are explicit expressions in indepen- 
dent variables x. and dependent variables U^, and is written for 

1 7 i 

CL 

the derivative 3U /3x i . Repeated indices of i (or j,k, ... etc.) are 
summed over the N independent variables; repeated indices of q: (or 
P, Y, ... etc.) are summed over the M dependent variables. We see that, 
in general, the numbers of the cf^ and D a are M 2 N and M, 
respectively . 

The inverse problem of the calculus of variations aims at deriving 
a Lagrangian density, Z [= Z(V®, u a , x^.)], that gives (2.1) as the set 
of Euler-Lagrange equations obtained by extremizing the action integral, 1 



through variations of the dependent variables, U°\ For an arbitrary 
£ U , x^), the Euler-Lagrange equations take the form, 



( t — * • • > Cfc=l, < « t ; M) , 


(2.3) 


where d/dx operates on both the explicit and the implicit x - 

i 

dependences through U?(x ) and U a (x.). 

1 J 1 

Equation (2.3) is a highly specialized form of (2.l): the coefficients 

and D are functions only of the derivatives of £ ( U^, u"^ x ) 

' i ' } i ’ * 

The necessary conditions for (2.1) to be a set of Euler-Lagrange equations 
will be obtained in Section 2.2.1 by elimination of j (u?, U a , x. ) from 
the expressions for and D^. More important, however, are the suffi- 
cient conditions on and D*^ so that a corresponding x (U^ x ) 

i v i’ y i ’ 

exists. In Section 2.2.2, we shall establish conditions which will enable 

us to solve the expressions of and iP for X (U^ TJ 01 x ) 

i ' i’ ’ i ' * 

2.2.1 Necessary Conditions 

In view of (2.3), X must be linear in the U? to give a set of 
Euler-Lagrange equations of the form (2.1). We write 


X=m“u? + H , • (2.4) 

where 7TL [ = ^(U^, Xj ) ] and ft [=> Jt(u a ,x i )] are functions of u a and 

x i • The Euler-Lagrange equations of this Lagrangian density with respect 

Qt 

to U then become 


!?L 3 4 V - at = 

/ 1 Sx i 


(2.5) 
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Comparing (2.1) and (2*5) gives , 


„CO 


au 1 


P 


3U a 


sm 01 

D a = _i 

9x i 


an 


au' 


.a 


( 2 . 6 ) 


The C?^ and D 01 must satisfy the following equations, obtained by 
eliminating the and Tl from (2.6), 


c op + c pa = 0 
i i 


ax i ' 3D? 9 o“ 


ac« atf v acP 

1 1 U 1 

au Y au 01 ai/ 5 


o . (2.7) ' 


2.2.2 Sufficient Conditions 

Given C?^(U^,Xj) and we require conditions sufficient 

to guarantee that (2.6) can be solved for the and Tl in terms of 

the and . Note that the uniqueness of solutions is not required. 

2 

The first equation of (2.6) represents at most M N equations for 
a 6 

the MN unknowns 71^, with U K as independent variables. For mv to 

exist, (2*7) must be satisfied. In particular, C?^ = -C^ a , so that 

the number of distinct equations covered by the first equation of (2.6) 

is reduced to MN(M-l)/2. These can be divided into N independent 

groups of M(M-l)/2 equations for each i. Consider in each of these 

qj 0 v 

groups the subset of equations relating and TIV . if we 

ry 

assume a form for we may then use the first expression in (2.6), 

which yields 



K 


m y 

_ sxn a 

3 ? = 

= an* 3 ’ 

1 ; 

3U a = 

au Y 


(2.8) 
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to solve for JTlj^ 
solutions of yn^ 


and subject to the constraint that the particular 

V 

and JT1J must satisfy 



1 


am? rf 
3? “i? 


( 2 . 9 ) 


This self-consistency condition requires that and 07°" satisfy 

the final relation of (2.7) • 

As a digression, it is of interest to note that this self-consistency 
condition is analogous to the well-known condition that magnetic fields 
are source-free. For the magnetic field, B, and the corresponding 
vector potential, A, we have 


* • B = 0 , 


B = VX A 


( 2 . 10 ) 


The analogy follows by taking a = 1, 0 = 2, Y = 3, and 


B = 


!c 2 3, C 3 1 ,cJ 2 | 

r ’ i 7 i } 



( a 5 a 1 

tdu 1 ’ BU 2 ’ QU 3 ) 

( 2 . 11 ) 


The foregoing argument is valid for any triplet from the set {«} . 

Therefore, with W® arbitrarily chosen, and the final expression of (2.7) 

satisfied, we can always use the first relation of (2.6) to solve for 
^ *Y 

7TL, JTI, etc. Having obtained this set of particular solutions, they can 
be substituted in the second expression of (2.6), together with the 

(X 

given D , to obtain M equations for the single unknown function ^(u 01 ). 

A solution for Jl exists if we have 



( 2 . 12 ) 
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(X 

Use of ( 2 . 6 ) indicates that this condition implies that the D must be 
related according to ( 2 . 7 )- 

The consistency relation of (2.7) for triplets such as 

reduces the number of independent equations in (2.6) from M(M-l)/2 


to (M-l) for each i . 

This follows 

since 

use of 

the equations, 

12 _ 

2 

ml 

r 13 


3JH 3 

i ” 2 

~ ~T > 

1 

= r~3 " 

■ Y > 

3U 

9U 


dir 

dv 


in the final expression of ( 2 . 7 ) will result in 


C 23 < 

1 " S u 3 


m i 2 2 

— i + A(U 2 ,U 3 ) , 

dV 


(2.14) 


where A(t 2 ,U^) is the constant of integration. Since and TR? 

2 3 

also contain arbitrary functions of U and U , because of (2.13 ) , 
it is always possible to adjust them to make £(U^,U^) - 0 . Symbolically, 
we can express the result that (2.14) follows from ( 2 . 13 ) as 


12, 13 - 23, 

12 , 13, l4 - 23, 2k, 3k, 


12, 13 , 14, 1M 23 , 24, . (H-l)M . 

Thus, by imposition of the third condition of ( 2 . 7 ), the first expres- 
sion of (2.6) is reduced to only M-l independent equations for each i. 

1 8 

By choosing TH. arbitrarily, JTV: (p > l) can then be obtained consis- 
tently by use of ( 2 . 6 ). 
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It is now clear that the relations expressed by (2.7) are the necessary 
and sufficient conditions for (2.1) to be a set of Euler-Lagrange equations 
of a Lagrangian density of the form given in (2.4). When expressions 
for CT^ and D satisfying these conditions are given, the procedure 
for deriving the corresponding Lagrangian density will involve only 

direct integrations. It should be noted that since TT1* is arbitrary, 

a 

the set of TIL and the function H are not unique. In addition to the 

05 05 

freedom in choosing 711., the set of JTl. and the function 71 are deter- 

i 7 l 

mined only to within arbitrary functions of the yi. . 


2.2.3 Linear Differential Equations 

As an exercise, we shall derive a Lagrangian density for a system 
of first order linear differential equations for which 


C7 


,ap _op , x 

: r = c r<* i> 


D a . + e“(x t ) 


(2.15) 


The sufficient conditions of (2.7) then reduce to 


ap _ 8a 
°x “ -*i 


dc 


Op 


i_ = d 0^ _ d P0£ 

■v ' 


5x. 


(2.16) 


where the final relation of (2.7) has dropped out. We shall choose 
to be the form, where, by (2.6), the are related according 

to 


c? = m“P - M ? 3 


For convenience, we shall impose the condition - -M^ a . 


n 0L 


With the assumed form of JTL, we obtain by use of (2.6) 



(2.17) 


(2.18) 
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Accordingly, It assumes the form. 


ft = - e¥ + n(* i ) 




»■ 


- d°PU 


( 2 . 19 ) 


Combining these results into ( 2 . 4 ) gives the Lagrangian density 


x = | c^u?tr p + - 2d Qp Ju Q U P - e 0 ^ 


4 * n 


(2.20) 


where n is an arbitrary function of the x . 

2.3 Second Order Differential Equations 

A system of second order quasilinear differential equations can 
always be transformed into a first order differential system by treating 
the derivatives of the dependent variables as new dependent variables 
(Courant and Hilbert, 1966) . However, this approach to the inverse 
problem is unsatisfactory for plasma problems because the introduction 
of new dependent variables is equivalent to introducing artificial 
degrees of freedom. Also, if we leave the number of degrees of 
freedom unchanged, the Lagrangian density for second order differential 
equations describing physical systems generally has a corresponding 
Hamiltonian density, X. Conditions sufficient for the existence of a 
Lagrangian density, £, will then also guarantee the existence of X, 
and allow it to be used in such applications as evaluating nonlinear wave 
coupling coefficients (sturrock, 1960a; Harker, 1970). 

For purposes of example, the system of quasilinear second order 
differential equations will be assumed to have the form 


13 



( i = 1 j . • ■ > 


N; <x,$ - 1, 


( 2 . 21 ) 


Ap! uf + B a = 0 
U ij 


M) , 


where u a is the set of dependent variables; the coefficients A 0 ^ 

ij 

B are functions of U.. U . and x. . and ITT and V a denote the 

i i i ij 

first and second derivatives of U 0 ^ respectively. Since U? . - U U . 




Ji' 


the matrix A7\ can be assumed symmetric in the indices i and j 


without loss of generality. 


2.3.1 Necessary Conditions 

Equation (2.3) ma Y be written 

a2j + + ^ ^L,o 

1J ou^tr^ 1 dU'^9x i oU Q: 

Comparing with (2.21) gives the necessary conditions for that system 
be the set of Euler-Lagrange equations of £ as 


op _ 


2ATT = 


S 2 il 


lj au a su p du tx 3u p 

i J J i 




5 2 £ + 5 2 £ 


bz 


1 su“su p &u?bx ;L au a 

Eliminating £ allows us to express these conditions in the form 

Ya 


hA^ 


- 

ij . Jk 


vP? 

ki 


3 A, 


ki 


K < 

a«“ 

J 

+ 

si?. 

J 

> 



A'P = A^ a 

IJ IJ 

y 

dB a 

*£„ gfe- 
au“ w x 

i 

) J aW 

A°^ 

A ij 

|7i +u 6 i-l 

l s 

.a_l 

SA°^ 
a 0Y + iJ 

aa?T , 
+ ^ _ 1 

/aV + 

i\*\ wJ 


au p J 

iJ su Y 

au a 2 ' 

\oU a 9U Y 

1 1 


aV \ 

iifiivV 

J i 


and 


( 2 . 22 ) 

to 


(2.23) 


(2.24) 


l4 



2.3*2 Sufficient Conditions 




We can assume the general form of the Lagrangian density as 


X « P(ljJ,U P ,x ) + ?J(U P ,x )uj + 0(u a ; x.) 


(2.25) 


where £ is responsible for A?^. Then it follows from the first rela- 
tion of (2.23) that 


2A °e = + 

i j J i 


(2.26) 


where only the set [U?'} are treated as independent variables. By 
interchanging a and p, we see that this differential equation can 
be solved only when = A ji> which is the second relation of (2.24). 

Equation (2.26) constitutes MN(WH-1 ) (N+l )/4 linear differential 
equations to be solved for For p to exist, consistency among the 

third order partial derivatives is required, i.e. 



(2.27) 


This implies that the first condition of (2.24) must be satisfied. 

In the cases where either a = p or i = j, p can be obtained 
by simple integration of (2.26). When both a £ p and i ^ j, (2.26) 
becomes an ultrahyperbolic differential equation with constant coeffi- 
cients (Koshlyakov, Smirnov, and Gliner, 1964). Its particular solutions 
can be obtained by Fourier transforms (Koshlyakov, Smirnov, and Gliner, 
(1964). The general solution of Q should include that of the homo- 
geneous equation. 
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0 


(2.28) 


- 

ai/^au^ 

i j 


il 

a °>? 


which can be solved by separation of variables (see Appendix A). We shall 
show below that this nonuniqueness in p will be restricted by other 
sufficient conditions. 

From the second expression of (2.23) and (2.25)> we obtain. 


\SU H 31 W " °~i 3U U x au^3U p 3U^3x i dU^ 


Pi 


(2.29) 


Since (2.29) is an identity in U^, \P, and x_^, the right-hand side 

OL Ot 

must be linear in IL . Taking the derivative with respect to th, yields 


a?“ a? p 


li - _i . / 1- + u Y 3_\ *l2 + ^ 

3U^ 3U a \ J 3U Y / 3 Uj 3U^ 3U a 3U^ 9%yiP 3U? 


(2.30) 


the left-hand side of which is antisymmetric in ot and p* The third 
expression in (2.24) guarantees that the right-hand side is also anti- 
symmetric in a and p. However, the fourth expression, which shows that 
(2.3l) is symmetric in i and j, is not sufficient to establish that 
the right-hand side of (2*30) is independent of U^. This requirement 
can be obtained by differentiating ( 2 . 30 ) with respect to 


r/|_ + u* *_ 


a "1 

[W k go« 

1 * * 

au p . 


.1 ^ 

J aujau] 


+ 




^3 p 

rY au Y 3XJ a 3U^ 3U a 3U?3U Y 3l/?3U Y 


. (2.31) 


1 J 


1 J 


This is more restrictive than the third necessary condition, and represents 
limitations on p in addition to the fourth. 



Additional conditions from (2,30) for to exist are 


a$^ Y as Ya 
au Y au a &P 


( 2 . 32 ) 


ana 


logous to the last expression of (2.8), where is defined by 


a 0 * 3 . 

X 


ZB 

au? 


l! - /&_ + u Y $-2— 

f ys'j J auV au“a«f 


(2.33) 


a 

The solution of (2.30) enjoys one arbitrary choice of just as the 

rv 

solutions of the first expression of (2,6) for for each i. But 

if ( 2 . 29 ) is to give a consistent solution for Q, additional conditions 
restricting the set of then follow from ( 2 . 29 % from the requirement 

that 



(2.3*0 


analogous to (2.12). We have 


/s_ , „y a_\/£s .a g iL.-\ + 

\^*i 1 SUV \au a bU^ QU P SU “/ csU^ SU a 



In summary, sufficient conditions for (2. 21) to represent a set of 
Euler-Lagrange equations of a Lagrangian density, £, of the form (2.25) 
are as follows: for a solution of ( 2 . 26 ) for $ to exist, A^ must 
satisfy the first two conditions of (2.24); for a solution of ( 2 . 30 ) for 
5\ to exist, p must be further restricted by ( 2 . 31 ) an ^ ( 2 . 32 )> while 
A?^ and Ep must satisfy the third condition of (2.24); for a solution 
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Q Qf 

of (2.29) for Q to exists g^* and B must be restricted by (2.35)- 
These conditions* as well as those of Section 2.2.2 for first order equa- 
tions* agree with the general forms obtained by Tonti ( 1969 b) by poten- 
tiality analysis of differential operators in function space (Vainberg* 
1964; Tonti* 1969 a). 

2.3.3 Nonlinear Differential Equations of the Second Hank 

As a demonstration of the foregoing results* sufficient conditions* 
and the corresponding Lagrangian density* will be obtained for the 
differential system 


ij 


ij 


+ a^Y 


♦ u k Y 


B a , 


b“ + + b<»W + bW + b°PW + b°P w 

11 11 ij 1 J 


(2.36) 


The coefficients, etc. and b a , etc., are functions of ; the 


.00 


a £7 are symmetric in i and j, the b^Y in p and y, and the 

b?\^ in (p,i) and (y,j). With the coefficients of (2.36), the equa- 
«J 

tions of ( 2 . 2 l) are nonlinear and of the second rank because of the 

presence of &? , a** b°** b ^Y and 

ij ijk ; * 1 ij 

Substituting ( 2 . 36 ) into (2.24) yields for the a 0 ^, etc., 

“ 1 


a°P . aP a , a 0 * 37 , a *** 7 , a 0 * 37 _ a P«V 

aj ij ij ij ’ ijk ijk ' 


a 0 * 37 + , a ? 70 + ,W 

ijk jki kij kij 


(2.37) 


and for the b a , etc., 
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>f + = iw 1 ) * + f 


Za°P y 

h cpy + b w _ a osv + lijk 

iJ ij ij 5x k 


We can now assume the general form for ^ as 


' 9 , S V 


3 - CW + ?( a S + a S v ) ,, ? D ? > 


and choose to impose 


^QPY _ P«V _ QYP 
ijk ^jik ~ ^ikj 


The requirement of (2.31) reduces to 


- *(* 2 ? ♦ -sr - ■£?) . 


while (2*32) becomes 


r - >r a * »r* - tr( a S Y + c + ‘if) ■ <*•*> 


Using (2.30) we find that p® takes the form 


= pfl^ + pfW , 


: 2 A 3 ) 


where and are to be obtained from the differential relations, 

t -v oft cfty 

pf - pf = tf - ^ , pfv . . |ii.j . (2-4 , 
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The conditions of (2.35) reduce to 


b°P - 




b OpY _ b pO^Y = 1 a_[ b 0f3Y _ 


2 | i 


1J 

SX 3 l 


(2.45) 


without further restricting the choices of p?^ and p?^\ 
The general form of G is 


Q = -bV + q ^xf^P + q^V^U 7 + q(x. ) , 


(2.46) 


where, according to ( 2 . 29 )> 






,c£Y 


/ atl c«pYv 

= - i(^ Y - It) - < pcff - <*■ 


47 ) 


and q(x^) Is arbitrary. 

We find that the f)?^ are determined from (2.4l) only to within an 

arbitrary curl tensor, e , where e. is the antisymmetric 

kmn ljn m kum 

unit tensor* According to (2.44) the p^ and p?^ are hoth deter- 
mined only to within an arbitrary tensor symmetric in a and (3, 

In the special case that the set of differential equations expressed 
by (2-21) are linear, the relevant sufficient conditions reduce to 



(2.48) 


These are equivalent to those used by van der Vaart ( 1967 ); who con- 
sidered the case N = 1 [Equation ( 16 ) by van der Vaart ( 1967 )]. 
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2.4 Differential Transformations 


In Sections 2.2 and 2.3, we assumed that the oth equation of the 
differential systems represented by (2.1) and (2.2!) is an Euler-Lagrange 
equation with respect to variation of U^. Now these sets of equations 
are transformable, for example by change of variables, to equivalent 
sets which no longer satisfy the sufficiency conditions. Solution of 
the inverse problem of the calculus of variations is consequently less 
restrictive than is suggested by the sufficiency conditions: it may be 
possible to convert seemingly non-Euler-Lagrange equations to Euler- 
Lagrange equation form by use of appropriate transformations. The 
problem this poses is how to recognize when such transformation is 
possible. Here we shall comment briefly on the transformations likely 
to be involved. We cannot provide a general solution. 

Among the range of possible transformations, those retaining the 
number of dependent variables unchanged include: (a) matrix transformation 
of the dependent variables; (b) matrix transformation of the differential 
equations, by use of integration factors, leaving the dependent variables 
unchanged (Davis, I 929 ) J and ( c ) differential transformation that raises 

the order of the differential equations. For a self-consistent system 

2 

of M differential equations^ Method (a) involves M functions of 

2 

as the elements of the transformation matrix; Method (b) involves M 
functions of x_^ as matrix elements, the dependent variables, and 
perhaps their first derivatives, while Method (c) involves only M 
functions of x_^, the new dependent variables, and their derivatives. 

In each of these cases, the number of functions at our disposal in 
general falls short of the number of sufficient conditions to be applied. 
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Consequent ly, we cannot expect that these methods will always be success- 
ful in transforming non-Euler-Lagrange equations to Euler-Lagrange equations. 

Despite their inability to guarantee successful transformation, 
methods (a)-(c) are of interest in dealing with equations of forms, such 
as those occurring in some physical and engineering problems, where the 
number of sufficient conditions is reduced. Well-known examples can be 
readily found in the cases of replacing electromagnetic field variables 
by potential variables (Goldstein, 195O), and velocity variables by 
Clebsch variables (Lamb, 1930), before the Lagrangian densities can be 
obtained. To illustrate their use here, we shall apply a linear differen- 
tial transformation to convert a first order linear non-Euler-Lagrange 
system to a second order Euler-Lagrange system in Section 2.4.1. Using 
the sufficiency conditions obtained in Section 2.3*3> we shall then 
determine the maximum number of dependent and independent variables for 
a successful transformation to be possible. As examples, a resistive 
transmission line, and a warm collisional plasma, will be studied in 
Sections 2.4.2 and 2.4.3, respectively. 

2.4.1 Differential Transformation of Linear Differential Equations 

Under the linear transformation, 

U°: _ t°^(Xj)v£ + s Q P( Xi )v® + R o: (x i ) , (2.49) 

the first order linear differential equations of (2.1) and (2.15) 
become second order equations, with coefficients of the form (2.36) 
given by 
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7 


a°e = i/c^T YP + dV + c 0 ^ 

ij 2 \ i j J i / ’ i j\8Xj/ i i 

b oe = cf (^y d 0 ^ , b“ = cf(^y d^ + e“ . ( 2 . 50 ) 

The sufficiency conditions imposed on these coefficients by (2*4-8) 
constitute MN(M-l) (JW*l)/4+ 1^(1^ l)/2+ M(M-l)/2 conditions to be 
satisfied by the M^(Nfl) unknown functions, and S^. The 

number of unknowns will consequently be no less than that of the 
conditions only when 

m(m-i)(n+i)(n- 2 ) < 4m . ( 2 . 51 ) 

When (2.5l) is satisfied, the equations obtained by substituting 

Ct£3 Q Qnj 

(2.50) into (2.48) are in general solvable for the T. p , S , and R , 

aB ccb a 

for the given coefficients , d , and e . The resulting second 

order differential equations will automatically become Euler- Lagrange 
equations, with the corresponding Lagrangian density derivable using 
the results of Section 2.3.3* 


2.4.2 Resistive Transmission Line 

The first order equations for a linear resistive transmission line, 
as shown in figure 2*1, take the form, 


31 

ax 


+ 




+ ^ + „ « 0 

ax 


(2.52) 


where I, V, C, L, and R denote the normalized current, voltage, 
distributed capacitance, inductance, and resistance, respectively. 
According to (2.1 6 ), (2.52) is not in Euler-Lagrange equation form. 
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Figure 2.1 Resistive Transmission Line 



[f,g] - {v\v 2 }, 


the procedure outlined in 


With {l,V} = [XI 1 , U 2 } and 
Section 2.4.1 may be followed to give 


C + c ' tt 

3x at 


|- + LC |r 

dx &t 


+ LC |r - RC 

ax at 

L — + L It ■ R 

ax at 


LJ 


(2.53) 


as an appropriate transformation. The corresponding Lagrangian density 
is obtained, through the use of ( 2 . 25 )> (2.39)> (2.43)> (2.44), (2.4 6), 
and (2.4-8), in the form. 


*-5 


4- 

f + 2 a£ 

U + L (| 

\ 2 ~ 

-) 

+ C + L If) 

\3x 

/ ax 

ax \a 

x / _ 

o x \3t at / 


L M 

tc* 

+ ai\ + b£ 

c(||f 

ax 

\ at 

at/ 

2 

\s t / 

RC 


(M + 

a£ \ _ Ju 

+ M\ 

2 

& 

\ax 

at / 

\ax 

at/ 


+ 2 — ' — + 
at at 


C ? 2 

' 2 RB 




(2.5*0 


The Euler-Lagrange equations of (2.54) can be shown to agree with the 

result of using (2.53) in (2.52). Because the number of sufficient con- 

ditions is less than the number of transformation coefficient s, 

and y ( 2 * 54 ) is but one example among an unlimited number of appro- 

0 

priate Lagrangian densities, 

} 

2,4 -. 3 Warm Collisional Plasma 

For small one-dimensional perturbations in a plasma with a homogeneous 
immobile neutralizing positive ion background, the macroscopic equations 
can be written as. 
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y 


y 


5 E ] en 

* x e 0 


= 0 


an 

at 


+ n 


1 ^ = 0 


0 d x 



+ mvv + VPq 



+ eE - 0 


(2.55) 

where m and e are the electron mass and charge; and P Q are the 

quiescent electron density and pressure; v is the effective electron- 
ion and electron-neutral momentum transfer collision frequency; y is the 
adiabatic index, and e Q is the permittivity of free space. The pressure 
term, \P q -|^ n^, results from assuming an adiabatic equation of state 
for the electrons, 


( p o +p ) (n o +n)_Y = p o n o Y • (2-56) 

In (2.55)> the dependent variables are the electric field, E, the 
electron density perturbation, n, and the electron drift velocity, v. 
Application of (2.l6) indicates that these equations are not in Euler- 
Lagrange equation form. They satisfy ( 2 . 51 ), however, so that the 
differential transformation defined in Section 2A.1 can be used. To 
reduce apparent complexity, it is helpful to rewrite ( 2 . 55 ) in terms of 
normalized variables, defined by 


E = 


eE 

mv cc 

t p 


n 


n 




X 


CJD X 

P 


T 7 = CD t 

p 


V 



(2.57) 


where [ = (Pq/hi) and ^ [ = (n^e^/me^)' 5 ''^] are the electron 

thermal velocity and plasma frequency, respectively. We then have, in 
place of ( 2 . 55 ), 
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From (2.48) - (2.50), and ( 2 . 58 ); it can be shown that are appropriate 
transformat ion , 
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should have the following elements: 
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°33 = ' h + 1 ’ 


(2.60) 


where a, b, c , d, and e are arbitrary constants, except that they should 
make (2.59) a reversible transformation. 

The corresponding Lagrangian density, obtained through using (2.25); 
(2.39)* (2.4-3); (2.46), and (2.48), may be expressed in terms of f, g, 
and h as, 
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( 2 . 61 ) 


The Euler-Lagrange equations of (2.6l) can be shown to agree with the 
result of using ( 2 . 59 ) and ( 2 . 60 ) in ( 2 . 58 ). 

2-5 Discussion 

In this section, we have considered the inverse problem of the cal- 
culus of variations for systems of first and second order quasilinear 
partial differential equations. The approach has been to compare the 
form of the Euler-Lagrange equations with that of an arbitrarily chosen 
set of equations. The resulting sufficiency conditions agree with 
general results obtained previously by the more abstract method of imposing 
conditions of potentiality of operators in a function space (Vainberg, 

19 64; Tonti, 1969)4 By restricting ourselves to equations of quasilinear 
form, we were able to determine the Lagrangian density explicitly, if the 
sufficient conditions are satisfied by the given set of differential 
equations. As examples, the results were applied to systems of first 
order linear equations, and second order nonlinear equations. 

The explicit formulation described here has led to some success in 
using differential transformations to convert non-Euler-Lagrange equations 
to Euler-Lagrange equation form by changing the dependent variables. The 
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examples on the resistive transmission line and warm collisional plasma 
served to show their importance. This operation, together with the 
other transformation techniques discussed in Section 2 , 4 , constitute 
powerful mathematical ways to define appropriate generalized (dependent) 
variables in a given inverse problem. This is complementary to the 
procedure of choosing suitable 'geometric’ and ’force 1 variables in 
physical variational problems, described by Penfield and Haus (1967)- 

Our discussion of transformation techniques suggests some significant 
problems for further study in the inverse problem of the calculus of 
variations. For instance, recognition of equivalent systems of formally 
different differential equations becomes important if the area of appli- 
cation of the inverse problem is to be enlarged. 

Although the mathematical approach is rigorous, the need to choose 
appropriate generalized variables, and to rewrite the equations in Euler- 
Lagrange equation form, reduces the practical value of our approach to 
the inverse problem in those physical situations in which intuitive energy 
considerations can be readily used to guess at, and establish appropriate 
Lagrangian densities. The latter approach relies on the experience that 
a physical system with well-defined energy can be described in terms 
of Hamilton’s variational principle. The trial procedure then involves 
assigning generalized variables, according to the number of degrees of 
freedom, before guessing a Lagrangian density that includes various 
forms of kinetic and potential energies. The trial Lagrangian is then 
checked by the relatively simple process of applying Hamilton’s principle 
to produce Euler-Lagrange equations that are the correct equations of the 
problem. A small number of trial-and-error exercises may then result in 
a suitable Lagrangian density. 
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In anticipation of the Lagrangian applications to plasma physics, 
to be discussed in Section 3 > we may point out two other problems* The 
first concerns the explicit distinction between dependent and indepen- 
dent variables generally assumed in the inverse problem; the variation 
principle established in Section 3 relies on an unorthodox arrangement 
which we term the T dual role' of the variables* This dual role can be 
understood to be some implicit dependence, in the form x[y(z)], with 
both x and y assuming the role of generalized variables. The exten- 
sion of the present results to incorporate this dual role is a problem 
that we have not had time to study in detail. 

The second problem is related to the so-called ’Lagrange (or Bolza) 
problem* (Rund, 19 66; Bliss, 1946). This deals with the variation 
principles whose solutions are to be determined under subsidiary constraints. 
This type of problem immediately introduces serious practical difficulties 
into the corresponding inverse problem, since there is no mathematical 
rule to determine which of the given differential equations are to be 
treated as constraint equations. In the corresponding physical problem, 
however, no such difficulty arises once the degrees of freedom are 
determined , 
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3. LAGRANGIAN AND HAMILTONIAN DENSITIES FOR PLASMAS 


3.1 Introduction 

In this section, appropriate Lagrangian and Hamiltonian densities 
will be established for the macroscopic model, in which the plasma is 

described by Maxwell's equations and moments of the Boltzmann equation, 

\ 

We shall first construct the Lagrangian density from energy considerations, 
rather than the formal mathematical approach of Section 2, and then 
verify its validity by applying Hamilton's principle to obtain the re- 
quired Euler-Lagrange equations. Next, the corresponding Hamiltonian 
density will be derived, 

Lagrangian densities have already been obtained for plasmas under 
various assumptions. In the microscopic model, in which the plasma is 
described by Maxwell's equations and the Vlasov equation, an. appropriate 
Lagrangian has been obtained by Low ( 1958) * The Lagrangian in a 
relativistically covariant formulation was discussed by Sturrock ( 1958 a), 
who pointed out a difficulty associated with the choice of variables: 
the calculus of variations distinguishes between the generalized (dependent) 
and the integration (independent) variables; the description of plasmas, 
however, involves the charged particle trajectory, which is conventionally 
treated in Lagrangian coordinates, and the electromagnetic field, which 
is conventionally treated in Eulerian coordinates. The distinction is 
illustrated in figure 3 * 1 ; which shows a particle trajectory in phase 
space. The Eulerian coordinates, (x, v,t), describe the particle with 
relation to a fixed set of axes, whereas in the Lagrangian description, 
the axes follow the trajectory and the particle location is denoted by 



ca 


Figure 3.1 In phase space, a particle trajectory ( ) 

be specified by its Lagrangian (initial) coordi 
nates, (x^v^O), or its Eulerian (present) 
coordinates, (x,v,t). The polarization vector, 
1* used b y Sturrock (1958), connects the real 
particle trajectory ( ) to a specified tra- 
jectory ( ). 




To resolve this difficulty, Sturrock introduced a field-like” 
expansion of the total Lagrangian with the particle position vector 
written as x + £, where g; is the polarization vector with respect to 
a specified particle position x. The case when f is defined by com' 
paring the two position vectors simultaneously is also shown in figure 3*1* 
Thus x becomes the integration variable for both the particle displace- 
ment, g , and the electromagnetic field. On the other hand. Low pre- 
ferred to modify the conventional formulation of the calculus of variations 
in order to incorporate the two types of variables into his Lagrangian; 
it is composed of two parts, one for the motion of the particles in a 
given electromagnetic field, and the other for Maxwell’s equations with 
given particle trajectories. In the application of Hamilton’s principle, 
this requires that the electromagnetic potentials be treated as given 
functions of particle spatial position, x , in the variations with 
respect to the particle trajectory, x (x^, v , t), and that x be 
treated as an integration variable in the variations with respect to the 
electromagnetic potentials. Since the particle spatial coordinates play 
the dual role of generalized and integration variables, while the 
electromagnetic potentials play the dual role of given functions and 
generalized variables, in what follows we shall describe this as the 
"dual role" approach, 

A Lagrangian analogous to that of Low ( 1958) and Sturrock ( 1958 a) has 

not yet been presented for the macroscopic plasma model. There have 

< 

been successful attempts at derivation for simplified models, often using 
Lagrange multipliers that necessitate the use of generalized variables 
not corresponding to the physical degrees of freedom. For example, 
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Katz (1961) has given an appropriate Lagrangian for a one-fluid, com- 
pressible plasma, with scalar pressure and adiabatic motions assumed. 

He obtained Maxwell’s equations and the equation of fluid motion by 
application of Hamilton’s principle, with the use of two Lagrange multi- 
pliers to incorporate the fluid mass and charge continuity conditions. Su 
(1961) has obtained a suitable Lagrangian for an inviscid, adiabatic plasma 
with infinite electrical conductivity. The resulting Euler- Lagrange 
equations include Maxwell’s equations, number density continuity equa- 
tions, and three other equations from which the correct force law can be 
derived. Because he used some implicit form for the pressure term in 
his Lagrangian, the number density, n, was also treated as a generalized 
variable. Newcomb (1962) derived the Lagrangian for a hydromagnet ic 
plasma, without using Lagrange multipliers, for the cases of scalar and 
axisymmetric pressure. The mass and magnetic flux continuity equations, 
and the plasma adiabatic pressure state equation, were treated as subsi- 
diary constraints that define the variation of the Lagrangian and result 
in the correct force law. 

The foregoing results (Katz, 196 1 ; Su, 1961; Newcomb, 1962) are 
characterized by the absence of the dual role of the variables used in 
their formulations. In contrast, our own approach in Section 3.2 will 
follow that of Low closely: the plasma cell trajectories and the electro- 
magnetic potentials will be used in dual roles. The Lagrangian thus 
established includes tensor pressure and elastic collisions. Following 
the technique of Newcomb, the first (continuity) and third (heat flow) 
moment equations are used as subsidiary constraints to show that the 
Euler-Lagrange equations obtained by application of Hamilton’s principle 
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are the second (momentum) moment equation and Maxwell’s equations. Our 
reasons for taking the "dual role" Low approach, rather than the polari- 
zation vector approach of Sturrock, will become clear in Section 4.2, 
where we shall show that the perturbation expansion necessary to treat 
nonlinear wave-wave interactions is thus greatly facilitated. 

After the appropriate macroscopic Lagrangian is established through 
the dual role approach, an important question arises: does this form of 
Hamilton's principle, synthesized at the expense of the explicit distinc- 
tion between the dependent and independent variables, form an acceptable 
basis for a plasma canonical formulation that agrees with the known 
plasma equations? If this can be answered affirmatively, the macroscopic, 
as well as the Low Lagrangian, will be acceptable in their total form. 

This constitutes our motivation to derive the appropriate plasma Hamiltonian 
and canonical mechanics through the dual role approach. Furthermore, the 
corresponding Hamiltonians will be readily subjected to perturbation 
approximations, analogous to those applied to the Lagrangians (Low, l 958; 
Sturrock, 1958a; Newcomb, 1962; Galloway and Kim, I971 ; see also Section 4). 
The resulting approximate Hamiltonians, although not to be derived here, 
should have important applications to nonlinear plasma problems (Sturrock, 
1960a; Harris, 1969 ; Harker, 1970) . 

Sturrock (1958a) and Newcomb (1962) have derived Hamiltonians for 
their quadratic plasma Lagrangians. In classical mechanics, if a 
given .Lagrangian is at least quadratic in the derivatives of the 
generalized variables, a corresponding Hamiltonian can always be 
obtained (Goldstein, I95O) * Obvious difficulties arise in the case of 
the macroscopic Lagrangian to be derived, in which the variables play 


35 



dual roles* Sections 3*3-2 to 3.3.^ will show how to modify the well-known 
Legendre transformation and Poisson brackets so that the dual role is 
embodied in a self-consistent Hamiltonian formulation. For completeness, 
a Hamiltonian is derived in Section 3-3*6 from the Low Lagrangian to 
facilitate a comparison between the microscopic and macroscopic models 
of plasma Hamiltonian mechanics. 


3*2 Lagrangian Density 

Newcomb ( 1962 ) has given the following Lagrangian density, 

= ‘ T N ” ^ = 2 pU ' U N = y-1 + 2jl^ , (3-1) 

where 3^ and denote the kinetic and potential energy densities, 

respectively; p is the mass density; u and P are the plasma velocity 
and scalar pressure; B is the magnetic field; y is the adiabatic index, 
and is the permeability of free space* Gravity has been neglected, 

and the hydromagnetic assumption of infinite conductivity, 

E + u X B = 0 , (3*2) 


has been made. 

Choosing x as the generalized variable, he applies Hamilton's 
principle to the total Lagrangian, 



V 


(3.3) 


The variations in p, p, and B, due to that in x, are defined by the 
constraint equations of mass continuity, adiabatic state, and magnetic 
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flux conservation, respectively, 

- y 

Pdx = constant , Pp = constant , B . da = constant , (3.1+) 

where da denotes an area element moving with the fluid* The resulting 
Euler-Lagrange equation is the appropriate force law in the hydromagnetic 
approximation, 


( 3 . 5 ) 

With the general form of (3*l) in tnind, the main part of £ will 
be obtained in Section 3.2.1 by integrating the Low Lagrangian in velocity 
space* Use of Eulerian coordinates in Section 3,2.2 enables us to 
establish the necessary corrections to complete this approximate Lagrangian. 
Section 3.2.3 confirms its validity by demonstrating that the corresponding 
Euler-Lagrange equations are Maxwell's equations and the force law for 
macroscopic plasmas* 

3.2.1 Macroscopic Approximation to the Low Lagrangi an 

The total Low Lagrangian, L > has the form 

L 

\ "Sj" a .lqf d , v _o t( $y "a* 4 + f d i:(~§ - fe) ’ 

V V ' 

0 v 

tm I v2 (^0- > () - Ot ‘ ) - vtXg, Vg, t) ■ A(*,t)J , (3.6) 

where the summation £ implies more than one particle species; f is 
the velocity distribution function; x and v are the particle position 


it <& + 


V ♦ (puu) + VP - (7 X B) X B = 0 
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and velocity at time t; and v^ are their values at t = 0 ; 

V denotes the spatial volume occupied by all the particles under con- 
sideration at time t; Vq denotes the spatial volume occupied by the 
particles of a species at t = 0 , as is schematically shown in figure 3 . 2 ; 
m and q denote the particle mass and charge; is the permittivity 

of free space; E, B, to, and A represent the electric and magnetic 
fields, and the scalar and vector potentials, respectively, where 


E = -Vcp - 


CSA 


B = V X A 


(3-7) 


Note that we have chosen to use V and V , as shown in figure 3 . 2 , to 

confine ourselves to a plasma temporarily enclosed in volume V at time 

t, as opposed to letting V =» «> > as in Low's and Newcomb's work. The 

double integral contains the part , the Lagrangian due to each particle, 

whose first term represents kinetic energy and second term represents the 

interaction energy of each particle with the electromagnetic field. The 

second integral describes the energy associated with the free space 

electromagnetic field. A careful discussion of the application of 

Hamilton's principle to has been given by Galloway and Kim (1971). 

The variation of L due to the variation in x (x^, v , t), constrained 

- 1 - 1 ^ ~Cr 

by the continuity of particles in phase space, 

f(x,v,t) dx dv - constant , ( 3 - 8 ) 

and with the electromagnetic field treated as a given function of x, 
gives the particle force law. The variations in cp and A, with x 
treated as the integration variable, yield Maxwell's equations. 
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Figure 3*2 Particles of a species that occupy a spatial 
volume V at time t are assumed to occupy 
a volume at t = 0. This approximation 

is acceptable when t is sufficiently close 
to t - 0* 


) 
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In anticipation of Section 3.2.2, the integration variables in L 

L 

may be transformed from the Lagrangian (initial) phase space (x , v ) 
to the Eulerian phase space (x,v). Applying the conservation of particles 
in phase space 


f <3o’ V “io d i - f ( x ’ v ’ t ) dx dv > 




(3-9) 


transforms ( 3 * 2 ) to 


l l= f dxJ l 
V 


h. -E L f (i'Z' t ){l 1,2 - w + Ti) + 


e o E 


b; 

2 H, 


0 


(3-10) 


Equation ( 3 . 10) applies to a collisionless (Vlasov) plasma, for which 
f is the smoothed velocity distribution function. The elementary volume, 
or cell, dx, must be much larger than the mean particle spacing to 
justify the smoothing, and much smaller than the mean free path to justify 
neglect of collisions. The electromagnetic field variables are treated 
as constants within dx, and represent the collective charge effects. 

We shall refer to this model, used by Low, as the "collisionless micro- 
scopic model” . 

The Lagrangian, L^, for plasmas in the ” collisional microscopic 
model 11 , is identical to that of ( 3 * 10 ), but with different interpretations 
for f, dx, E, B, cp, and A: the distribution function, f, is now fine 

'V /n; r*,/ 

grained; the cell, dx, roughly speaking, becomes smaller than the mean 
particle spacing and larger than the size of the particles; the electro- 
magnetic field variables include the microscopic particle self-fields 
effective in collisions. 
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Despite their similarity in form, and are not equal in 

value. The difference can be determined by relating the quantities in 
the two models by 


CP= ' P M +CP R ' + *R ' etC - 

f(x,v,t) = f (x,v,t) + f (x,v,t) , ( 3 . 11 ) 

rvi 1’1 /V XV 


where subscript R denotes quantities which average to zero in a cell 

of the collisionless microscopic model. Substituting (3.11) in ( 3 . 10) 

yields the macroscopic approximations, L' and L ', of L and L , as 

T L T L 



where Tr represents the trace of a tensor; the macroscopic density, n, 
drift velocity, v , and pressure tensor, P, are defined by 



and the size of dx corresponds to that of the collisionless microscopic 
model . 


y 

/ 


kl 



If Hamilton 1 s principle is applied to the correct macroscopic 

force law cannot be obtained as one of the Euler-Lagrange equations* The 
discrepancies arise in two ways* First, there is no term in that 

accounts for the energy associated with the heat flux, q, defined by. 


q = 


m 

2 




x,v. 


t) 


(3.14) 


To perform the variation of the term Tr p/2, one then has to assume 
q negligible. Second, the Euler-L a grange equation resulting from varia- 
in x (x , t) is found to be correct only if the sign preceding the term 
Tr P/2 in is reversed. Clearly, it would be convenient if we 

could prove that £" -2 Tr P. This would involve a statistical theory 

of plasma fluctuations (Harker and Crawford, 1973)* which may be con- 
sidered as elastic multiparticle collisions, with the added complication 
of elastic two-particle collisions. We shall take a much more direct 
approach, however. Since we already know that the corrections should 
contain -E Tr P, and may result from including energy contributions 
from the plasma fluctuations and particle collisions, we shall obtain 
them from energy considerations, 

3*2.2 Macroscopic Potential, ]} 

We postulate that, for each particle species, the corrections can 
be considered as a macroscopic potential energy density. If, that has 
two parts: a random (collisional and f luctuational) potential, V , 

(when particle collisions are sufficiently frequent) and a heat transport 


potential, 1? , 

H 


defined by, 





l R " AV f d * 

AV 


dv f(x,v,t)F. 


* .x_ 


V 


V*Q 


( 3 . 15 ) 


where, as shown in figure 3*3; AV is a macroscopic volume with size 
much larger than the particle mean free path and a boundary moving with 
a fixed number of particles of a species; the random collisional force, 
F p , and position of the particles, x^, are defined with respect to the 




macroscopic force, F M , and the position, R, of AV> and Q is the 
total density of heat energy transported across a macroscopic cell as it 
moves along its trajectory, x (x ,t) , 


Q = j dt q[x(x Q ,t),t] . 

/•V./ 1 I 'N-i' 


(3.16) 


To obtain an explicit macroscopic expression for XL, consider the 

K 

following equation which appears in proving the virial theorem in classi- 
cal mechanics (Goldstein, 1950b), 


— I dxjdv f(x,v,t)mv • x - j dxldv f ( x, v, t ) [mv^ 
AV AV 






X] , (3.17) 


where we have used the fact that the number of particles of a species in 
AV is constant. Rather than averaging over a long time period, as in 
classical mechanics, we shall consider the approximation over a large 
number of particles within AV. Expressing (3.17) in terms of the macro- 
scopic and random quantities yields 


d 

dt 


J 

AV 


dv f m(v 
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Figure 3.3 The size of AV is assumed to be much larger 
than the particle mean free path so that the 
same particles of a species are contained in 
AV for a time period much larger than the 
mean particle transit time across AV. This 
assumption of large Av size enables us to 
justify the necessary correction to the 
v- integrated Low Lagrangian. It will be seen 
in Section 3*5 that the large size assumption 
is not necessary for the Lagrangian of (3.22) 


Since 


where v is the random particle velocity with respect to v i 
v = dR/dt and F = mdv /dt, (3.18) simplifies to 

r-vJVl 


I dxldv f m v_ . x. = AV Tr P + f dxldv f F • x . ( 3 . 

dt J ^ ^R ^R ~ J ~J ~ 

AV AV 

The left-hand side of (3.I9) involves d(AV)/dt. In the case when this 
is zero, e.g., in a stationary plasma, combination of (3.16) and (3*13) 
then gives the expected expression, 


IL = Tr P . ' ( 3 . 20 ) 

R ^ 

This agrees with a result obtained by Gartenhaus (1961), who showed that 
the collisional interaction potential of a stationary gas is proportional 
to the particle thermal energy. 

Combination of (3.1 6) and ( 3 . 20 ) gives a guess for at least 

for the case when sufficiently frequent collisions can be assumed, 


Z" . - 21 ) , U r Tr P + V*Q 


(3.21) 


Substitution of ( 3 * 21 ) in ( 3 * 12 ) then gives the total Lagrangian as 


L = J dx Z , 

V 

where the subscript M will be dropped from here on. 


| Tr 


P - V-Q - qn(cp M - 
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( 3 - 22 ) 


Firsts we should point out that the assumption of a stationary 
plasma in obtaining (3-20) can be relaxed for the purpose of obtaining 



L in (3.22), if we let V tend to infinity. This can be seen by making 
AV -* dx in ( 3 . 19 ), and integrating dx over V -* « } making the left- 
hand side automatically zero because no particle now crosses the boundary 
of V. Second, the assumption of sufficiently frequent collisions requires 
that the size of AV be much larger than the particle mean free path. In 
the following, it will be seen that this requirement is overly restric- 
tive since the macroscopic quantities in ( 3 * 22 ) are well-defined in a 
cell whose size need only be much larger than the mean particle spacing. 
Further discussion of this point will be included in Section 3.4. 


3.2.3 Application of Hamilton’s Principle 

We shall now test the validity of (3.22) by application of Hamilton’s 
principle. Note that L contains the generalized variables x (x n ,t), 

A (x,t), and cp (x,t). 

Variation with respect to A, and cp: The variation of L with 
respect to A and cp is comparatively simple since x is now treated as 
an integration variable. For variations 6A and 59 , which vanish on 
the boundary of V and at times t^ and the Euler-Lagrange equa- 

tions corresponding to 

( 2 

I = J dt L , (3.23) 

h 

take the standard form (gchiff, 1968 ), 
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(3.24) 
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where the functional derivatives of L with respect to A and cp, in 
Eulerian coordinates, are defined by 


gL _ &X 

a [" ax 

gL _ gZ 

a 

ax 

gA SA “ 

r+-/ 

ax i [a(5A/ax i )_ 

* ^ a?' 

ax. 

_a(a«p/ax i )_ 


Using (3*7) in (3*22), and then applying (3*24), leads to the Euler-Lagrange 
equations, 

1 5>E 

— VXB^e^^+2 qnv , S V . E = E qn . (3-26) 

|J.q rv vJ 0^” U 

These equations, together with (3«7 )j form the set of Maxwell equations. 

Variation with respect to For variation, the Lagrangian is 

simply. 


E f E 
L = 1 dx X , 

V 

Z = I nV D 2 ' I Tr Z " 7 ’,S, -qn ( Cp “^D * > (3-27) 

where the superscript E signifies Eulerian coordinates, and the integral 
of the elect romagnetic field energy in (3*22) has been dropped because it 
has no effect on the resulting Euler-Lagrange equation (Hill, 1951 5 
Galloway and Kim, 1970* The Quantities ^v^, n, Tr P, and V«Q 
are to be considered as implicit functions of x via subsidiary con- 
straints. [Otherwise, the resulting Euler-Lagrange equation will be 
a zero identity (Hill, 195 ^) •] These constraints are the first and 
third moment equations derived from the Boltzmann equation (Braginskii, 

1965 ), 
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(3.28) 
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dv 


- , + v . yf + — . v f 

$t ^ M dt v M 




where C represents the contribution of elastic collisions. These are 
the continuity and heat flow equations (Braginskii, 1965 ), 
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dt 
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2 dt («P) + |TrPV 


V +P: W + V 

rj) 



( v ~ 2 Z'Xi)^ c 


(3.29) 


The dyad notation will be used for dot products of adjacent vectors and 
tensors. The continuity equation is equivalent to the mass conservation 
law, the first expression of ( 3 . 4 ), used by Newcomb (1962), but the heat 
flow equation is not in the form of a conservation law such as the second 
expression of ( 3 . 4 ). 

Virtual displacement method : A modified version of the virtual dis- 

placement method used by Serrin (1959) and Lundgren (1963) will now be 
employed to derive the variations of L from (3.29). This method consists 
of introducing an arbitrary parameter, e, into n (x,t;e), etc. The 
corresponding nonlocal and local variations are then defined as, respectively, 


6'( ) “ (jo** fixed, e = 0) , 

&( ) = fixed, e = 0) . (3*30) 

For our problem, it is more convenient to use t in this role, without 
introducing e. 
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The nonlocal and local variations are defined in figure 3,4, from ‘ 
which we can write 


6 / n(x^t) = n'(x',t) - n( x, t) . 

^ ' 

6 n(x,t) = n'(x,t) - n( x, t) , ( 3 - 3 l) 


with n' denoting the density after the virtual displacement of 



(3.32) 


Variation in v : A schematic definition of 6 'v„ and 5v is 

shown in figure 3 . 5 > from which we obtain immediately 


. > / dx* dx d? 

6 v = v'(x , t) - v__(x, t) = — = 

~-*D ^ --D ^ * dt dt d t ^ 


dj 


'So-St -i-*D - 


(3.33) 


where (3.32) has been used to obtain the first expression, and Taylor 
series expansion of (x. , t) at x has been used for the second 

expression. 


Variation in n : We shall now use the first expression of (3.29) 

to obtain the virtual displacement in n. First, we perform a real 
displacement, Ax ( 39 1 ) , of a plasma cell in a short time period, At, 
with the particles of the species in the cell conserved. This operation 
is illustrated in figure 3.6, which shows both Ax and F. Applying 
the first expression of (3.29) along the path, Ax, we obtain. 




Figure 3.** The definition of nonlocal and local varia- 
tions in n due to the virtual displacement 
£(x, t)* The nonlocal variation, 6 *n, is 
defined by comparing n * with n for the 
same cell, while the local variation, 6 n, is 
for the same coordinate before and after 

the virtual displacement. The trajectories 
of the same cell before and after virtual 

displacement are denoted by and , 

respectively. 
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Figure 3-5 The definition of nonlocal and local variations 
in due to the virtual displacement 5(x.t). 

The nonlocal variation, 6 v , is obtained by 

'mJ 

comparing v ' with v_ for the same cell, while 

'vD ^1} 

the local variation, 6v , is for the same coordi-- 
nate x. The trajectories of the same cell 
before and after virtual displacement are denoted 

by and , respectively. The virtually 

displaced trajectory that passes through (x, t) 
is denoted by — • — * — . 


5i 



/ 



Figure 3.6 The approach of a real displacement, Ax, performed 
in a short time period. At, to the virtual dis- 
placement, as At diminishes to zero. Since 

(3.29) describes the changes in n and Tr P/2 
along the path Ax, the virtual displacements in 
n and Tr P/2 are then described by (3*29) as a 
special case At -* 0 and A x -* 
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- 0 


(3.3*0 


An Ax 

At At 


An = n / (x + Ax, t) - n(x, t) 


Since Ax and At are arbitrary, we can let Ax -♦ § when At 0. 

Hence by definition. An, the change of n in the real displacement, 
approaches the nonlocal variation 6 ; n. By cancelling At in (3.3*0 
before approaching the limit At = 0, we obtain 

An -» 6 = -n V • j* 

The local variation in n is obtained by subtracting £ . Vn 
giving 

6n = -V . (§n) 

Variation in (Tr P/2 + V . q) ; The foregoing procedure can be applied 
to any equation linear in d/d t and such as the second expression 

in (3.29). Caution must be exercised, however, for those terms which do 
not contain Vp explicitly. From (3.29), we have along the displacement 
Ax , 


( 3 . 35 ) 

from 6^1, 

( 3 . 36 ) 


\ A(Tr P) + | Tr P V • (Ax) + P g :V(Ax) + (V . q)At 


= (f | v C d zj At - m | 


V - ( Ax)c dv 


( 3 . 37 ) 


When At 0, the term (mJv^C dv/2) At will vanish if the macroscopic 
displacement hx does not significantly alter the microscopic random 
processes within the cell* The term (V. q) At does not vanish: the 
identity 
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(3.38) 


V • q 



d_ 

dt 


( V . Q) + Vv • VQ 


} 


implies the relation. 


(V.jq)At = A(V • Q) + V(Ax):VQ . (3.3 

Substitution of (3.39) into (3.37) before taking the limit At = 0, and 
& — & t gives the following nonlocal variations; 



Tr P + 

''Vm' 



I -* 1 v&'S-I-J. 


m v C dv 


(3.40) 


The local variation corresponding to (3.40) is obtained by subtracting 
% ' V ( Tr P/2 + V . Q) , 


Tr P + V . Q 




m v C dv 


(3.4i) 


The Euler- Lagrange equation with respect to x: We are now in a 

— ■ - • — 

position to derive the Euler-Lagrange equation from L E by variation with 
respect to x. The starting point is the integral of L E from time 

(3.42) 

Because the integration variables are x and the time t (for fixed x ) 
local variations are to be used. The vector ?(x, t) is assumed to vanish 
on the boundary of V, and at times ^ and tg. Using (3.2?), the 
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variation of (3. 42) becomes, 




dt 6L 


E 


6L E = Jd* jfinjjjf v D 2 - q(cp . A)j + n 5 ^ • (mv D + qA) - 6^ Tr P+ V . q|J . (3.43) 

Substituting (3.33), (3.36), and (3.4l) for 6v , 5n, and 6(Tr P/2 + V. Q) , 
in (3.43), integrating by parts those terms involving (3^/5t) and V§, 
dropping any surface integral that involves § on the boundary of V, 
and dropping any initial (t^ and final (t g ) terms that involve £, we 
obtain from (3.43) by straightforward manipulation the following result. 


' 6-l e = 


J *s 1 • + io ■ %) - v - 


P- qtilvtp + - 


) 


+ qn( VA • - Vjj • VA) - + qA)[|| + V . (nvj + Jm v 


c dv 


(3.44) 


Use of the dyadic relation. 


v X ( V X A) = VA . v - v_ - VA , (3 4 K ) 

and substitution of (3.7) and the first expression of ( 3 . 29 ) in (3.44), 
finally gives the Euler- Lagrange equation as 

d V D /• 

' m dr +7 ^ =qn ^^ D x |) + jmvCdv , (3.46) 

which we recognize as the macroscopic force law, with general pressure 
tensor and elastic particle collisions (Braginskii, 1965) . 
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In the above derivation. Q is not restricted any further than by 

its definition in (3.l6). The form V*Q in the integrand of L is 

E 

not sufficient to make Q disappear from the final form of 5L , since 

Q is in general nonzero on the boundary of V. In arriving at (3.^-)> 

E 

however, it can be shown that Q appears in the integrand of 6L only 
in the form, V ♦ . VQ ) , and hence drops out since £ vanishes on 

the boundary of V. 

3*3 Hamiltonian Density 

When the independent (integration) variables of a Lagrangian, L, 
are distinctly different from the generalized variables, and L is qua- 
dratic in the time derivatives of the generalized variables, the 
Hamiltonian, H, can be obtained by straightforward application of the 
well-known Legendre transformation (Goldstein, 1950) . The canonical 
Hamilton system of equations is then obtained by varying each of the 
generalized variables, and their conjugate momenta, separately. 

For L of (3-22), however, the dual role of the variables calls for 
additional care in formulating the correct Legendre transformation. If 
the momentum conjugate to x is defined as 

3e = f^ n sn ( m ^) + q £) > (3-^7) 

by the use of (3.25) in the Eulerian coordinates. 7] „ will be dependent 
on because jc is the integration variable. This violates the re- 

quirement of the Legendre transformation that 6jT] E and 6x (= are to be 
independent of each other. A way to avoid this difficulty is to trans- 
form the integration variable from x to x^. for that part of L 
describing the particles, so that the integration is effectively performed 
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at ts 0, The dual roles played by A and 9 require that their 
dependence on be only through their explicit dependence on x (x^, t) . 

The momenta conjugate to the electromagnetic potentials must be defined 
only in the Eulerian coordinates (x, t) (see Section 3.3.2). 

The validity of the foregoing modifications to the Legendre trans- 
formation will be established in Section 3.3.3j where the Hamilton equa- 
tions will be shown to agree with the Euler -Lagrange equations obtained 
in Section 3.2.3. After making similar changes in the definition 
of the Poisson brackets, the energy equation is obtained correctly, 
proving that the canonical formulation is self-consistent. 

For completeness, a Hamiltonian corresponding to the Low Lagrangian 
will be obtained in Section 3. 3. - 6 ,... The resulting total microscopic and 
macroscopic energy conservation equations will be compared to help clarify 
the meaning of the macroscopic potential, 0. 


3.3.1 Lagrangian Coordinates 

In deriving their variational principles, Katz ( 1961 ) and Newcomb 
( 1962 ) used the transformations between x and x n extensively. These 
transformation relations include, for x (x_, t), the Jacobian, J, the 


transformation matrix* x_, * and its inverse matrix, 

i j 1 


K, 




j - 


&x 




bx. 

1 

x. . = 5 — 
1J ox 


Oj 


bx . 

K. . = 

bx ± 


(3.48) 


whence it is easily shown that 


bJ 


5x. 


= JK. 




3 1 


K . x. . = 6 , 

ki ij 3 k 


a(JK,J 


Ji' 


3x 


Oj 


= 0 , dx = Jd^ , (3.49) 


\ 
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where 6 is the Kronecker delta with respect to the Cartesian 
coordinate indices i and j, and the summation convention has been 
implied. The initial value of n(^t) is n^^x^), By conservation 
of particles along the cell trajectory, we have 


ndx =s 

n O d 35o * 

n o = 

Jn 

Under these transformations, 

the Lagrangian 

> L, 

of (3 

lT f 

dx. X + fdx/ 

V* 

_ B 2 \ 

— ' J 

v o 

—Us J — \ 
v v 

2 

2 [ 1 o) 


becomes 


(3.50) 


X S = I n o * 2 - qn o (cp -£’^ - J (| ^ q) • (3.51) 

The time derivative, x (= dx/dt), used above is taken with fixed x . 

~0 

The x, _de P end e n ce of A and cp in the integrand should be retained, 
in keeping with the dual roles of these variables. 

The Euler-Lagrange equations with respect to A and 9 need not 
be explicitly rederived in the Lagrangian coordinates { xq ? t) , They 
should agree with those in the Eulerian coordinates since Hamilton’s 
principle is invariant under transformation of integration variables 
(Courant and Hilbert, 1953). The variation of L with respect to x, 
however, should be re-examined, since x is no longer the integration 
variable. To do so, only the part L , corresponding to a particle 
species, 

L s“J d £) £ s ’ (3.52) 
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has to be considered. The standard form of the Euler-Lagrange equation 
for x is then 


d_ 

dt 




(3.53) 


In Lagrangian coordinates, the functional derivatives of L with 

s 

respect to x and x have the forms 








(3.5M 


which may be compared to those for SlL/jA^ etc. in Eulerian coordinates 

/“w 

in (3.25). 

The only part in X whose partial derivatives with respect to x. 

s i 

and x are unknown at this point is Tr P/2 + V .Q. To obtain them, 

1 J 7 

we use the heat flow equation, but in a different fashion from that of 
Section 3.2.3. In Lagrangian coordinates, (3.29) becomes 


-i-^-(Tr P) + Tr p '?.x+P:7x+7 



2v • x)c dv 


(3.55) 


where the gradient, V, is to be read as K 5/dx^ . Replacing V. q 
by the dyad identity (3.39) converts (3-55) to the form 



(3.56) 
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where the third formal expression is obtained by treating Tr P/2+ V . Q 

as a uniquely determined function of x, .(t), x.(t), and t. Thus, 

i J ^ 

by equating the coefficients of x. . and x., we obtain 

ij i 7 

3 (l \ i SQ 

■w- — k Tr P+ V.Q = - - Tr PK..-P K.. -t-K, , 

lj ~ ~/ 2 ~ ji xk jk 3x. jk 7 


|^- (| Tr P + V . = -Jm v. 


C dv 


a_ 

3t 


(l s)=Jf 


c dv . (3.57) 


These relations enable us to find expressions for iil^/Sfx. by deriving 
d£ s /§x^ and d£ g /3x_ from (3*5l). After some algebra, we obtain 



C dv 


- J 


K 


jk 3x 


Oj 


P + 
ik 


K.. 

Ji 



(7 . Q) - K 


jk 3x, 


Oj 





Zl , • . 

g* = v m ^ + qAj 


(3.58) 


By replacing 3/3 xq^. by 5/3x^, we find that the two terms involving 

Q in the first expression in (3-58) cancel each other. It is now 
straightforward to use (3.4$)), ( 3 - 50 ), and (3.58) in (3-53) to show that 
the resulting Euler -La grange equation is identical to (3.46), except 
for a multiplicative factor, J. 
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3.3.2 Modified Legendre Transformation 


Since jx is no longer formally used as an integration variable in 
( 3 . 52 ), the momentum variable conjugate to x for each particle species 
can be properly defined as (Schiff, 1968 ) 

8fL 

^ n 0 Cmx+ qA) , (3.59) 


where the second expression of ( 3 . 58 ) has been used. 

The momentum variable, a, conjugate to A, must be defined 
differently. It is reasonable to let a have exactly the same dual 
role property as A has in L. The way to assure this is to make 
(*>*)* rather than (*^*t), the independent variable in the definition 
of a. Thus, we have — 



( 3 . 60 ) 


which, by using (3.22), may be put in the familiar form (Schiff, 1968 ), 



(3.6l) 


A momentum variable conjugate to 9 does not exist, because L does 
not involve Scp/9t. 

Since the conjugate momenta of (3-59) and ( 3 . 61 ) are defined in 
Lagrangian and Eulerian coordinates, respectively, the Hamiltonian must 
be defined by adding contributions separately from the plasma variables 
(*>3,) } and th e electromagnetic field variables (A, 0,9). it is given 
by the following modified form of Legendre transformation 


6l 



( 3 . 62 ) 


H = Ef 


SA 

- a- l » 


0 


which is assembled in a similar form to the Low Lagrangian, L in 
(3.6), and the macroscopic Lagrangian, L in (3*5l)« Substitution of 
L from (3.5l)> and elimination of k and SA/dt by use 0 f ( 3 . 59 ) 
and (3.6l), reduce H to the form. 


"= Ej + \ *Mr 

v L U 


V X A j ‘ 


2 k 


- a . vcp 


0 


0 


M = 0 

s 2 mn 


1 31- i n oil 2+ q V + J (l 


Tr P + V 


(3.63) 


Alternatively, H can be transformed to Eulerian coordinates by use of 
the last expression of ( 3 . 1 * 9 ) and the first expression of ( 3 . 50 ), 


1 

T| 

2 mn 

- - qnA 

J ^ 


= J dx H , 

V 

•] 


2 1 I a 2 |VXA|' 

+ o Tr P+ V • Q+ qnqp + < ±~+ 

~ * 2e_ 2u. 

0 


a • vcp . (3.64) 


By use of (3*59) (3.6l), we have from (3*6^) 

p ) 


2 

H \^(mn 2 1 \ e 0 E B 2 

* = 2^\T V D + 2 Tr f +v *Q + WV) + + §£ + e oS • VC P , (3.65) 


which is the familiar form of energy density for a plasma, plus the term 
7 • Q- The latter represents all of the energy density lost through 
heat conduction prior to time t, by a cell of a particle species, as 
it moves along its trajectory. 
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It is of interest to compare H with that given by Schiff for a 
quantum mechanical system of charged particles in an electromagnetic 
field [Schiff, 1968 , equation (57-^)] 


H_. = ldx H 


QM 




QM 


2 | 7x A| 2 

*QM = b - *»a)y| 2 + V Q j¥j 2 + q<p|¥| 2 + ~+ Q ,,~ - O • *P , (3.66) 


2e, 


2u 


0 


now written in MKS units. In (3.66), ¥ and are the Schrodinger 

wave function and potential, respectively, and -iiiV is the particle 
momentum operator. With the understanding that |y| 2 is equivalent to 
particle number density, a term-by-term correspondence is observed 
between (3.66) and (3.64), with the potential energy density [ ¥ | 2 

corresponding to Tr V/2 + V • (3, in agreement with the interpretation of 
the negative sign that the latter has in the macroscopic Lagrangian. 


3.3.3 Hamilton Equations 

The approach of Section 3*3.2 has been plausible. To establish 
that the formulation is self-consistent, we should verify that the 
Hamilton equations derived from H agree with the Euler-Lagrange equa- 
tions derived from L. 

Using (3.59) and (3.60), the Euler-Lagrange equations, (3.53) and 
(3.24), become 


HL . j0L 

gx & > M. - at ' w - 


(3.67) 


The variation of L of (3.5l)j 
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E J <&>(§ • s* 1 ' ■ i) ♦ ^ I H ’ < 8 - 68 > 


can be put in the form 


6l ‘EJ + 


%] . 

( 3 . 69 ) 


where (3.59)? ( 3 . 60 ), and ( 3 . 67 ) have been used, and 6 ($0 is the 
local (nonlocal) variation defined in Section 3.2.3. The variation of 
(3.6?) consequently becomes 

5h = -EJ %(ii-s-io' 5 'a)- f • (3.?o) 

v o 

Since the variables in H to be separately varied are x, Tl, A, a, and 
9? it is convenient to write &H as 



where the functional derivatives of H are analogous to those of L 
given in ( 3 . 5*0 and ( 3 . 25 ). 

Comparison of (3.70) and (3.7l) shows that the correct forms of the 
Hamilton equations, for macroscopic plasmas described in terms of variables 
with dual role, are 


• 3h _ jzSH 

~ “ 3T] » at ■ m > 

which are identical to ( 3 - 59 ) and ( 3 . 61 ), and 


(3.72) 


6k 





(3.73) 



3h 


3x 

. ^ 





The last two expressions of (3.73) lead to (3.26) by straightforward 

differentiations of (3.64) before using (3.6l). The right-hand side of 

the first expression of (3.73) can be written as -dtt /5x. + 

(d/Sx )(BJi /3x ) . Its first term can be derived from (3*63) as 

u J s 1 J 



- 2- (Tl.-qn-A.) 
m J 0 J 


5x. 

i 


+ "-o §7 ‘ J j mv i 


C dv 


(3.7*0 


where the second expression of (3*57) has been used* Its second term 
can be shown to be 




BP 


~JK 


Ik 


jk Bx 


0j 




> 


(3.75) 


where the first and third expressions of ( 3 . ^9 ) j and the first expression 
of (3.57^ have been used to obtain the first result, and the relation, 
a/a Xi ~ b/bx^y has been used to obtain the second result. It is 

then straightforward to show that the first expression of (3.73) is equiva- 
lent to (3.46) by the use of (3.75), (3.74), (3.59), (3.45), and the 

relation, x = v . 

~ '-"-D 

3.3.** Poisson Brackets 

Poisson brackets are convenient in writing down the time derivative 
of a physical quantity in a system that is formulated in terms of a 
Hamiltonian* They must also be modified on account of the dual role of 
the variables* Consider in general a physical quantity G, which is a 
functional of x* T), A, <x <P, and their first derivatives* Without 

i-s S r~wi r-^/ 
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loss of genera lity, we may write 


G= Jdx§=£ j **Q% + J d ^§EM ' ^ 3 ' 76) 

V V Q V 

where is a function of A, a, cp, and their first derivatives, only, 

and contains terms that must involve x, Tj^ their first derivatives, 

and Tr P/2 + V *Q. The time derivative of G can be obtained by 
accounting for the dual roles played by the variables, x, T\j O, 
and cp, in G . This requires that 

1-El <*§ . (3.77) 

v o v 

where, in the right-hand side, the first time derivative is defined with 
A, a, and cp treated as given functions of x, and the second time 
derivative is defined with x and T| treated as constants. Use of the 
Hamilton equations, (3-72) and (3.73), transforms (3.77) to 

E$ %% I'fe-la. > (3-78) 

s o s 

where the plasma and electromagnetic field Poisson brackets, and the 

functional fluxes, C J and T, are defined by 

—EM 


£=t?*EM 


x,T) 


§>» 


'a, a, cp 
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U >, I f . llG 3 h 2 g 3 h \ 

f 8s ' Ms !v.n J * W ’ 


x>lL 


'0 


S, 


j - f h J ?g .Is * 9 . . 2”+ *p\ 

\a,<p" ) % *'* % *>*) ’ 

rZ7 r+*S V 


58 s ■ 8§ S 

*sj 3x. . X i + 3(3T] /3x ) ' 

1J v 'i Oj 


-? 


r a§ ag as a? 

®ij a(aA./a x .) at a( aoj.yax ) at a( a^/ax ) at 

* j ^ j j 


aa. 


(3.79) 


In the integrals of (3.78), and S are the bounding surfaces of 

Vq and V , and do^ and da are the elementary surface areas, with 
vectors taken along the outward normals. The partial time derivative, 
SG/at, takes care of any explicit time-dependence that G may have in 
addition to its implicit time-dependence through x, T|, A, a, and cp. 

3.3.5 Energy Conservation Equation 

Equation (3.78) can be used to formally separate the time deriva- 
tive of any physical quantity into the bulk and the surface contributions, 
and the plasma and electromagnetic field terms. As a demonstration that 
the foregoing formulations are self-consistent, the total energy con- 
servation equation is derived by substituting H for G in (3.78). 

We obtain 



5_ + 



( 3 . 80 ) 


where the plasma and electromagnetic energy flux densities, J and 

~S 

3^, respectively,, are given by 
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(3.8i) 


A more familiar form may be obtained from (3-80) by use of the third 
expression of (3.^8) and (3.^9) > and the second expression of (3.50). 
After some manipulation, we have 


Ux - 2(7 • Q dx)l + dx(V • 3) = 0 , 

9a 

J = 2(P • v + q) - - ^ X B - e E § . (3.82) 

^ ^ ^ f^o ot ~ 0^ 3t N 

By use of the identity. 


It (“*> - ■fe( 7 'I D ) - (3.83) 

and the differentiation scheme of (3.77)j we can further transform (3.82) 
to the following form in Eulerian coordinates, 


|- ( Jt - 2 v • Q) + 7 • S 0 , 

ot ^ ~ 



1 

2 


Tr 


P 


+ p 


^o + . 




EX B - 


0 


at 


>05^) 


(3.84) 


Since U - Z is the energy density in the plasma at time t, (3.8l), 

to which (3,84) is equivalent, is the correct expression of energy con- 
servation in the mixture of Eulerian and Lagrangian coordinate systems. 
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It is 


( 

also possible to derive (3.84) by evaluating 3Ji/3t from ( 3 . 65 ), 
and then using (3.20), (3.46), and Maxwell's equations. 

The derivation of (3.84) completes our development of a self- 
consistent Hamiltonian description of a plasma in the macroscopic 
approximation; the difficulties associated with the dual role of the 
variables in the description of the field-particle system, originally 
pointed out by Sturrock ( 1958 a), have been resolved. Since the Low 
Lagrangian also assumes 'dual role, the microscopic Hamiltonian descrip- 
tion of a Vlasov plasma can also be established in a similar fashion, as 
will be shown briefly in the next section. 

3.3.6 Microscopic Hamiltonian Density 

By the use of the Low Lagrangian in (3.6) and (3.10)> the conjugate 
momenta can be defined in forms analogous to (3.59) and ( 3 . 60 ) as 


2Sl t 


bt 


J&SvXc’V - W = H s "£ + ' 






a( ?c, t) ji( QA/^t) a(3A/dt) 3x 


d 


ax. 


3 ( o A/dx^t) 


= e Q vcp + 


3A 

at ' 


(3.85) 


where the second expression for C follows because -t does not involve 

3 / 3 x„ and a/3v_ terms. 

'* vJ J 

The corresponding Hamiltonian is defined analogously to H as in 

(3.63), 





X + 


J 


dx a 


3A 

Ft ' l l 


( 3 . 86 ) 
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After some manipulation* this takes the form 


\ = f( 2 o'V h + J 


Vx A 


2e o 2|a o 


- a • V<p 


h = b li- q Af 2 + qc P • 


(3.87) 


The Eulerian form of H can be obtained from (3.87) by use of (3*9); 


\ = j d £ “l < 


» L “Dj d ~ t fei> t )(l v2 + <&) 


s E 2 

+ “i“ + |^ + € oS* vcp • (3.88) 


With the definitions of (3.85); the Euler-Lagrange equations of L T 


are 


2SL. . %l Ba 

L _ ^ L 

~ * * M ~ at 


j&x 


* l l 

29~ 


= 0 


( 3 . 89 ) 


These can be shown to represent the microscopic particle force law and 
Maxwell's equations. Taking the variation of H , and using ( 3 . 89 ), 

Lj 

the Hamilton equations are easily obtained as 

2 h 3 a 3 h . SfH da dn, dH r 

i k k r _ k ~ L L _ 

~ " JfC ' 3t aa ' * " " 3x » at = " 2A 1 j8cp “ 0 

(3.90) 

Since H does not involve the derivatives of x and C with 

■Lr 

respect to and it suffices to consider a physical quantity. 
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\ 


G, for plasma with the same simplification. We may write 


C= j*dx§ (x,^; a,A,cp) 

V 

= d 3) j% ^Zo’Zq* § s ( Z>& + J d £ § m ( ^ 9) • (3.91) 


0 


The expression for dG/dt analogous to that of (3.78) is then obtained 
as 


f - f°'\l p + 


J 


da • F 


-EM 


(3.92) 


where the microscopic Poisson brackets, and j > are defined by 




g, h ! = S|s\ , 

l *em J ~\*£ 5 Q *P ht f 


_* _ j~ . ^ ^ . as , 39 as 

emj at aovax ) at a(&*/ax j at a(a<p/ax.) 

J J J 


(3.93) 


The energy conservation equation in the microscopic plasma model is ob- 
tained by replacing G by H in (3.92). We obtain 

Ju 


f!i+ f 

dt J 


da • ? L = ° , 


J. = - i- £ X B - V E If . 
[~Iq 9t rs* 0-^ Si 


( 3 . 94 ) 


By use of (3.88), the following is obtained from (3.9*0 
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Figure 3 . 



( 

\ ^ / 



d cr 

P 4 dcr, 

IV v 

q*dcr 


(b) 


7 Hie conceptual difference between the cells of the (a) 
microscopic and (b) macroscopic models in a collisionless 
plasma. The arrows in (a) indicate the motions of the 
particles that carry the boundary surfaces of dx with 
them. In (b)^ dcr is the surface area of dx with an 
outward normal. The cell; dx, in the macroscopic 
model, can only move with the macroscopic drift, v 
The particle thermal motion then gives rise to momentum 
and energy transfer across dcr. 
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4r (dx n ) + dx(V • 3 ) = 0 
at ~ JL 


(3.95) 


A comparison of the energy conservation equations in the microscopic 
and macroscopic models, (3-95) and (3.82), respectively, can be found 
in Section 3.^.3. 

3.3*7 Entropy for Plasmas 

In this section the Lagrangian and Hamiltonian formulations have 
been developed, so far, for the purpose of obtaining self-consistency and 
consistency with the well-known equations describing plasmas in the 
macroscopic model, e.g. the force law, Maxwell 1 s equations, and the 
energy conservation equation. Now, the dual role approach has resulted 
in a canonical formulation of plasma dynamics in which the plasma and 
electromagnetic field variables are all treated as canonical variables. 

It is proposed that these canonical variables may also be effectively 
used to express other plasma quantities of interest, such as the appro- 
priate plasma entropy. 

The question of whether the Gibbs or the Boltzmann H-f unction is 
the correct definition for a many-particle system with arbitrary inter- 
particle forces, has been discussed by Jaynes ( 1965 )* He concluded 
that an appropriate definition of entropy must include all of the degrees 
of freedom that we intend to use for the physical system. In this 
section, the number of degrees of freedom equals the number of generalized 
field-like variables. Thus an appropriate expression for entropy in the 
macroscopic plasma model can be obtained only in terms of all of the 
generalized variables and momenta. Since the Boltzmann H-function, in 
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its hydrodynamic form (Braginskii, 1965 )? 


S =2jdx S n , S = | 'En (Tr P) - | In. n + S Q , (3.9^) 

is defined without including the degrees of freedom associated with the 
electromagnetic field variables, it clearly does not represent the total 
entropy of a plasma in the macroscopic model. 

Consequently, the Boltzmann H- theorem is likely to be violated in 
plasmas whenever the process involves an appreciable change in the electro- 
magnetic field energy. This conclusion agrees with recent results of 
Jaynes (1971), who has shown that the necessary condition for violating 
the Boltzmann H-theorem is that the initial kinetic energy associated 
with particle velocities, both drift and random, be greater than the 
equilibrium kinetic energy. The difference is then converted into 
potential energy, associated with particle interaction forces, both 
collective and random, during the evolution towards equilibrium. 

An important question which we have not had time to pursue is the 
nature of an appropriate expression for plasma entropy based on the 
canonical models established in this section. 

3.4 Discussion of the Macroscopic Potential, V 

In this section, we have been mainly concerned with a rigorous 
verification of the Lagrangian and Hamiltonian formulations for a plasma 
in the macroscopic model. However, in order to obtain the correct 
macroscopic Lagrangian, a heuristic argument that assumes sufficiently 
frequent particle collisions was used in Section 3.2.2. 
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This assumption was in turn found to be overly restrictive because 
the macroscopic Lagrangian and Hamiltonian, verified in Sections 3*2.3 
and 3.3*3* respectively, are equally applicable to a collisionless 
macroscopic model when the moments of the Vlasov equation are appro- 
priate. Therefore, it is of great interest to find a reasonable 
explanation for the necessity of introducing V of (3.2l) also for 
the collisionless macroscopic plasma model. We shall proceed by con- 
sidering the relation of V to viscosity and heat conduction first. 

3.^.1 Relation of V to Viscosity and Heat Conduction 

According to (3.2l), U represents twice the thermal energy density 
of a particle species plus the energy density lost due to heat conduc- 
tion by a macroscopic plasma cell, as it moves along its trajectory. 

The relation of V to viscosity and heat conduction can be further 
demonstrated by considering the particles of a species in a macroscopic 
cell, dx, that has a size much larger than the mean particle spacing 
and follows some trajectory in the plasma. By the use of (3,83), (3.29), 
and (3*38); the total time derivative of U within dx is given by 


dx -2P:Vv + 


lm(v - 2v 


v )C dv + 


3Q 


The first term in the right-hand side brackets represents twice the time 
rate of heat generation in dx. due to viscosity and mechanical compression 
(Braginskii, 1965 ). The second term is twice the time rate of heat 
generation in a particle species resulting from collisions with particles 
of other species. The third term contains twice the rate of energy 
gain due to heat conduction. 
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If (3-97) is integrated over V, before summing over all species of 
particles, the collisional term will drop out because of conservation of 
energy and momentum in elastic collisions. Equation (3*97) then becomes 


dt 


f 1*-* \ 

C _-,r / n 

J d i(E v) 

- J d ,:Z [- 2 £ ; % + 7 • ( V i * St - 2 &) 


( 3 . 98 ) 


The phenomena of viscosity and heat conduction, which are usually regarded 
as dissipative, are being compensated for. The total energy of the plasma 
in the macroscopic model, as denoted by H in ( 3 .63) “( 3 . 65 ) , is con_ 
served through the inclusion of -V in £' of (3.12). 

Li 

3.^.2 Loss of Particle Discreteness in Applying Macroscopic Approximation 
Equations (3.97) and (3.98) are applicable also in the case when 
collisions are negligible, e.g„ in the macroscopic approximation 
of a Vlasov plasma with the size of dx much larger than the mean 
particle spacing. In this case, the argument of Section 3.2.2 that leads 
to the expression for U is no longer appropriate because the assumption 
of frequent particle collisions within dx is no longer valid. 

The following question consequently arises: why is it still necessary 
to include 1) in X/, in the macroscopic approximation of a Vlasov 

li 

plasma, to obtain the correct macroscopic Lagrangian, L, of (3.22), 
now that U R in ( 3 . 15 ) can no longer be considered as the energy 
associated with particle collisional interactions? This question suggests 
the following argument; that it is actually the loss of particle dis- 
creteness in applying the macroscopic approximation of (3.12) which must 
be compensated by the correction expressed by ( 3 . 15 ). 




In the collisionless microscopic model, the trajectory of each 
particle is, in theory, to be solved. When the Low Lagrangian is 
written, the cell dx is used for the purpose of summing the contribu- 
tions of particles within dx. Owing to the assumed particle continuity 
in phase space, the bounding surface of dx in this model is considered 
completely flexible; if dx is cubic at time t, we can, in theory, 
always deform the shape of dx to include the same particles at some 
later time* This definition of dx is illustrated in figure 3.7(a), 

I-V 

which implies that, in the microscopic model, as long as the particle 
trajectories are solved, there is no momentum and energy transfer 
across the flexible boundaries of dx. In the macroscopic model, however, 
particles of the same species are indistinguishable. Only the spatial 
coordinates of the macroscopic cell, dx, can be used to identify the 
plasma, as shown in figure 3.7(h). Although moving with the drift velo- 
city, v_, the macroscopically smooth boundaries of dx are penetrated 
by particles due to thermal motion* This then gives rise to the well- 
defined outward macroscopic momentum and energy transfers V ■ P dx and 
V • q dx, rexpectively, from dx. Furthermore, since elastic collisions 
conserve momentum and energy, they do not affect these macroscopic trans- 
port phenomena. In terms of these momentum and energy transfers, the 
dynamics of the cell, dx, and that of a cell, AV, introduced in 
Section 3.2.2, become identical* The macroscopic potential energy of 
( 3 . 15 ) must correspondingly be introduced for a collisionless macroscopic 
plasma. 
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3*^.3 Definitions of Plasma Cell Boundary 


That the microscopic and macroscopic models define the boundary of 
a spatial volume differently can also be seen by comparing the micro- 
scopic and macroscopic energy conservation equations, ( 3 * 95 ) and ( 3 - 82 ), 
respectively . We see that the difference between the microscopic and 
macroscopic energy flux densities, in ( 3 - 9*0 anc * 3 F in ( 3 . 82 ), 

respectively, is the flux density, £ (P-v_+q), associated with 
particle thermal motion. Mathematically, this flux density is absent 
i n J T ) because the microscopic Hamiltonian, H T in (3.87). does not 
contain any derivatives of x and C with respect to x„ and v^. 
Physically, this difference between (3.95) and (3.82) is attributable 
to a difference in defining the microscopic and macroscopic plasma cell 
boundaries . 

Consider the time derivatives of H in (3-88) using two different 

Li 

definitions of the volume, V. First, define V as a stationary volume. 
We have, 



Z | d *J d z[S v 


2 bf 

at 


+ q 





(3.99) 


which can be reduced to the macroscopic energy conservation equation, 
( 3 . 82 ), by use of the Vlasov equation. Maxwell's equations, and integra- 
tion over velocity space (Van Kampen and Felderhof, 1967 ). Second, 
define the boundary of V as moving with the enclosed particles, so that 
the Liouville theorem can be applied. We then have 




d y. f(x,v,t) 


d 2of 


d *o 


) = constant , ( 3 . 100 ) 
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because of ( 3 * 8 ). The total time derivative now becomes. 


dt 


£ 

:= s f -H* m ^*d? +q ft tp( ^ t )] + ft i d 4“V + fc +e ^ ,vcp ) ' (s - 

t r tr \ * 


101) 


which can easily be reduced to the microscopic energy conservation equa- 
tion in ( 3 * 95 ) by use of the particle force law, Maxwell’s equations, 
and integration over velocity space* 

Since these energy conservation equations, ( 3 * 95 ) and ( 3 . 82 ), are 
also derivable from their corresponding Hamiltonians, ( 3 . 87 ) and ( 3 . 63 ); 
which are in turn derived from the corresponding Lagrangians, (3.6) and 

( 3 . 5 1 )* it is seen that the difference between the Low Lagrangian, L_ , 

L 

and the macroscopic Lagrangian, L, is to reflect the difference in 
defining the microscopic and macroscopic volumes (or cells) and their 
boundaries , 


3.^.^ Relations Among the Variational Principles of Various Models 

The arguments of Sections 3-2.2 and 3.^.1 - 3.^.3 are helpful for 

the purpose of understanding the macroscopic Lagrangian, L, and potential 
energy density, U. The relation between these viewpoints can be under- 
stood from figure 3.8, which compares the plasma variational principles 
of the collisional and collisionless microscopic models defined in 
Section 3.2.1, and the macroscopic model. Procedure 1A involves dropping 
the random terms defined in (3.1l) to relate to L^, while 

Procedure IB involves only smoothing within the macroscopic cell, dx, 
as described in (3.12), to obtain from L^, by keeping the quadratic 

random terms. Procedure 2 reflects the difference in definitions of 
the plasma volume, V, or cell, dx, in the microscopic and macroscopic 
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Models: 


Euler-Lagrange 
Equations 
Subsidiary \ 
^constraints/ 


Approximate 
Liouville’ s 
Laws 
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Colli sional 
Microscopic 


Collisionless 

Microscopic 


Macroscopic 



Figure 3.8 Relations among the variational principles of various plasma models* Here 
L t and H t represent the total Lagrangian and Hamiltonian, respectively, 
in the collisional microscopic model; and represent the (Low) 
Lagrangian and Hamiltonian, respectively, in the collisionless microscopic 
model; H.P. represents Hamilton’s principle; L.T 0 represents the Legendre 
transformation, and E.C.E. represents the total energy conservation equation. 



Boltzmann 

Equation 













plasma models, as discussed in Section 3-^-3, with regard to the energy 
conservation equations in the two plasma models. In the case of a 
collisional plasma, Procedure 3 represents the heuristic approach of 
Section 3,2.2 that led to the appropriate expression of V in (3-2l). 

For a collisionless plasma, Procedure k represents the argument of 
Section 3.U.2 that discusses the loss of particle discreteness when we go 
from the microscopic model to the macroscopic model, and the dynamic 
equivalence between the macroscopic cells of the collisionless and 
collisional plasmas, so that the introduction of U is seen to be 
necessary also for the collisionless macroscopic model. 

All of these arguments should be considered less rigorous, for the 
purpose of proving the validity of L of (3.22), than the rigorous 
verification by application of Hamilton’s principle in Section 3-2.3* 

These different viewpoints do, however, converge to the conclusion that, 
in general, the macroscopic potential energy density, U, is a result 
of the combined contributions of plasma fluctuations (which is equiva- 
lent to elastic multiparticle collisions), elastic two-particle collisions, 
the loss of particle discreteness, and a redefinition of the plasma 
cell, dx, in the process of making the macroscopic approximation. 

3.5 Discussion 

The principal contributions of this section verification of an 
appropriate Lagrangian density in Section 3-2.3 and Hamiltonian density 
in Section 3-3-2 for the macroscopic plasma model, including a pressure 
tensor, heat conduction, and elastic particle collisions. In Section 3-3-2, 
it was shown to be necessary to introduce a macroscopic potential energy 
density \) (=* Tr P+ V • <3 ) . In order to obtain the Hamiltonian density. 
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it was necessary to modify the conventional Legendre transformation to 
conform to the dual role of the generalized variables. Low’s Lagrangian 
(Low, 1958) shares the same property; the modified Legendre transforma- 
tion was used successfully to derive a corresponding Hamiltonian in 
Section 3.3.6. The modified form of the Hamilton equations and Poisson 
brackets presented in Sections 3.3.3 and 3.3.4 may have further interesting 
implications in the canonical formulation of plasma dynamics. 

The plasma model used in this section is, of course, a macroscopic 
approximation to the microscopic model described by the Boltzmann 
equation, or the Vlasov equation when collisions are negligible. The 
Euler-Lagrange equations, (3.2 6 ) and (3.46), consequently have a one- 
to-one correspondence with those of the Low Lagrangian. The continuity 
of particles in phase space, (x,v), is used by Low to constrain the 
variation with respect to the particle trajectory, x t) > while 

here the first (particle continuity) and third (heat balance) moment 
equations are used to constrain the variations with respect to the 
plasma macroscopic cell trajectory, x^, t) . That only two moment 
equations are sufficient here does not imply that the mathematical system 
is closed at the second order moment equation. What we have done is to 
appiy Hamilton’s principle to a mathematically open system of equations: 
the Maxwell equations, and the moment equations. The latter are not 
truncated, since jq appears in the heat flow equation, and can only be 
obtained by use of higher order moment equations. In practice, trunca- 
tion will be accomplished by making some arbitrary assumption about q, 
typically q = 0. 

Whenever the scale length of variation of the macroscopic quantities 
becomes comparable with the mean particle spacing, the results presented 
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in this section for the macroscopic model becomes questionable* In addi- 
tion, the assumption of elastic collisions neglects excitation, ioniza- 
tion, and recombination phenomena, and electromagnetic radiation 
associated with charged particle collisions (bremsstrahlung) . Never- 
theless, the variational principle developed here is appropriate to a 
multicomponent plasma, including neutral components. 

In later sections, we shall apply our results to a number of plasma 
problems involving perturbation expansions* While these are straight- 
forward for v and n (Newcomb, 1962), difficulty arises in connection 
\ 

with the term Tr P /2 + V • Q in X. Inspection of the heat flow equa- 
tion in (3.29), reveals that it is not self-consistent: even with an 
explicit form given for C, we have to solve for all of the elements of 
£ as well as the components of qj these depend on higher order moment 
equations* Even when truncation at a specified moment of the Boltzmann 
equation is assumed, the heat flow equation is not generally equivalent 
to some constant of motion for a particle species* For the case in which 
q and C are neglected, approximate expressions for P in terms of 
the perturbation in ' x are generally impossible to derive mathematically 
(see Section 4 . 2 * 2 ). It is not surprising that only problems in which 
collisions between different particle species are neglected, and adiabatic 
processes are assumed to occur with either a scalar or axisymmetric 
pressure, have been successfully studied in terms of perturbation expan- 
sions. Such problems form the subject of Section 4 * 
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4. THREE-WAVE INTERACTIONS: THEORY 


4.1 Introduction 

In this section, we shall show how the macroscopic Lagrangian ob- 
tained in Section 3 can be applied to the description of nonlinear 
wave-wave interactions in homogeneous plasmas. After developing a 
perturbation expansion of the Lagrangian, the method of averaged 
Lagrangians will be employed to derive general coupled mode equations. 

In Section 5* these equations will be specialized to a number of cases 
involving parametric amplification of ion-acoustic waves for which 
experimental data are available. 

The use of Lagrangians in describing nonlinear wave phenomena 
was considered by Sturrock in 1961. He showed that in a conservative 
distributed system the time-averaged Poincare invariants, which he 
later extended to the case of field variables (Sturrock, 1962 ), are 
equivalent to power-balance relations of the type well-known to electri- 
cal engineers as "Manley-Rowe relations" (Penfield, i960). These results 
were extended by Whitham (1965) who averaged the Lagrangian in such a 
way as to remove rapidly varying terms, but conserve the slow variations 
in amplitude, frequency, and wave vector characteristic of a wave train 
in a weakly nonlinear medium. Vedenov and Rudakov (1965) used this 
approach to describe the interaction between ion-acoustic and Langmuir 
waves in a plasma. Using Low's Lagrangian (Low, 1958) and its pertur- 
bation approximation, several cases were examined by Suramlishvili, 
who derived the wave coupling coefficients for interactions between Langmuir 
a ^d ion— acoustic waves (Suramlishvili, 196 ^-), between one transverse 
and two longitudinal waves (Suramlishvili, 1965), between Alfven and 



whistler waves ( Suramlishvili, 1967 )) between Alfven and ion-acoustic 
waves (Suramlishvili, 1970), and among magnetoacoustic waves 
(Suramlishvili, 197l) • A variety of general results have been established 
by Galloway (1972) and Dysthe (197^)^ using the averaged Lagrangian 
method. In particular, the former demonstrated that the energy conservation 
equation may be used to derive the Manley-Rowe relations, and the 
coupled mode equations describing wave-wave interactions. He has also 
shown that the wave energy and energy flux follow directly from the 
energy and energy flux terms quadratic in the perturbations, while the 
nonlinear wave coupling coefficient follows directly fran the term in 
the Lagrangian cubic in the perturbations. Similar general results were 
employed by Boyd and Turner ( 1972b), who used Low's Lagrangian to examine 
the generation of longitudinal plasma waves by two high frequency electro- 
magnetic waves in a warm field-free plasma, and the interaction of 
three electromagnetic waves in a cold magnetized plasma. This approach 
will be followed in Section 4.3. 

Generalizations of the averaged Lagrangian method to nonlinear plasma 
phenomena other than wave-wave interactions have been made. Using a 
hydromagnetic Lagrangian derived by Newcomb (1962), Dewar has shown that 
the interactions between the wave and the slowly varying background 
plasma may be derived by use of Hamilton's principle (Dewar, 1970) . 

Dougherty (1970) has obtained ray tracing and coupled mode equations, and 
demonstrated the conservation of wave action in a relativistically 
covariant formulation. Four-wave interactions, self-action effects, and 
sideband decay phenomena have been treated in the paper by Dysthe (1974) 
mentioned above. Derivations of the wave kinetic equation have been 
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supplied by Suramlishvili ( 1964 ; 1965) and Galloway (1972) ; the latter 
has also provided a description of quasilinear wave-particle inter- 
action. A Lagrangian theory for nonlinear wave-packets propagating in 
a collisionless plasma has been developed by Dewar (1972) which describes 
the nonlinear frequency shift, Landau damping, and modulational sideband 
instabilities . 

It will be clear from the foregoing remarks that the averaged 
Lagrangian method has been utilized to discuss many of the significant 
nonlinear phenomena occurring in plasmas. It has the merits of concise- 
ness and efficiency in the analysis, and gives considerable insight 
into the physical mechanisms involved. 

4.2 Perturbation Approximations to L 

Prior to its use in Section 4 . 3 ^ we shall expand the Lagrangian, 
(3*22), in terms of perturbations in the generalized variables, x, cp 
and A, up to the third order. Similar to the well-known approach for 
the macroscopic Lagrangian (Newcomb, 1962), the perturbations in n, 

Vp, and Tr P /2 + V • CJ will be considered as due only to those in the 
generalized variables. Denoting the perturbed and unperturbed Lagrangians 
by L and L, respectively, we write 

L ' = L + ^ + L 2 + L 3 + . . . , L. = J dx I. , ( 4 . 1 ) 

V 

where denotes the Lagrangian ith order in perturbations. The qua- 

dratic Lagrangian, L^>, will provide the first order equations that 
determine the linear wave properties, and the cubic Lagrangian, L-,> will 
lead directly to the nonlinear coupling coefficients among the waves. 
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where (x, <p. A) and ( §, cp_, A 1 ) are the unperturbed and perturbation 
variables, respectively. Figure 4.1 illustrates the definition of 
which is analogous to the particle displacement vector defined by Sturrock 
(1958a) and Low (1958), and the macroscopic cell displacement vector 
defined by Newcomb ( 1962 ). The quantities (9> A) and (cp^> A^^) are 
considered as functions of the Eulerian coordinates consistent 

with their dual roles discussed in Sections 3.1 and 3.2.T. 

The perturbed Lagrangian differs from the unperturbed Lagrangian 
(3.22) due to two factors: first, the presence of and A^ at 

a given and second, the use of the integration variable x' rather 

than x. Thus the perturbed Lagrangian may be written as 


1 ' = | 6x' l' , 


2 - \ I* P' 

d ~ 


- V -Q'-qn^cp'- * A'j] 


S 0 E " 2 B ' 2 

+ ^ 2 ~-% ' < 4 - 3 > 


where the summation applies to the particle species. In (4.3), and what 
follows, primed and unprimed quantities are functions of x' and x, 
respectively, unless otherwise noted. 

Imposing the condition that the particles are conserved in the 
process of applying perturbations, implies that. 
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Figure 4.1 Hie definition of perturbation § in the macroscopic 
plasma cell position from x to x \ The perturba- 
tion is performed in such a way that the number of 
particles of a species in the cell is fixed. The 
perturbed and unperturbed cell trajectories are 

denoted as and , respectively. The cell 

velocities before and after perturbation are denoted 
by and (£%*)* respectively. The velocity 

of the cell that is at (j^t) after perturbation is 
denoted by v^x, t) . 
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while the Jacobian, J, of the transformation from x to x ' is 

I d *'l 



(4.5) 


The integral in (4.3) can consequently be written in terms of an 
x-integration as 



4.2.2 Expansions of v', n and P * 

The terms that involve v_T, J* A ' E and B ^ are easily expanded 

by the use of (4.2), (4.4), and (4.5), and the Taylor series expansions 
at x. However, expansion of the terms, Tr pV 2+ V • Q in (4.6) 
introduces considerable difficulty* We shall confine ourselves to the 

case of scalar pressure and adiabatic processes satisfying the equation 

/ 

of state 



(4.7) 


for each particle species, with Y being the adiabatic index 



Expansion of v^: Taking the time derivative of the first expression 

of (4.2), we obtain 


~D = + ^ > 


(4.8) 


where % denotes d5[x(x_, t), t]/dt. Then the nonlocal expansion of 


^ is seen to be, 


v_ ~ v_ + v_ + v n + v 0 + 

'*'-'1 J rvj 


V =r f! 


Za -3 ‘ ' 0 


/k rV\ 
*7t 


If we define the local expansion of v'(x,t) as 

1 — ^ 


V '( X, t) = v + v Tn + V + V _ + ... , 

~DW > ~D ~L1 ~L2 --'"LB J 


(4.10) 


it can be shown by Taylor series expansion of Vp(x/, t) at x that 


X 

!i = Zu + 2 jr ( I ' 7)1 ~ l , i-i 


As a result, we have 


JLi-£-£-*b i I L2 = -£' v i + • ^- 12 ) 


Expansion of n^ : Combination of (4.4) and (4*5) gives the particle 
conservation law along § as (Newcomb, 1962 ) 


n J = n 


[4.13) 
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With the Jacobian, J, expressed in its exact expansion, 


% 

Sx 


1 + <7. 1+ l[(v 5 ) 2 -7|:7|] + El , 

‘ 5 ( v -i) 3 -! ( v + 3 h • v l) 


) 


the nonlocal expansion of n * becomes 


n ss n 


n i + n 2 + n 3 + 


n, = -nV • § 
X ~ 




By use of a relation between nj^ and n Li analogous to ( 4 ,ll), the 
expansion of n' then takes the form 


n'(x,t) = n + + ... , 

”li » ”L2 “ T\ 


The results of ( 4 .l 6 ) agree with Sturrock’s generalization of the 
Lagrange expansion in the form (Sturrock, 1960b) 


f-l ) 1 

n Li = i! V k * ‘ ' ' ^k * ‘ ‘ = 1 > 2 , 3 ; i = 1*2, ...) , 

while n L1 and n^ have been derived by Newcomb (1962). 


(4.i4) 


(^.15) 

local 


(^. 16 ) 


(^. 17 ) 
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Expansion of P The nonlocal expansion of p' can he obtained by 
substituting (4.15) and the second expression of (4,l4) in the first 
expression of (4.7), 


P = P + P 1 + P 2 + P + ... , 


Pj.-YPfVI) , p 2 = |p( Y (v._|) 2 + 7|:V|] , 

P, = (V • 5 ) 3 + ^ {V • 5)(V§:V?) + 4 • ??)1 

W r d <"0 ^ j r</ N /C I 


/I. *. Q \ 

J 


By use of a relation between P i and P analogous to (4.1l), we obtain 
from (4.18), 


P'(x,t) , p + P L1+ Pl2+ ... , 


P T . = -YPV • § - ? ♦ VP 


2 

P L2 = P [l~ < 7 *i) 2 + I YJ • V(V • ?)] 

+ [l ‘ 7 I + Y ( 7 ‘ J)?] “ 7P+ | J?: wp . 


(4.19) 


Because n j_ t y an< * ^L3 no ^ needed to derive the second 

order expansions of the perturbed force law (3-46) and Maxwell's equations 
(3*26)j these local expansions of v^, n 7 and P were not given 

in (4.12), (4. 16 ), and (4.19), The nonlocal expansions, v~, n^, and 

'■'O 3 

however, are included in (4.9)* (4.15), and ( 4 .l 8 ) because they will 
be used in deriving L^, 
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4,2*3 Expansion of the Lagrangian 

The expressions for the of (4,l) can be obtained by substitution 

of (4.8), (4.l4), (4.l8), and the second expression of (4.7) in (4.6). 
Then, for the unperturbed Lagrangian density, we have 



(4.20) 


where the term P/(y-l) is a result of generalization from ^ 

by replacing 5/3 by Y (Newcomb, 1962 ) . For the perturbation Lagrangian 
densities, we have after using integration by parts. 


£ rE[ nm vi + P 7 -?-i n (V ?-v<p-|.-a) + qn V(£i + i - ™,}] + s o£'£i " ^ ~ £1 * 
* 2 E Jl “ 5 2 - I [(Y-iX ’-?) 2 + + \ |£:w<p) 


2 2 

e„E. b. 


♦ + t Ii : ™) + i-(ii + l’™)] + -¥- - ir 0 S ' 

(v ’5> 3 + TT + I 

- qn|| + 1 5 |§|: WVcpj 

+ “■(;v(§ ll :TO £i + 5 + + 1 1£ ;VV *)] I • 


(4.21) 


The Euler-Lagrange equations of L, from £ of (4.20), with respect to 
x, A, and cp, are the unperturbed force law and Maxwell's equations. 
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ft 

mnv D + VP - qn( E + X B) = 0 , 

1 

^ vx r e oK ~E qn ^ = 0 * e o v ‘£'L qn = ° • ( 4 - 22 ) 

Use of these shows that L^ from £ of (4.21), is identically zero. 

The Euler- Lagrange equations of L ^ from £ of (4.21), form the 
system of linear equations in j, A , and cp , describing wave propagation. 
We shall not derive them at this point, but consider the Euler-Lagrange 
equations resulting from variation of L 2 + Ly describing nonlinear 
wave-wave interactions. 

Variation in 9^: Taking the variation of h n + L^ with respect 

to <p yields 

Vil + E q [ 7- ( n |) - \ 77; (”^)] - 0 • (*>.23) 

Local expansion of the second expression of (3.26), and use of the last 
expression of (4.22), gives 


VS Ll + "L2 ) = 


By use of (4.l6), it is easy to demonstrate that (4.23) and (4.24) are 
identical. 

Variation in A^: Taking the variation of with respect 

to A^ yields 

1 S ~1 

7 X ~1 ~ e 0 - 2Z q [ n i“ “ 7 *( n IS) + | =° • (^*25) 
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Local expansion of the first expression of (3.26), and use of the second 
expression of (4.24), gives 


— V X B - 

n 0 ~1 


Se. 
5 0 St 


1 -E^Ki 


n v + 
Ll-4) 


nv 


12 


n Ll^Ll 


n L3^0 


)■ 


0. (4.26) 


Use of (4.12) and (4.l6) shows that (4.25) and (4.26) are Identical. 

Variation in Taking the variation of with respect to 

£; yields 
*</ 

mn‘i- (y-l)VP V • §- V£ • VP+ yp|(y- 1)(V*5)9( V*§) + V(v§:VE;) + V£*V(V*§) 

+ 7p|(y- 1)(7«§) 2 + V§:V§ + (y-l)( V*£) V§*VP + VF*V£*VP 

- qnl E + § • VE + v X B - § • VB X v + ? X B ) 

IV'-'X ^ <-v4 } *<> ~ **Ju ^ f 

- •> n [s- ™i + |lS ; ™) x ^o-(£i + i-'®) x i]- 0 • < 4 - 2 V 

The local expansion of (3.46) at x, with 7 *p replaced by VP, 
does not agree with (4.27)- The discrepancy occurs because the former 
is inconsistent with the nonlocal expansion of L * of (4.3)- The 
appropriate expansion of the force law must be constrained by the particle 
conservation law, i.e. by comparing the first expression of (4.22) with 
c the force law 


" dv * n 

_ m '3r + v ' p ' - «” '<£' + Id' x £')J - 


0 


(4.28) 
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which describes the same cell of particles in its perturbed motion. 
The nonlocal expansion of (^.28) is 


mn| + VP V ♦ §+| (V -£) 2 - | (V£:V£)] 


- qn 


[i + ?•?+§ - 1 ( 5? i' 7 |.)] ( VP L1 + l' TOP+ '® L2 + ! , '’ w u + iiS :VWP ) 


(ii + r™ + £ , ’Si + lil ;TO £) - q ”(i >< £ + ^o x £i'S'™ x iD 


+ § X B - ?*VBX ? - §‘VB X V - i 5?: WB X v 

<*s> rv ~ rw> £_ 


)”° * 


(ll.29) 


where (4-. 8), (4,l4), (4.15), and the first expression of (4.22) have been 
used. Introduction of (4.19) establishes that (4.28) and (4.29) are 
identical. A similar force law expansion is examined in Appendix B. 

4.3 Nonlinear Wave Coupling Coefficients 

We shall now use to obtain the coupling coefficients for 

waves in a homogeneous, stationary, two-component magnetoplasma . The 
general expression will be derived by the averaged-Lagrangian technique. 

4.3.1 Averaged Lagrangian 


It is convenient to introduce the normalized quantities 



me e 

e 


O) 

^ = to 3 £ = to - £ 3 £ = £ 3 T = “pe* 3 ^*30) 

pe pe 
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where the subscript s denotes electrons or ions; electron 

plasma frequency, c is the free space speed of light, and is the 

Debye length normalized against c/a> . Then X^, and X, from (4.21) 

p€ C. J 

may be written as 

- °[l k 2 - r v i 2 < v, i> 2 ] + ii% - lii'Ui x + -* v ’ii 


+ I C . 2 - r - ie’ 5 + IVU. X - » 7 'ie 


1 • 2 1 t 2 

+ 4 v§ + a - ^ |v x fi , 


.*3 = ^ (ViiU^:^.) - £ (2-YjK^) 3 ] + A l -'S-i 1 -li l £ 1 :W§ 


v e 2 [r (Y-4)( v ie :7 A) -5 s ( 2 -Y e )(v-4> 3 ]- i e -'SAe + f AA, ;OT} ’ C*- 31 ) 


2 

where the values of and ha ve been normalized against n e m e c * 

the pressure terms have been reduced from their counterparts in (4.21) 
by integration by parts, and we have assumed the background charge to be 
neutral (q.n. = en * q = -e) and the dc VVA and Vcp to be zero. 

The averaged forms of £ p and £^ ma y obtained using the Fourier 
transformation 

c (X,T) = — ~ |dK fdn c (K,Q)exp i(fiT-K-X) , (4;32) 

'* e ~ (2it) J ~ ~ ~ 

where C (K,fj), K, and Cl can vary slowly with the normalized space 
and time coordinates, X and T, due to weak nonlinear wave- wave inter- 
action. If the scales of these slow variations, AX and AT, satisfy 
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ATj 0 -q'| » n , | AX’ ( K-K ') | » jt , ( 4 . 33 ) 

then we have (Schiff, 1970) 

J dX | dT exp i[(CWT)T- (K-K0*xJ^(2n ) 4 6(0-0') 6 (K-K') . (4.34) 

AX AT 

The action integral of (4.3l), 


I i = J‘“J dTf l (1 = 2,3) , 


(^.35) 


can be rewritten by use of (4.32) and (4.34). For i g , we have 

I = fdO fdK A< 2 > , A < 2 > . A i(2) + A e ^ 2 ) + a F ^ 2 ^ 

2 J J ~ K,n ' K,Q %Q \n + \ a 

= CT [i" Uil 2 -^ VViJ 2 ] + in 4‘5*-| Q ii-(4 x 4) + ^(K*ii) * 
A KyQ “ 2” lie I ~T J^ 2 ~ iQ 4 ‘if + | Q i e ' ( 4 X I^ ) " 


l 2S ) -5te»‘«S 2 -J* 2 lal 2 , 


2 1 ,, 2 1 „ 1 2 


(^• 36 ) 


where the Coulomb gauge (V-fl = 0 ), and the symmetry, C(-K,-Q) = C*( K ,Q), 
have been assumed. The latter reflects the requirement that t (x,T) 

be real, and has also been applied to $, and & . 

~x 

To reduce I 3 to its simplest form, it will be necessary to use the 
properties of linear waves, which can be obtained from I We shall 
therefore establish them before proceeding further. 
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4,3.2 Linear Waves 


By taking the variations of Ig with respect to ^(K^fi) and 
(1 (K,fi), we obtain the following linear equations 

a[n 2 £. - y.v 2 K(K-cJ|- i qa+ ik$ + ±nt. x n = o , 
fi 2 c - y v 2 k(k*c ) + ina - ik$ - i qc x n = o , 


K 2 tf + 0(K$ - fifl) - iQ(£. - C ) = 0 


(4.37) 


By introducing the normalized electric field, 

eE, 


e = i(K§ - fifi) = 


cd m c 
pe e 


we can obtain from (4.37) 
iK«e 


5 = 


K 


./ JfK’ftX 

1 I _ /NJTU /-V 1 

~ _n r k 2 / ' 


C = M -ft , 

~e 


(4.38) 


c I . (4.39) 

C7 


Here, the polarization tensor, Bfl , is defined in a Cartesian coordinate 


system with jQ^ in the positive z-direction. Its elements are given by 
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are defined 


where * 
as 



denotes the complex conjugate, and Q . Q . and A 

e* x' 


n i-r ’ 


(*.4i) 


Using (4,39) in (4.37) yields the familiar results (Allis, Buchsbaum, 
and Bers, 1963 ) 


£*e = 0 , d = nn - i + k , 


I 


- U - 

~ e 


M, 

a 



( 4 . 42 ) 


where Jt is the plasma equivalent permittivity tensor. Nontrivial solutions 
for (4.42) exist when K and 0 satisfy the dispersion relation, 

det D ^ ° . (4.43) 

For these values of Q(k), the electric field polarization vector is 


e ee 


(D d 
yy xz 


- D D ) 
xy yz' 


(D D - D D ) 
xx yz yx xz' 


(D D 
xy yx 


-D D )j 
xx yy | 


(4.44) 


We shall postpone discussion of the various linear waves described by (4.43) 
until Section 5, where Langmuir, ion-acoustic, and whistler waves will be 
treated in applications of the results to be derived in this section. 

^•3.3 Wave Coupling Coefficients 

We now return to the derivation of wave coupling coefficients from 

^3* We first separate into its component wave amplitudes and 

unit polarizations, 
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■ (^. 45 ) 


£(K,Q) - T )% > 

r^j 

u 

where £^(X, T) has slow X- and T-dependences because of nonlinear 
K v ~* ~ 

wave interaction. Use of (4.45), (4.32), (4.34), (4,35); and (4.39)* 
in (4.3l) then gives 


I 0 = i 1 dKdK 'dK" SZ't" k vv -„// , 

3 I ^ ^ Kh. ft 

s 

a kk'k' - E %[ A T {&£'■£") + A I (w'.i')] - 

/W /^w< f ' 

F, S 

A J(k. K K") . V*(k-M S -e)p5 (* 

( 2 -v s )(i''-“ s ’-i')(i5'-“’ 4')] - 


\ 

r 


a i(m'^') - ?f (•*' V) (S'S* V) 


where t . , t , e ", and £ " are given by 
1 6 ~ 


2 ^ U ^ ^ U 

Tj « I/O ; T e = -1 > e'=e R , > & = ^ ' 


K 




in I~, the integration, j, over K, K , and K is restricted by the 

3* - s ~ ~ 

synchronism conditions, 


K + K ' + k" = 0 , n( K) + Q '( K ') + cf ( K" ) = 0 

— • ■ ~ * Av ~ ' 


(4.46) 


(4.47) 


(4.48) 
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The summation in A * n is assumed to cover all permutations of its 

KK K 

j-vyv^ 

+ ff ” 

following terms with respect to K, K , and K . By use of the frequency 
synchronism condition of (4.48), we find that the second term of 
A^(K, K K y/ ), when summed over all permutations, will make zero contribu- 
tion to A rk , k „ . 

It is noted that k * n is composed of two types of nonlinear inter- 

KK K 

action energies. One is associated with the thermal motion of the charged 
particles, represented by A^( K, K K*) . The other, Af(K, K k") 


r ,. s " a/.. ,.s * ,\ / 


i = w e‘iw 


..s' A/.' 

;(*•* )/“ 


, is the field-particle interaction contri- 


r 

bution associated with the current, charge density, and electric and 

magnetic fields. It consists of those terms in of (4.3l) that are 

multiplied by the unperturbed charge density, q n . 

s s 

The quantity (A^ + c * c *)> with c.c. representing the complex 


conjugate, can be interpreted as the rate of energy transfer to the (-K) 
wave due to nonlinear interaction between the K * and K" waves of unit 
amplitude in the normalized electric field (Galloway, 1972). Because of 
this energy transfer, the energy of a single wave is no longer conserved, 
as in the linear case. Instead, we have (Galloway and Kim, 197l) 


s(!t \ + fx ‘ Ik) = J d *' d £" ee ' e " a kk V 


(4.49) 


where S signifies that the synchronism conditions of ( 4 . 48 ) are satisfied 

The wave energy and energy flux densities, 1L and T , can be obtained 

K ~K 

by applying the Legendre transformation to £ , but retaining the original 
dependent variables $, and (1 , before applying the averaging process 

They are given by, 
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with * denoting the complex conjugate. When only three discrete waves 
are involved, the right-hand side of (4.49) reduces to a summation symmetric 
in K, k ', and k" . We have, by substitution of (4.50), 

^(Ir + ~k ’ tx) | & 4 * V'(!t + V * lx) l e 1 2 

- 3 k"(It + V • lx) |«f - <*'«' ‘bY • f 4 -? 1 ) 

The first two equations are the Manley-Rowe relations expressing the 
conservation of wave action^ in three-wave interaction 

II 

( Louisell, i 960 ) . ( 

The coupled mode equation can be obtained from (4*51) by cancelling 
ft , ft , and ft in (4. 51 ), 
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3 U /|_+ u u . e u = _ c * *v w 

P \ ST ~p bX/ p psi si 9 

£ (fe + 3ii • fe) i i = ° P si ®r «p - 
>: (if ■ * £ • fe) < - v - 


(4.52) 


where U, V ? and W denote the waves, and the following substitutions have 
been assumed^ 


C 

psi 


r uvw ^ u'vw 

kkV - a -k,kV 




y 


K = 
~P 


-K 


y 


K 

~s 


K' 


y 


K. - K* 

rv->l /~v> 




U 







(4.53) 


with the last two relations following from (4.50) . The subscripts p, s, 
and i denote the pump, signal, and idler waves, respectively. This 
nomenclature is commonly used in describing parametric amplifiers, where 
a large amplitude pump wave, £^, is applied to amplify the signal and 
idler waves (Louisell, i 960 ). 


^•4 Parametric Wave Amplification 

In Section 4.3 we have derived the general coupled mode equations for 
nonlxnear three-wave interaction under three important assumptions: first, 
that the waves are undamped; second, that the background plasma is homo- 
geneous, and third, that there are no additional nonlinear wave processes, 
such as wave-particle interaction. We shall now consider spatial solutions 
to the coupled mode equations, and assess the effects of relaxing some of 
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these assumptions . In particular, some results with damping taken into 
account will be considered. They will be of use in Section 5 f° r com- 
paring theoretical predictions with experiments on parametric amplification. 

4,4.1 Parametric Amplification Without Wave Damping 

The solutions of the coupled mode equations in (4.52) are well-known, 
and will be briefly discussed in this subsection. We proceed with the 
assumption of a large pump wave amplitude. 

Is I » |S L |S. |: In this case. S can be considered as a constant 

ipi 's i P 

because the right-hand side of the first expression of (4.52) is now 

negligible. The remaining two equations become linear in S g and S^, 

Their solutions are then subject to the instability criteria well-known 

in the linear wave theory concerning absolute (temporal) and convective 

(spatial) growth (Sturrock, 1958b; Briggs, 19^4; Derfler, 1967 ) . This type 

of analysis has been extended by Harker and Crawford ( 1989 s ) anc * 

Van Hoven (l97l). In particular, they have shown that, (a) if 3 3, >0 

and > 0, & s and are convectively unstable with both 

temporal and spatial growth possible, (b) if >0 and < 0, 

they become absolutely unstable, and evanescent in space but amplified 

in time, (c) if 3 3, < 0 and U U >0, they are evanescent in time 
si si 

but amplified in space, and (d) if < 0 and ij^li < 0, they vary 

like beat waves and no persistent amplification in space or time is 
possible* Cases (a) and (c), in which the signal and idler wave group 
velocities are in the same direction (l^lK > 0), are usually termed 
co-flow or forward scatter, while Cases (b) and (d), in which the signal 
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and idler wave group velocities have opposite directions (ft ft < o) 
are usually termed contra-flow or back scatter. The distinction 
between the two types of parametric amplification is illustrated in 
figures 4.2(a) and 4.2(b). 

The temporal (assuming space-independence) and spatial (assuming 

time-independence) parametric amplification rates, IV and K , 

’ 0 O’ 

respectively, are given by (Louisell, i 960 ). 


r 0 - VW 


1/2 


l C p8ill«pl 


(4.54) 


These expressions can be obtained from the last two relations of (4.52) 

by assuming that £ g and £. vary as exp^T) or exp(K Q Z) in a one- 

dimensional approximation in the z-direction. Since imaginary r Q or 

K q amounts to a small nonlinear frequency or wavenumber shift, respectively, 

for the signal and idler, the results of (4.54) are seen to be consistent 

with those by Harker (1969) and Van Hoven (1971). 

If we assume that the initial (or boundary) value of £. is zero, 

(4.54) breaks down because £^ can no longer be an exponential function 

in T (or Z) . We then find from the right-hand side of the second expression 

of (4.52) that the initial slope of £ in T (or Z) is zero. Also, 

from the last expression of (4.52), the starting behavior of £ is seen 

1 

to be T 0 .Z (or K 0i T), with 


r oi = 


Voi 


3 i 


(4.55) 


When [£. j becomes comparable to |ej due to this linear growth, the 
exponential growth of (4.54) is then applicable. The spatial behavior 
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Figure 4,2 The distinction between (a) the co-flow (forward scatter) 
and (b) the contra-flow (back scatter) cases of parametric 
amplification. The solid lines indicate the dispersion 
curves near the (K, Q) of the pump, signal, and idler waves. 

In the co-flow case, the group velocities (ll ** dQ/dK.) of 
the signal and idler have the same sign, while in the 
contra-flow case, their group velocities have opposite signs* 
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of and 6 ^ in each of the above mentioned four cases were presented 

in a summary paper by Barnes (1964). 

| £ | ~ 1 £ | ~ | £. | : In the case of parametric growth in £ and 

P S 1 s 

£^, the results of (4.54) and (4.55) can be considered accurate only near 
the boundary where, or soon after the initial time when, the pump wave is 
being injected. These results break down as soon as |£ g | and |e.| 
become comparable to l^pl' this case, the temporal solutions to 

(4.52) must be put in terms of elliptic integrals involving the wave 
actions 3 Is | 2 , 3 1 [ft | 2 , and 3, \l. | 2 . The solutions of these wave 
actions, when all the 3 ’s have the same sign, are illustrated in 
figure 4.3 (see for example, Sagdeev and Galeev, 1969 ) . Several features 
of this figure should be noted: first, that jC^j can no longer be con- 
sidered as a constant; second, that the interaction process is reversible 
and has a nonlinear period T^, which can be shown to be roughly propor- 
tional to the inverse of the maximum value of \& I, and third, that the 

p ’ 

Manley-Rowe type of wave action conservation laws, derivable from (4.51), 

3 p^p | 2 " 3 S |£ S | 2 = constant , 3 g | | 2 + 1 t ± | 2 = constant , (4.56) 

are satisfied. Finally, it should be noted that if all the ll's 
have the same sign, the spatial solutions in this case are similar to the 
temporal solutions shown in figure 4.3. To complete the analogy, it is 
only necessary to replace T and the .l|(>| 2 's in figure 4.3 by Z 
and the corresponding 5U|&| 's, respectively. Also, (4.56) becomes, 

V^pl £ pl - = constant , V l) s ! S s| 2+ 3 i 'U i | | 2 = constant , (4.57) 

which is the Manley-Rowe type of power-balance formula (Penfield, i 960 ) . 
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T„ tor Z„) 


T (or Z) 


Figure 4. 


3 The solutions to (^.52) when |& p |> l&Jz and f I are 
comparable in magnitude are in terms of the wave action, 
for temporal behavior, and wave power transfer, 
3U|ft| 2 , for spatial behavior* (or Z^) is the 

interaction periodicity time (or spatial length) of 
wave interaction. 
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A situation in which the wave amplitudes increase to infinity in a 


finite time period occurs when 
from (4.52), 


3 < 0 and 
P 


3,5. > 0. We then have 
s' 1 


1 3 

n 


de 


p 1 dT 


= c e e. 

s 1 


de 

3 — 

s dT 


= c e.e 

1 p 



= c 


e 1 

s p 


( 4 . 58 ) 


where C = C . , Now & , £ , and £. can be made real in value by 

1 psi' p s 7 1 J 

choosing n/2 as the phase angle difference between the pump and the 
signal and idler waves in (4.52) (Davidson, 1972). From (4.58), we see 
that, whenever £ p , £ g , and simultaneously have the same sign, the 
wave amplitudes grow without bound. As illustrated in figure 4.4, we 
have for the special case £ p = £ g = (Davidson, 1972), 


£ (T) = & (0) 

P P 



c & (0) 

p 



T J 


(^. 59 ) 


which becomes infinite in a time period T = |3 j/C ft (0). This is the 

exp 1 p 1 v 

simplest example of the so-called explosive instability (see for example, 
Coppi, Rosenbluth, and Sudan, 1969 ) . 

This infinitely large amplitude cannot occur in practice, however, 
because of the limiting effects of linear wave damping, and other nonlinear 
wave processes, which are omitted in the foregoing model of a wave- 
triplet, On the other hand, the spatial solutions to (4*52), when 
3 'll < 0, and 3 U , 111 > 0, have a similar form to the temporal 

P P s s 11 

solutions shown in figure 4*4, with the explosive instability length 

being Z = I 3 li l/C ft (0). In the next subsection, it will be shown 
exp P P p 3 

that the inclusion of weak wave damping will drastically change the 

behavior of ft , ft, and ft from that briefly indicated above* 

P s x 


no 




Figure b.k The behavior of | & | in the explosive instability of 

P 

three-wave interaction, occurring when 3 is different 

in sign from 3 and 3, for their temporal solutions, 
s 1 

or 3 \\ is different in sign from 3 ll and 3.11. 
p p s s xi 

for their spatial solutions* 
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4.4.2 Wave Damping 


The effects of weak linear wave damping can be included phenomenolo- 
gically in the coupled mode equations (4.52) (Sjolund and Stenflo, 1967 ), 


1,(4= + K ‘ 4? + r = -c* . e e. , 

p\ oT ~p ax p / p psi s 1 9 

1(4? + , u*4r+r)£=c , e*e , 

s\oT ~s oX s/ s psi 1 p 9 


JL. 

“i\3T 


u . 




a 

ax 


T 

A i / 


0 

ly . 
1 


psi 


* 

o 

s p 


M. r 




where the r represent the normalized damping rates of the pump, 

signal, and idler. Due to the presence of T . , |fl | must be larger 
than some threshold before £ and e. can be amplified. This threshold 

S 1 

field, C for amplification can be obtained by setting d/dT = 0 = 3 /Bx in 
the last two expressions of (4.60) . We then obtain 


£ 


th 


< s . 3 i r . r i> 


psi 


(4.6i) 


after cancelling £ g and £^ . When the pump field is not far above the 
threshold, we may obtain the initial parametric amplification rate for 
waves propagating parallel to Z, in a fashion similar to that used in 
obtaining (4.54) (Louisell, i 960 ), as 


2 2 
1C 2 = K 0 K th 

V^i 


(4.62) 


where K. Q is given by (4.54); ja g , 1 ^, and are defined by 
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r s r i 

= T ’ * TT ' 

s i 

2 2 

and we have assumed (u + uij » 4( K. - K ). 

x s i' ' 0 th' 

For the case of a large amplitude undamped pump wave, the spatial 

solutions for linearly damped signal and idler waves have already been 

presented in the context of nonlinear optics ( Bloembergen, 1965 ; Bobroff, 

1965 ). The stability analysis, in line with the approach by Harker and 

Crawford (1969a) and Van Hoven (l 97 l), has been extended to this case by 

Bers, Chambers, and Hawryluk (1973). Here, we wish to extend the analysis 

to account for a damped pump wave. A similar problem has been studied 

by Ohnuma and Hatta (1970), and by Porkolab and Chang ( 1970) . 

By making the phase of t different from those of t and t by 

p s i 

rt/2 in (4 P 6o), we have for the one-dimensional case (Davidson, 1972), 




ce.e 

1 p 



c t l 

s p 


(4.64) 


where C = |^p S ^| ? an d are real in value, and we assume an 

exponentially decaying, large amplitude pump wave 


£ p = £ p (o)ex P (-n p z) (|e p ) » |e g | ; |e i |) . (4.65) 

Case 1 . u 4 u = n = |i: The simplified case where u. = u. = li = u 
P s 1 p 'g “ 

has been treated by Ohnuma and Hatta (1970), and a further simplification 
to second harmonic generation by a large amplitude ion~acoustic wave has 
been studied by them and by Litzenberger , et al. ( 1972) . For purposes 
of comparison between theory and experiment in Section 5.2.3, we shall 
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consider solutions for parametric amplification when p ^ p = u = u 

P s p 

We then have as solutions of (4.64) , for the case of co-flow (figure 4.2a) 


= i 


ft = ft (0) cosh M + ft. 


/V 

(V2 -| 

(k 3 

\ s s 

J sinh M 


r fu 3 \V2 1 

\ = ^8^(0) cosh M + 8 s (0) ^~ ~j sinh MJ exp(-pZ) , 


, v r 1 c e (o) e (o) 

M(Z)=M 0 [l-exp(-nZ)J , M 0=— ( W V2 ’ < 4 - 66 > 

^p lll S J s ll iV th P 


where ft g (0) and ^(0) are the values of £ g and ft^ at Z = 0, and 
the second expression for M Q is obtained by using (4,6l). 

By evaluating d£ /dZ and dft./dZ at Z = 0, we find that ft 

» 1 g 

and ft. grow in the vicinity of Z = 0 according to (4.6l) and (4.62). 

For Z such that P p Z » 1, ft^ and ft^ both decay by linear damping. 

This behavior is illustrated in figure 4.5 and is drastically different 
from that shown in figure 4.3 where damping was not included. In the simpler 
theoretical model used by Ohnuma and Hatta (1970), in which p p = p, the 
following assumptions were made, 

. . 8 (0) /3 \l/2 

M(z) « 1 , e s (°) > -Jj ( y ) ’ (4.67) 

0 \ s / 


where the second condition is obtained by making dftVdZ >0 at Z = 0 
according to the second expression of (4.66). It then follows that ft 

i 

has its maximum value at Z Q = -tn2/^ - o.7/n pJ a result obtained 
without assuming p = p. 


il4 




Figure 4.5 The schematic behavior of £ and £ according 

s i 

to (4.66), where wave damping of and 

has been included. When ll = 4 . and under the 

p 7 

condition ( 4 . 67 ), the peak location of 

becomes independent of & (0). 
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Case 2. p 4 M- 4 M<. - The solutions for £ and £. must now be 
p s i si 

expressed in terms of Bessel functions (Porkolab and Change 1970)* We 

shall consider only the case where \i > 0 > 0), u < 0 ('ll < 0), 

p p s s 

and > 0 (*lL > 0), i.e, the contra-flow or back scatter case (Barnes, 
1964; Bloembergen, 1965 )* as illustrated in figure 4.2b, relevant to 

examples in Sections 5*3.2, 5* ^*3* and 5 . 5-^. Under these conditions the 

solutions of (4.64) for £ and £. can be written as 
' s' 1 
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where J (x) and N (x) are Bessel functions of the first and second 
v v 

kinds (Abramovitz and Stegun, 1965 ) • 

Simplifications of (4.68) are possible under various conditions of 
interest for Section 5: 


X « 1: In this case^ &. and reduce to 

I T i 1-olSil l M ol (\\\K\ 1/2 r 1 

i £ i ( °i es(o) ^I'w) sinh L (v_i vj 

X exp [-(n p + n s + n t ) fj , 

...X M 2 1 KI/MVT ) 


' £ i (0) l^-TM^TTl + e s (0) 


s <°> 


x exp -in + u s + nJ 


a s (0) \s-iwnp 


l» 0 l ] l M 0 l 

^wrS hi0) — 

;( 1J p + "s + “i) l] . > 


v \ Vi 


sinh( vH Z) 


where j and are defined by 

lx Is 


X 1± = exp£(v-l)n p z] , X ls - exp(vn p Z) 


v » X: We now have for t. and £_ the expressions 

x s 


£,(0)/ \v-V2 i , v 

(^l) ’ & s^*2ij ex P(^i Z ) 


i e 


+ - e i (0) 


w-3/21 p 

/ \ 1 

1 expl- 

K + A z \ 


s e 


PP {pi ) V ^ e i (0)X 2 S j ‘A-**) 

_-4“~ - e s (0) (v) 7 exp [-(^ P + w i) z ] ’ 
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where X 2i and X 2g 


are defined by 


X 2 i " 


1l 3 

S ’ S 


v 2 s “ 2 v 1 i 3 
\ i i 


i\i » s \ v2 


and e has been written for exp(l) = 2 . 7 l 8 ... . 


X > v » 1 : In this approximation, £. and £ reduce to 

1 s 


e, - &, (O)cos M + £ ( 0 ) (!^\ ain m! *vX/„ + „ \ 2 

A L A s \ Vi / J ’“"L r s ' "V 2J ' 

/|1I |3 \!/2 1 t , -I 

e s ^ £ s (°) cos m + sin m ex p -|ii s + ui i j | 


III 1 3 \l/2 


£ s ^ [ £ s (0)cos M + Sln M ex ^-K + ^i) I 


We note that. 


M(Z) = l- 0 l[l - exp(- Mp Z)j - |M 0 | 4p Z , (4.75) 

where the second expression for M(z) is appropriate when u Z « l i n 

P 

(4.7°), (4.72), and (4.75), the quantity |m q | has been defined in (4.6g). 
It will be seen in Section 5 that the behavior of £ g and £. according 
to (4.70)-(4.75) is drastically different from that indicated in 
Section 4.4.1 where linear damping is neglected. 

Case 3. Noncollinear Propagation : So far we have assumed that 
the group velocities are collinear. In preparation for Section 5 . 5 . 4 , 
we now consider the case when (U p , 1 ^, and ]L are no longer collinear. 
The coupled mode equations (4.64) become 


118 




with 7 and Y. being the angles of K and li. with respect to the 

s i — s ~i 

positive Z direction. Substituting 

e s (X,Z) = e s (Z) exp(-p sx X) , & ± (X,Z) = e t (z) exp(-n. x X) , 

r s r s 

'Vr sl “' |, s * >win cosY s 

s s 

(^. 76 ) can be reduced to 

v„(!z + ^sz) . C ? t (Z) e p (Z) exp[( Msx - P lx )x] , 

Vlz(?E + •*!») V Z) - c £ s< Z) £ p< Z > “ffKx ‘ 0*] • ( “' 79) 

It will be seen that, along a path with X constant, (4.79) equivalent 
to (4.64), but with different wave coupling coefficients. Along an 
appropriately chosen Z-axis, where X « 0, e.g. the central plane of a 
plasma slab, the form of (4.79) becomes identical to that of (4.64). The 
solutions, (4„68)-(4.75), become the solutions of (4.79) when lJ s , ll , p , 
and p i are replaced by U gz , U iz , p gz , \i ^ 2 , respectively. 
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4.4.3 Effects of Other Nonlinear Processes and Plasma Inhomogenei ty 


In addition to linear wave damping, other nonlinear wave processes 
and plasma background inhomogeneity may also significantly alter the 
behavior of the interacting wave- triplet . Here we shall only briefly 
discuss the conditions under which the solutions of and given 

in Section 4.4.2 are still acceptable in practice. 

Other Nonlinear Processes : The competing nonlinear wave processes 

likely to occur in the presence of a large amplitude pump wave are 

primarily in the form of modifications of the background plasma. While 

causing them, the pump wave may have a significantly different behavior 

from the exponential decay assumed in (4.65), thus rendering the solutions 

for 6 and £. given in Section 4.4.2 inaccurate, 
s 1 

(i) Particles trapped by an electrostatic pump wave : In a collisionless 

2 2 

plasma, and assuming K ' « 1, this process is important under the 
condition (see for example, Davidson, 1972), 


% » % » r L ’ 



(4.80) 


where is the electron bounce frequency when trapped in a trough of 

the pump wave, and T is the normalized linear Landau damping rate, 

Li 



(4.8i) 


It is well known that, under the condition (4.80), the electrostatic wave 
decays like exp(~r L T) only for T < 1/Qg, but executes modulational 
oscillation when T > I/O . This behavior is illustrated in figure 4.6. 
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Figure 4.6 The behavior of an electrostatic pump wave amplitude under 

the influence of trapped electrons ( — — ) in a collisionless 
plasma under condition ( 4 . 80 ). In the absence of electron- 
trapping, the wave decays according to linear Landau 
damping as exp( ~r T) ( ). 

±j 
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The results of Section 4.4.2 are therefore appropriate only for T ^ i/q^ 
or for Z ^ 1/f^ll when considering spatial evolution. 

The first condition in (4.80) is necessary to restrict the size of 
I so that the electron density perturbation is much smaller than the 
quiescent density (Davidson, 1972), and is consistent with the assumption 
of weak perturbation used in Section 4.2 in approximating the Lagrangian. 
The second condition in (4.80) is required to allow the trapped electrons 
to make many bounces in a wave trough before they become un trapped due to 


wave danrninp-. 


Therefore whpn +h*=> 

—J w wxxxuxunui J-Xlg 


is much large r- 


than r c should replace P^ in (4.80). Following the same reasoning, 

T c should be replaced by the total measured pump wave damping rate T 

P 

which may be a sum of T c and the nonlinear damping rate that 

accounts for all the nonlinear effects on the pump wave. We then have 
from (4. 80), 


"b = ( K l e P l) 


(4.82) 


as the condition under which the trapped particle effects on an electro- 
static pump wave are negligible. 

(ii) Plasma heating : Substantial background modification in the 

form of plasma heating by an incident wave can occur in a collisional 
plasma due to collisional randomization of the ordered perturbation in 
electron motion. It has been shown that the electron temperature* T * 
increases from its quiescent value* T Q * roughly according to 
( Ginzburg, 1970)* 


= 1 + 




V e Q ' 


(4.83) 


122 



2 2 2 

when Q » v cd * with v being the electron-ion and electron-neutral 
e pe' e ° 

2 2 

energy transfer collisional frequency. When |fi^| » it is seen 

that the corresponding temperature increase can be substantial, thus 
changing the wave coupling coefficients and linear wave properties. 

Thus, in the presence of a large amplitude pump wave, the electron 
temperature should be determined from (4.83) before using the results of 
Section 4.4.2. 

(iii) Density Modifications : Particle trapping and plasma heating 
deal with the modification of the plasma background in velocity space. 
Modifications in configuration space can occur when the large amplitude 

wave is inhomogeneous in amplitude, as is the case in (4.65). This results 

2 

in an equivalent potential energy, (normalized to ...me ), for each 

charged particle (Ginzburg, 1970), 


|e I 
$ = LsL 

N o 2 


l<y o) r 

p exp(-2u Z) 

CT P 


(4.84) 


This nonlinear effect becomes important when is comparable with the 

2 

average electron thermal energy, /2. So the condition under which 
the background plasma density is not significantly altered by an inhomo- 
geneous large amplitude pump wave is then. 


2 J2 


ft f « v ri 

p ' e 


(4.85) 


over the plasma region where j-i^Z ~ 1. The effects of background 
inhomogeneity on the three-wave interaction process will be briefly 
discussed below. 
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Plasma Inhomogeneity : Hie size of the plasma region, AL^, in which 

nonlinear wave interaction occurs, may be limited either by boundaries 
in the case of a uniform plasma, or by satisfying the synchronism condi- 
tions (4.48) in the case of a slightly inhomogeneous plasma. 

(i) Bounded homogeneous plasma : In this case, we require 

KAL q » 1 , KAL q > 1 , (4.86) 

to see significant nonlinear wave amplification within AL^. If there 
is a departure AK (= K-K — K ) from synchronism, then the parametric 
growth will be substantially unaffected only if 

AK • AL q < it . (4.87) 

This conclusion may be reached by considering the exact expression for 
the integral in (4.34) (Schiff, 1970; Phelps, Van Hoven, and Rynn, 1973 ) 

f 0 / AK * AL n\ 

l L dZ ex P<> AK * Z > - ak Sin (^SI^j > (*.88) 

which has the significant value of roughly unity only when (4.87) is 
satisfied. Note here that we have neglected the reflection of waves from 
the boundaries. 

(ii) Weakly inhomogeneous background : In this case, the quantity 
AK is a function of AL, where AL is the distance from the point of 
perfect, synchronism. Then (4.87) can be used to determine the width, 

AL q > of the region over which synchronism is sufficiently well satisfied 
for parametric amplification to occur. Assuming linear dependence of AK 
on AL, we then have from (4.87) 
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AK = aAL 


0 *. 89 ) 


4L 0 < 





which agrees with a result by Harker and Crawford (1970). 


^•5 Discussion 

The purpose of this section has been to apply the macroscopic - 

Lagrangian obtained in Section 3 to the description of parametric wave 

amplification phenomena in a homogeneous plasma. In Section 4.2, the 

Lagrangian of ..(3.22) was expanded in terms of the perturbations, S, A, 

and Cp , in plasma cell position x and the vector and scalar potentials 
1 ~ 

A and In Section 4, 3, the averaged Lagrangian, technique was applied 

to obtain the energy densities and energy flux densities of the linear 
waves, and the nonlinear wave coupling coefficients. We specialized in 
Section 4.4 the analysis of the three-wave coupled mode equations to 
parametric amplification, and demonstrated how to solve them in the presence 
of linear wave damping. The results will be used in Section 5 to make 
comparisons between theory and experiment. 
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5. THREE -WAVE INTERACTIONS: APPLICATIONS 


5 « 1 Introduction 

Nonlinear wave-wave interaction represents the lowest order correc- 
tion to linear plasma wave theory , It must be described and understood, 
together with the phenomena of nonlinear wave-particle interactions, which 
have been ignored in this thesis, before effective progress can be made 
towards understanding such significant plasma problems as turbulence and 
anomalous transport (Kadomtsev, 19^5) * In view of its importance, one 
would expect to find a variety of theoretical predictions for interactions 
between various types of plasma waves, and a corresponding series of 
experimental investigations designed to test the validity of the theory. 

This desirable stage has not been reached, however. Although there 
have been many analyses of wave-wave interactions, there have been few 
experiments carried out to test them rigorously. A representative sample 
of the literature is given in Table 5*1. The list is not intended to be 
exhaustive, but illustrative. Most of the nonlinear wave coupling 
coefficients have been derived by the conventional iterative method. This 
approach is conceptually simple, but generally involves greater algebraic 
complexity than the Lagrangian method. 

Although all of those results listed could be rederived by the 
averaged Lagrangian method, by use of an appropriate Lagrangian density, 
only a small number of the analyses were carried out by this method. 

These may be categorized by the plasma models used. 

Cold Plasma : Using the cold plasma Lagrangian, Galloway and Crawford 

(1970) have obtained the coupling coefficients for transverse waves 
propagating at arbitrary angles to the static magnetic field. As men- 
tioned in Section ^J-.l, a similar study has been presented by Boyd and 
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Table 5.1 Some derivations of wave coupling coefficients using (a) the 
averaged Lagrangian method and (b) the iterative method, with 
(c) related experimental work. The abbreviations used below 
represent, C; cold plasma, M: macroscopic warm plasma, and 
m: microscopic warm plasma. With respect to the static 
magnetic field, we have, || : parallel propagation, -ff- : quasi- 
parallel propagation, Jl: perpendicular propagation, and 
X: oblique propagation. With respect to the directions of the 
phase velocities, we have, ® - 0: col linear waves, and 0^0: 
noncollinear waves. 


Interacting Waves * 


References; Theoretical Assumptions 


Three ion-acoustic (b) Litzenberger and Bekefi ( 1969 ); M, 0=0 

waves 

(b)(c) Ohnuma and Hatta (1970); M, 0=0 

' ( b) ( c) Litzenberger, Mix, and Bekefi (1972); 
M, © = 0 


Two electron plasma 

(a) 

Suramlishvili (1964); m, © 4 0 

waves and one ion- 
acoustic wave 

(a) 

Kim (1972) ; M, ® = 0 


(b) 

Oraevskii and Sagdeev ( 1963 ); M, © 4 0 


(b) 

Lee and Su ( 1966 ); M_, © ^ 0 


(b) 

Gratzl (1971) ; m, © 4 0 

Three electron plasma 

(a) 

Galloway and Kim ( 1971 ); m, © = 0 

waves 

(a) 

Kim ( 1972) ; m, -ff - } © ^ 0 


Two whistlers and one 
ion-acoustic wave 


(b) Forslund, Kindel, and Lindman (1972); 

M and m, |j ■ 

(b) Porkolab (1972); m, X, 0^0 

(c) Porkolab, Arunasalam, and Ellis (1972) 
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fable 5,1 (cont.) 


Interacting Waves References; Theoretical Assumptions 


One long wavelength 

(a) 

Kim ( 1972 ); M 

ordinary wave, one 
electron plasma wave, 

(b) 

Nishikawa ( 1968) ; M 

and one ion-acoustic 

(b) 

Lee and Su (1966); M 

wave 

(b) 

Harker (1971); M 


(b) 

Dubois and Goldman (1965); 


(b) 

Silin (1965); m 


( c ) 

Stern and Tzoar ( 1966b) 

One extraordinary 

(b) 

Porkolab ( 1972 ); m, X, 


wave, one elec- 
tron plasma wave, 
and one ion-acoustic 
wave 


One cyclotron harmonic (b) Gratzl (l 97 l); in, ©^ 0 
wave, one or two 
electron plasma waves, 
and one ion-acoustic 
wave 


Two ion cyclotron* or 
two Alfven waves, 
and one ion-acoustic 
wave 


(a) Suramlishvili (1970); m, X , ©^ 0 

(b) Lee and Kaw (1972); M, m, || 

(b) Hollweg ( 1971 ); M, || 

(c) Belcher and Davis (l 97 l) 

(c) Dubuvoi and Fedyakov ( 1968) 


Two (or one) lower 
hybrid waves and 
one ion cyclotron, or 
magnetosonic, wave (or 
two ion-acoustic 
waves) 


(a) 

(b) 

(b) 

(b) 


(e) 


Suramlishvili (1971); m, X, 0^0 

Kindel, Oku da, and Dawson ( 1972 ) ; 

m, X, 0^0 

Fidone ( 1973 ); m, X., 0^0 

Karney, Bers, and Kulp ( 1973 ); M, X, © 4 0 

Hooke and Bernarbei ( 1972 ) 
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Table 5.1 (cont.) 


Interacting Waves 


References; Theoretical Assumptions 


Two cyclotron harmonic 
waves and one ion- 
acoustic wave 


(b) Tzoar ( 1969 ); m, X , ©^ 0 

(c) Chang, Porkolab, and Grek (1972) 
(c) Keen and Fletcher (1971) 


Three electromagnetic or 
right-handed circularly 
polarized waves 


( a ) 

(a) 

(a) 


(b) 


Galloway and Crawford (1970); C, X., ® / 0 
Boyd and Turner (1972b); C, X / 0 / 0 
Kim (1972); m, -fK ® / 0 
Harker and Crawford ( 1969 b); C, -fH ® / 0 


One transverse wave (a) 

and two longitudinal , v 

(electron plasma and ' ' 

ion-acoustic) waves 


Suramlishvili ( 1965 ); m, ® / 0 
Kim (1965) ; m, -H-r, ®/0 


Two whistlers and one 
Alfven wave 


(a) Suramlishvili ( 1967 ); m, X, 0/0 
(a) Harker, et al (1974); C, +-, 0/0 


Two circularly polarized (a) 
waves and one electron / % 
plasma wave ^ '■ 

(b) 

(b) 

(b) 

(b) 

. (c) 


Kim (1972); m, -H-, 0/0 

Boyd and Turner ( 1972 b); m, 0 = 0 

Harker and Crawford (1970); M, || 

Sjolund and Stenflo ( 196 ?) ; M, II 

Montgomery ( 1965 ); M, || 

Kim, Harker, and Crawford ( 197 1 ) ; 
m, "4b y ® / 0 

Stern and Tzoar ( 1966 a) 


Two ordinary waves 
and one extra- 
ordinary wave 


(b) Etievant, Ossakow, Ozizmir, and Su ( 1968 ) ; 

C, x, 0/0 

(c) Cano, Fidone, and Zanfagna (l9?l) 


Two electromagnetic (b) Boyd and Turner (1972a); m, x, 0/0 

waves and one 
cyclotron harmonic 
wave, or upper 
hybrid wave 
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Table 5.1 (cont.) 


Interacting; Waves 

Three cyclotron 
harmonic waves 


References; Theoretical Assumptions 

(a) Kim (1972); ® ^ 0 

(b) Harker and Crawford (1968); © 4 0 

(c) Porkolab and Chang (1970) 

(c) Chang and Porkolab (1970) 
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Turner (1972b). The excitation of Alfven waves by nonlinear interaction 

of whistlers in cold plasmas has been studied by Harker, Crawford, and 

1 

Fraser-Smith (197^). 

Microscopic warm plasma : Using the Low Lagrangiaii (l958)> 

Suramlishvili ( 196 ^; 19^5; 19^7 ; 1970; 1971) has presented a series of 
studies concerning the nonlinear interactions between electron plasma 
and ion-acoustic waves, between one transverse and two longitudinal 
(electron plasma and ion-acoustic) waves, between Alfven and whistler 
waves, between Alfven and ion-acoustic waves, and among magnetosonic 
waves. Galloway and Kim (l97l) have obtained the coupling coefficient 
for three collinear, longitudinal electron plasma waves in this model. 
Boyd and Turner (1972b), in a study mentioned above, have considered 
the interactions between two transverse waves and one electron plasma 
wave in a warm plasma. A comprehensive study by Kim (1972) has presented 
the coupling coefficients, in the microscopic warm electron plasma model, 
for interactions among three plasma waves, among one plasma wave and two 
circularly polarized waves, among three circularly polarized waves, and 
among two plasma waves and one circularly polarized wave, all of which 
propagate quasiparallel to the static magnetic field. For the case 
of perpendicular propagation, he has obtained the coupling coefficients 
for interactions among three longitudinal cyclotron harmonic waves, and 
among one longitudinal and two ordinary cyclotron harmonic waves. 

Macroscopic warm plasma : In the work just mentioned,* Kim (1972) 

has also examined wave interaction processes using the macroscopic 
Lagrangian of Section ^-,2. His results include the coupling coefficients 
for interactions among electron plasma and ion-acoustic waves with 
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parallel propagation, and among an ordinary wave with perpendicular pro- 
pagation, an electron plasma wave with quasiparallel propagation, and 
an ion-acoustic wave with quasiparallel propagation. 

Most theories of wave-wave interaction have been carried out assuming 
the plasma to be homogeneous. In the experiments with which they have 
been compared, the plasmas have often been strongly inhomogeneous, 
rendering the wavenumber synchronism condition ( 4 * 48 ) difficult to 
realize experimentally* In these cases, comparison between theory and 
experiment will be impossible, unless the theory is improved by taking 
background inhomogeneity into account (Kino, i 960 ; Larsen, 1972). In 
mitigation, it should be pointed out that many of the experiments were 
aimed at studying the anomalous absorption of large amplitude waves, and 
the subsequent heating of the plasma (Porkolab, et al., 1972; Stern and 
Tzoar, 1966 b; Dubuvoi and Fedyakov, 1968 ; Hooke and Bernarbei, 1972; Chang, 
et al*, 1972; Keen and Fletcher, 1971), rather than to produce ideal 
conditions for testing basic plasma theory. 

Some experiments have, however, been carried out in effectively 
homogeneous plasmas with synchronism in both frequency and wavenumber 
realized. In particular, second harmonic generation caused by an ion- 
acoustic wave propagating along a plasma column have been observed by 
Li tzenberger, ot aj,( 1972) . Three interacting ion-acoustic waves have 
been studied by Ohnuma and Hatta ( 1970 ). Three interacting cyclotron 
harmonic waves have been studied by Chang and Porkolab (1970). Parametric 
excitation of ion-acoustic waves by whistler waves has been observed 
by Porkolab, et al .( 1972) * To explain the observed wave behavior, it 
was found necessary to include the effects of wave damping. With the 
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exception of the case for harmonic generation by ion-acoustic waves 
(Litzenberger et al., 1972 ), only qualitative agreement between the 
solutions of the coupled mode equation and the waves observed has been 
achieved (Ohnuma and Hatta, 197^; Chang and Porkolab, 197^? Porkolab, 
et al 1972 ). The purpose of this section, then, is to establish 
more nearly complete quantitative comparisons between some of these 
experiments and the theory developed in Section 4. 

In this section we shall specialize the general results of Section 4 
to two experimental situations. The first will be that of Ohnuma and 
Hatta (1970), for which we shall consider the collinear excitation of an 
r ion-acoustic wave by two other ion-acoustic waves (Section 5*2), and by 
two longitudinal electron plasma waves (Section 5*3). The second will 
be that of Porkolab, et al. (1972), for which we shall consider the exci- 
tation of a collinear (Section 5.4), and a noncollinear (Section 5 , 5 ), ion- 
acoustic wave by a large amplitude whistler, 

5,2 Nonlinear Interaction of I on -Acoustic Waves 
5,2.1 The Experiment and Interpretation 

The experimental set-up used by Ohnuma and Hatta ( 197 O) is shown 
schematically in figure 5*1* Their plasma had the typical parameters: 

n =" n. «“ 10 9 /cra 3 . T =“ 5 X 10^ °K , T, =* 300° K , 

e i e i 

7 ^ 

P *“'9 in Tor r , v ** ^.6 x 10 /sec , v. 3.0 X lO'Vsec , (5.l) 

y 7 en 'in 

where the effective electron- and ion-neutral momentum transfer collision 
frequencies have been calculated from (Ginzburg, 1970) 


133 




Figure 5.1 The argon (ct =“ 7*^ X 10 ) plasma column used by Ohnuma 
and Hatta ( 197 O). The pump (K^Q) and signal (k'‘,Q‘‘) 


ion-acoustic waves were excited with grid G . These 


waves and the idler (k } Cl") ion-acoustic wave were 


detected with grid G 


2 ‘ 






f 


en 


8.3 X lO'' (na 2 ) n T 1//2 sec" 1 (T in °K) 

c U c 0 


v. 1.1 X 10 5 (rra 2 ) 

in 


/T. \l/2 

i "him” ) ^c' 1 (T in °k) , (5-2) 

\ i / 1 


with n (cuT ) being the neutral density and m^(gm) the ion mass. 

The averaged electron-neutral collision cross-section, ( Jta ) “* 

8.5 X 10 -1 ^ cm 2 , has been obtained from the data given by Brown ( 1966 ), 

2 «- 1 ^ 2 

and the ion-neutral collision cross-section, ( jta ) . ^ 3.4 X 10 J cm , 

2 

has been assumed to be roughly 4(rta ) . 

For nonlinear interaction among three ion-acoustic waves, a large 

amplitude pump, £ (K,Q) and a small amplitude signal £ were 

7 p ~ 7 7 s 7 

excited by a grid immersed in the plasma. These waves and the parame- 
trically amplified idler, were then measured by a moveable grid 

(figure 5.1) . A typical example of the measured results is reproduced 

in figure 5 * 2 , which shows that £ and £ (when sufficiently weak) 

s p 

decay exponentially, but that £. is first amplified substantially from 

x t 

noise before it decays. The measured decay rate of £ g was found to be 

roughly the linear ion-acoustic wave damping rate due to ion-neutral 

collisions. But the decay rate of £^ for z ^ cm was found to 

increase with the exciting voltage, V , applied to the grid when 

V > V ~ 2V. This observed behavior of £ is reproduced in figure 5*3. 
p th p * 

and may be attributable to a combined result of wave-particle, and other 
wave-wave, interactions on the large amplitude pump wave, £^. 

Interpretation of the observed evolution of £± : It was found in 

the experiment that (~ 2 cm) is roughly independent of V , and 

P 

hence also independent of £ (0) if we can assume £ (0) cc V . To 

P P v ' P 

explain this observed feature of £ ± , Ohnuma and Hatta ( 1970) used a 
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Figure 5 . 


6 27 4 6 8 

z (cm) 


2 A typical observation by Ohnuma and Hatta (1970) of inter- 
action among collinear ion-acoustic waves. The exciting 

grid for the pump, & , and signal, & , is located at 

P s 

z = 0. The peak of the idler, is located at 
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set of coupled mode equations with identical wave damping rates, i a e. 

u = \i = u |_i. Assuming (4.67), they were able to show that 
p s i 

Zq (*“ 0.7/Hp = 0.7/M-) is independent of £ (0), as was shown in the 

paragraph containing (4.67) in Section 4.4.2. 

As indicated in figure 5-3, however, |i increases substantially 

with £ (0) [assuming £ (0) « V ] . Therefore, the theoretical behavior 
P P P 

of £ is more suitably described by assuming \i ^ p = p. = p in 
* p s x 

the coupled mode equations (4.64), whose solutions for this case are given 
by (4.66), Furthermore, because £ p is damped nonlinear ly, it is 
reasonable to assume that £ p (0) > Ihen, f rom (4.66), we have 


M(Z q ) 


e (o)u 

= — 2 - 


0 


& th |a p 


(5.3) 


showing that the first condition of (4.67) fails. Therefore, the apparent 
independence of in £ p (0) was not satisfactorily explained by the 

formula Z^ — 0.7/^. 

Since we can only assume £ (0) » £.(0), use of (4.66) then shows 

5 1 

that Zq should satisfy 

t “ h | M o[ 1 - exp( ' M p Z 0 ) ]j • (5 - 4) 


To determine Zq from (5.4), it is necessary to evaluate |_t, and 
evaluate M Q and |a p as functions of £ p (0)/£ th . This calls for the 
examination of the linear properties and the nonlinear coupling coeffi- 
cient of the ion-acoustic waves. 
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5.2.2 I on -Acoustic Waves (i) and Coupling Coefficient: 6=0 


In this case we have K = K i , where i is the unit vector in 
the positive z direction, along which the static magnetic field is 
directed. The polarization constant, M , of (4.40) then reduces to 


s s 

M = M 
xx yy 


_2 _ 2 
a -n 

s 


M 

zz 


D 


s s* s s* 

M = M = M = M 
xz zx yz zy 


= 0 


s s* 
M = M 
xy 


-in 


n(ct-n/) 


2 2 2 
D = fir - V V ■ K 

S s s z 


(5.5) 


For the experiment, we can assume 

v 2 » ctV. 2 (t » t.) , v 2 k 2 « l , (5.6) 

e, x e x 7 ez 

__Q 

where the second condition follows because ™ 2.9 X 10 from (5*l)> 

and K ~ 50 for a typical experiment wavelength of 0.21 cm. 
z 

Under the assumption of a low frequency longitudinal wave. 


e 1 = {0,0,1} , Cl « 1 , (5.7) 

■ 

the approximate ion-acoustic wave (i) dispersion relation can be obtained 
by using (5-5) a nd (5*6) in ( 1 4-.42) and (^.^3). For v e K z « i-t is 

0 2 ~ v 2 k 2 , v 2 ^-U 2 , (5.8) 

a z 7 aaee 7 

where V is the ion-acoustic speed. As illustrated in figure 5*^-j the 
a 

wave coupling synchronism conditions are automatically satisfied over the 
linear portion represented by (5-8). 
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Figure The dispersion curve for ion-acoustic waves (i) propagating 

parallel to the dc magnetic field under the conditions of 
(5.6). The synchronism conditions for nonlinear coupling 
of three ion-acoustic waves (p, and i) are always 
satisfied over its linear portion. 


lkO 



The corresponding ion-acoustic wave action, group velocity, and 
coupling coefficient can be obtained by using (5*5)"(5-8) in (4.46), 
(4.50), and (4.53), yielding 
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C . 
psi 


2 2 
Y V OK 
e e z 


A HI 

a -k kV 

z z z 




U 1 -! H- 

~K tz K 
z z 


y 1 

y 2 V \ k'k" 
'e e z z z 


(5.9) 


i e 

The significant contributions in C . come from the A„, A,, and 

psi I 7 I 7 

A® terms of (4.46). Using (5«9)> the coupled mode equations (4.52) can 


be reduced to the form derived by Litzenberger and Bekefi ( 1969 ). 


5.2.3 The Dependence of on £^(0) 

To determine the dependence of on £.^(0), we now have to deter- 

mine the temporal damping rate of an ion-acoustic wave due to ion-neutral 
collisions. This can be shown to be roughly v in /2 (Ohnuma and Hatta, 1970) 
The plasma parameters given in ( 5 . 1 ) show that the ion Landau damping 
rate (Ginzburg, 1970), 



is much smaller than v /2 , and negligible. Using (5.l), (5.2), (5.8), 
(5»9 )j and (4.63), we have the following values of the normalized 
parameters 


v ^2.9 X 10 
e ^ 


-3 


■ 1 . 4 x 10 ^ 


r ^8.3x 10 




^ 5.9 


(5.11) 


where we have taken Y = 

e 


5/3 as a convenient value. 



It remains to obtain the dependence of (J. /p and on £ (0)/£ 

n Q r\ ' ' f 4- "U * 


th 


From figure 5.3, we first plot the corresponding values of p /p and the 

P 

pump wave exciting voltage, V , as shown in figure 5.5(a). Assuming 

P 

£p( 0 ) “ figure 5 » 5 ( a ) indicates the approximate relations between 
Pp/p and and and P p /p, 


P t (0) 

^ -0.23 ^ + 0.77 , 

th 


M 


0 


'.. 3 - 3 . 3 ^- , 

P 


(5.12) 


where the second relation has been obtained from ( 4 . 66 ) by use of the 

first relation. Combination of (5.12) with (5.4) then determines the 

dependence of on £ ( 0 )/£y^, which is as shown in figure 5 . 5 (b). 

It is seen that z Q (— 2.5 cm) is in good agreement with the observed 

value of 2 cm, and its dependence on £ ( 0 )/£^ according to ( 5 . 4 ) is 

much weaker than that according to =- 0 . 7 /p. , which is obtained under 

the conditions (4.67) used by Ohnuma and Hatta ( 1970 ). Thus we conclude 

that, for 6 p( 0 ) > & th > and with an externally applied signal wave 

£ s » apparent independence of Z Q from £• ( 0 ) observed by 

Ohnuma and Hatta is likely to be a result of the increase in p with 

P 

V°) • 

Finally, it is of interest to examine the magnitude of £ for 

th 

the typical pump field strength indicated by Ohnuma and Hatta, which is 
Ep ~ 32 V/cm, corresponding to a ratio of rf to dc charge densities of 
roughly 1/20. For the typical pump and signal frequencies of 200 and 
120 kHz, we have, by use of (4.30), ( 5 . 1 ), and (5.8), 


n - 7.0 x 10 


n' - 4.2 x 10 



K' - 30 . (5.13) 
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Figure 5^5 ( a ) The increase of [x with the grid exciting voltage, V 

for the pump wave, as deduced from figure 5-3> and their 
approximate relation. (b) The dependence of the peak 
location, V of the idler f y on the boundary values 
of the pump, £^(0), according to (5^) and the relatio; 
Z. 0.7/d used by Ohnuma and Hatta ( 1970) • 



Then use of (4.6 l), (5*9)j (5.1l)> and (5.13) gives C th =“1.3 X lO - ^, 
corresponding to a pump threshold of 4.2 V/cm. This value confirms that 
the assumption £^( 0 ) > in connection with (5.3) and (5.4) is 

appropriate. 

5-3 Excitation of an Ion-Acoustic Wave by Two Electron Plasma Waves 

We will now examine the possibility of exciting an ion-acoustic wave 
by two electron plasma waves in the experimental plasma discussed in 
Section 5.2.1. Our purpose is to determine whether the nonlinear ly excited 
ion-acoustic wave should be observable under practical conditions, so that 
a measurement of the coupling coefficient can be carried out. 

5.3.1 Electron Plasma Waves (p) and Coupling Coefficient: 9=0 

Use of (5«5) and ( 5 . 6 ) in (4.43) and (4.44), with the assumption of 
^ ^ 1 j yields the dispersion relation for the electron plasma wave (p), 

C? - 1 + V. V , e P . { 0 , 0 , 1 ) ■ (5.14) 


This P wave dispersion relation is plotted, together with that of the 

ion-acoustic wave (i), in figure 5 >6, which shows that the interaction 

between the electron plasma (pump and signal) and ion-acoustic (idler) 

waves is of the contra-flow type according to the definition shown in 

figure 4.2(b). Because the slope of the P curve at V K ~ 0.1 is 

e z 

much larger than V . we have approximately K n 2K 

a z z z* 

The electron plasma wave action and group velocity can be obtained 
by using (5.5), ( 5 - 6 ), and <5.l4) in (4.50), giving 
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- r v v 2 

^2 e e ft 
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(5.15) 
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Figure 5 .6 Dispersion curves and synchronism conditions for longitu- 
dinal electron plasma (pump and signal) waves (p) and 
ion-acoustic (idler) waves (l). 
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The coupling coefficient can be obtained by using (5o5)-(5,8) and (5.1M 
in (4 0 46) and (^. 53 )* yielding 


C . 
psi 


P p I 

-K K'K" 
z z z 
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Y V V 


( 5 . 16 ) 


e e z 

Equations (5.15) and (5.1 6 ) may be used in (4.54) to obtain the temporal 
amplification rate, for the ion-acoustic wave as 



n"|e | 2 

' p 1 

4 q"Y v 2 


e e 


( tr 


1 <“7 \ 

A i ) 


which agrees with the result given by Kadomtsev ( 1965 ), who considered 

the case Y =1. 
e 


5.3*2 Evolution of and £ 

1 s 

As an example, let us consider an ion-acoustic wave excited with 

frequency 400 kHz and wavelength 1.1 cm, so that “ 1.4 x 10" 3 and 

K z ~ 100. The corresponding pump wave number is then K =“ 50. The 

z x 

Landau damping rate, for this electron plasma pump wave can then 

be obtained from ( 4 . 8 l ). It shows that 14 <0 is much smaller than the 

Li pe 

damping rate due to electron-neutral collisions, v / 2 , whose value is 

en 

given by (5.2) (Ginzburg, 1970) . Thus we have r P ^ 1.3 X 10~ 2 , By use 
of this result, and ( 4 . 63 ), ( 4 . 6 l ), ( 5 . 9 ), ( 5 - 15 ), and ( 5 .l 6 ), the para- 
meters in ( 4 , 68 ) and ( 4 . 69 ) fake the values 
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-|i ** 21 
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4 ± 5.9 , Z < 1.8 


v = 1.1 
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th 


2.1 X 10 


-4 
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2.1 X 10" 




0.1 . ( 5 . 18 ) 
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Since X £ |Mq| « 1 , the approximate solution of (4,70) can be used. 

When ( 5 . 18 ) is substituted, (4.70) becomes 

e i ~ [ e i(°) [ ex P(" 2 -l z ) “ 1 0" 3 exp(2.1Z)] + 0,021. e g (0)sinh(2.1Z) | exp(-3Z) , 

[e g (0)ex p (23Z) - 0.0019 e.(0)sinh(23Z)jexp(-3Z) , (5.19) 

-4 

where we have assumed ft (0) ^ 5,8 X 10 , corresponding to a field 

strength of 17 V/cm. 

The idler equation in (5*19) shows that if ft (0) < 50 where 

ft ( 0 ) = £ in is the thermal noise level for some acceptable measuring 

frequency bandwidth, no observable ft^(Z) above G^ n is possible. 

Assuming a measurement bandwidth of 2 kHz, the thermal fluctuations in 

the electric field of the ion-acoustic wave branch may be estimated roughly 

to be 1.5 X 10 ^ V/cm near 400 kHz (Bekefi, 19 66 ), To excite to 

above this level near Z ^0.2, i,e, 3-5 cm from the pump source, with a 

17 V/cm pump field, will require that £ g (0) > 2,5 X 10 , corresponding 

to a signal field strength of 7-5 mV/cm, According to the signal equation 

in ( 5 - 19 )j this requirement is equivalent to having ft ( 0 , 3 ) > 10 \ i.e. 

s 

injecting a 3 V/cm contra-flow signal wave at 2^5 cm. 

It remains to determine whether these parameters represent an acceptable 
experiment to demonstrate the interaction between collinear electron plasma 
and ion-acoustic waves. According to (4,82), the condition to allow trapped 

g 

electron effects on the pump wave to be neglected is ft^^3-^X 10~ , From 
(4,83), the condition to allow electron heating effects due to the pump 
wave to be neglected is ftp ^ 2,6 X 10 Also, from (4,85 ) > the 
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plasma will not become significantly inhomogeneous due to spatial variation 

of the pump wave amplitude, provided that & < V ~ 2.9 x 10~^ With 

p e 

-4 

the value of ~ 5 X 10 , it is seen that, first, the assumption 

fip « ©xp(-)ipZ) is inappropriate because of the trapped-particle effects, 
thus rendering (5-19) also inappropriate, and second, according to (4.83), 
the background electron temperature may become roughly 100 times the 
value given in (5.1) when the pump wave is injected, thus significantly 
changing the linear wave properties and nonlinear coupling coefficient 
used in this subsection. Also, due to the t.pnrnprafnra fu- 

thermal fluctuations in the ion-acoustic wave branch may become much 
stronger than the previously indicated value of 1.5 x lO”^ V/cm. As a 
result, we conclude that it will be difficult to demonstrate the interaction 
between electron plasma and ion-acoustic waves in the experimental plasma 
defined in Section 5*2*1. 

5*4 Excitation of an Ion-Acoustic Wave by Two Whistler s; Collinear 

Propagation “ 

This nonlinear process has been observed in an experiment by Porkolab, 
Arunasalam, and Ellis (1972). A detailed quantitative comparison between 
the theoretical results of Section 4 and the observed parametric amplifi- 
cation process in this experiment will be given here. It will be shown 
that there is wide disagreement between theory and experiment. The compari- 
son will suggest that important factors were overlooked in the theoretical 
model, and indicate possible directions in which the theory might be 
improved* 
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5.4.1 Experimental Observations 


The experiment by Porkolab, et al.(l972) was carried out in a 
steady state magnetized helium plasma column, which is illustrated 
schematically in figure 5*7* The plasma has the following typical 
parameters: 

n = 11 . ** 3 X 10 12 /cm 3 , T =“ 5 eV , T. =“0.1 eV , 

e i e 7 i 

6 6 

P = 2 mTorr , v 9-7 X 10 /sec j V. =“ 4.9 x 10 /sec , (5-20) 

7 en • m 

where v and v. have been determined by use of (5:2)* The effec- 
en in 

2 - 16 2 

tive electron-neutral collision cross section, (ora ^5 X 10 cm , has 

been obtained from the data given by Brown (1966), and the ion-neutral 

o -15 2 

collision cross section, (ita ) ^ 2 X 10 "cm , has been assumed to 

be 4(ita 2 ) . 

A large amplitude whistler pump wave ; with a typical frequency of 
2.45 GHz and a wavelength of 1.6 cm, was injected into the plasma with 
the dc electron cyclotron frequency set at 3*68 GHz. Their normalized 
values are then 

•fi^0.l6 , K ^ 1.4 , Q ^ 0.24 . (5.21) 

z c 

The pump wave variation was observed to be nearly exponential with a 

damping length of roughly 20 cm, corresponding to a normalized damping 

length of |i ~ 0.0l6. A typical interferometer output is reproduced 
P 

in figure 5*8, for a frequency of 7*5 MHz, i.e. 0* ““ 4*9 X 10~ . It shows 
that the excited ion-acoustic (idler) waves consisted of two wavelengths. 




Figure 5.7 The experimental set-up used by Porkolab, et al.(l972), 
in which ion-acoustic (idler) waves parametrically ex- 
cited by whistler (pump) waves were observed. The pump wave 
(jK; Q) was injected into the magnetized plasma column by a 
slow wave structure (SWS). The spatial behavior of the idler 
wave (K",ff) was measured by use of two rf probes (P 1 and , 
narrow-band-pass filters (f), and an interferometer (int.). 
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Figure 5 . 


0 2 4 6 

AXIAL POSITION z(cm) 


8 The interferometer output showing the ion-acoustic 
(idler) waves at 7.5 MHz, i.e* ^*4*9 X 10 
obtained by Porkolab, _et al* (1972)- The approxi- 
mate behavior of the longer wavelength idler wave 
for z ^ 2 cm is shown dashed, and was obtained 
by subtracting the shorter wavelength idler wave 
from the trace* 
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Additional measurements indicated strongly that the longer wavelength 

wave propagated at 0" =* 20° and the shorter at e" =“ 50°, where e" 

is the angle between k" and the dc B field. The wavenumbers were 

found to satisfy k' =“ -K^ and k" — 2 K^ in the first case, and 

K' ^ -k" and |.k"| » K in the second case. We shall consider the 

first case under the assumption q" ==“ 0 ° here, and defer discussion 

of the second case with 0" =“50° to Section 5 . 5 , 

It should be pointed out here that only the ion-acoustic wave in 

the second case was found to satisfy the dispersion relation o" =* v k" 

a ' 3 

so there is a reasonable doubt that the excited idler wave in the first 
case can indeed be considered as an ion-acoustic wave defined in an 
infinite homogeneous plasma background. We shall proceed by assuming 
that it is, and see whether the idler wave behavior shown in figure 5.8 
agrees quantitatively with the corresponding theoretical results. 

5^.2 Whistler Waves (R) and the Coupling Coefficient: 0 = 0 

Use of (5.5) and ( 5 . 6 ) in (4.43) and (4.44), with the assumption of 
fi c ~ 1 » a, yields the dispersion relations and polarization vectors 
for the right-hand (R) and left-hand (l) polarized transverse waves. 
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Q 2 - 


Q 

(flTQ) 

C 


2 1/2 [ 1)Ti)0 ] 


( 5 . 22 ) 


where the upper sign is associated with the R wave. For Q < n , the 

<r 

R wave corresponds to the whistler, whose dispersion relation is plotted, 
together with that of the ion-acoustic wave ( 1 ), in figure 5 . 9 , which shows 
that the interaction between the whistler (pump and signal) and ion- 
acoustic (idler) waves is of the contra-flow type. Because the slope 
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Figure 5*9 Dispersion curves and synchronism conditions for 
whistler (pump and signal) waves (r) and ion- 
acoustic (idler) waves (i). 
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we have approximately 


of the R curve at ^ 1 A is much larger than V , 

K." ^ 2K — -2K Since this condition was observed experimentally, the 
collinear approximation to the wave interaction process is seen to be 
appropriate. 

The whistler wave action and group velocity can be obtained by using 
(5*5), (5.6), and ( 5 - 22 ) in (4*50), giving 
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(5.23) 


The coupling coefficient can be obtained by using (5.5)-(5.8) and (5.22) 
in (4.46) and (4.53)> yielding 


r - r RRI - a LRI 
psi C K k'k" ~ -K K'K" 
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e e z 
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(5-24) 


where we have used the relation, e^. = -S_k * In ex P ress ion, the 

2 2 z z 

factor l/Y e V e k" is due to the factor 1/D" in of (5.5). If 

the ion-acoustic wave (i) were replaced by an electron plasma wave (p), 

this denominator would be replaced by f /' 2 - Y V 2 k " 2 =- n" 2 — l bv use 

(5*5) and (5*6) . The resulting coupling coefficient, A is 

— K K K 
z z z 

then identical to that obtained by Harker and Crawford ( 1970 ) for the 
interaction between two right-hand polarized waves and an electron 
plasma wave. 


5.^*3 Evolution of & and & 

1 s 

To determine the behavior of ft , we need to obtain the whistler 
and the ion-acoustic wave damping rates. The whistler wave is damped by 
electron-neutral collisions (T q ) and,, and by cyclotron resonance (T^) when 



Q is close to fi . They are given by (Ginzburg, 1970) 
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pe/ 2f2(Q-Q )^+ 0 
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r r “ * V K exp 


e z 
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. V K 
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. (5.25) 


Use of (5.20) and (5.2l) in (5-25) and (5.10) shows that the damping of 


the whistler and the ion-acoustic waves is predominantly due to electron- 
and ion— neutral collisions, respectively. Then the first expression of 
(5-25), and the relation f 1 “ V ln /aD pe ^ & ive 


r R “ 6.7 x io" 5 , r 1 2.5 x io 


( 5 . 26 ) 


R 

Use of (5.20), (5.21), and (5.23) gives the group velocities U =“ 7.2X10 

and U 1 =. V .^..4.7 X 10” 5 .. .Use of these results and (5.26) .in (4.63) then 
a 

gives 

H -- 9 . 2 X io -4 , |U ~ 0.53 , (5.27) 

S 1 

where u <0 because U « -i . Note that I fj- 1 is much smaller than 
r s ~z s 

the observed spatial damping rate for the whistler pump^ u =* 0 . 0 l 6 ; 

ir 

showing that the pump wave is heavily damped nonlinearly. 

It is also necessary to determine the values C . and g, before 

psi th 

the theoretical behavior of & can be predicted. Using (5.2l), (5.23), 
(5.24), ( 5 . 26 ), and (4.6l) then gives 

c , =- 1.8 x 10 6 , e,, 1.0 X 10~ 5 , (5.28) 

psi th 


where the corresponding unnormalized pump threshold is 1 6 V/cm. We 
expect this to be smaller than & ( 0 )j which can be determined 

approximately from the input power and cross-sectional area of the plasma 
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at the pump wave exciting structure (Lisitano, Fontanesi, and Sindoni, 

2 

1970 ; see figure 5*7). Taking 0.8 kW and 20 cm , respectively, and 

s 

using (5.23), we find £ (0) ^ 7.2 X 10 , which is smaller than t 

P th 

of (5.28). Therefore, if both v and v. given in (5.20) were not 

en in 

overestimated, the theory predicts no parametric excitation of the ion- 
acoustic wave. 

To remove this discrepancy, let us assume that = T =“ 3.4 X 10"^ 

s 9 

as opposed to (5.26). Then the threshold field is reduced to 

g 

^th ^ 7 X 10 . Use of this 6^, the reduced values of p 
(^ -4.6 X 10 ^), (5.26), and the observed value of (=* 0.0l6), in 

( 4 . 69 ) then gives 


v ““ *7 t |m q | “ 2.0 . ( 5 . 29 ) 

Because v » |m q | ^ X/ the approximate solutions of (4.72) can be used. 
With ( j U g | 3 g /U^5^) ““ 0.l4 for this case, (4.72) becomes, 

t ± (z) =“ ^( 0 ) [exp(- 0.53 Z) - 0.0038 exp(- 0 . 0 l 6 Z)J 

+ O.OO 85 £ g ( 0 ) [exp( - 0 . 0 l 6 Z) - exp(- 0.53 Z)J . ( 5 . 30 ) 

The observed ion— acoustic (idler) wave has a nearly monotonic variation 
from 0 to 6 cm from the pump wave exciter. According to figure 5.8, the 
ratio of to the apparent noise amplitude, t changed approximately 

from =1.5 at Z = 0 to = 2.0 at Z =“ 19 . It is clear 

that this observed £ behavior is substantially different from that 

indicated by ( 5 . 30 ). 
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The discrepancy between theory and experiment for this case may be 

due to the effects of the finite plasma column size, the noncollinear wave 

propagation (i.e. Q" — 20 °), and some other nonlinear processes, which • 

were neglected in the foregoing theoretical model. In a more realistic 

model, the parameters used in (4.68) and (4.69) may be modified by these 

effects. It will be noted, in particular, that because p^ » p^ » | p g | 

in the case of interest here, fi^Z) depends most strongly on p. 

(through v) . If we reduce, p^^ by a factor of 4 from (5.27) to p. =*0.13 

and I\ = 0.6 X 10 with p g ( = -4.6 X 10 ) already reduced by a 

factor of 2, then we have .£ (0) = 2& th , v =4.7, and |m q | =2.0. Using 

these parameters, (4.72) can be used to estimate £ and £ as 

1 s 

e 4 (Z) = e i (0) jexp( -0 , 13 z) - 0.061 exp( -0.016 z)] 

+ 0.038 & s ( 0 ) jjexp( - 0 . 0 l 6 z) - exp(- 0 .l 3 Z)J , 

£ s (Z) = £ s (°)exp ( 4.6x 10 _ 4 Z) . (5.31) 

Equation (5>3l) gives a good match to the observed £ with £.(0)-£ (0):£ 

i 1 ' s w in 

given approximately as 1:80:1. It is also seen from ( 5 . 31 ) that for 

Z ^ 19 , we have £ g (z) = £ g (0) ~ 80 ^ n * To determine if this modification 

in p. is appropriate, it would be necessary to measure £ (z)/£ in 
^ s in 

the corresponding experiment. 

5*5 Excitation of an Ion-Acoustic Wave by Two Whistlers: Oblique 
Propagation 

In Section 5 * ^ it was shown that the observed behavior of the idler 
ion-acoustic wave with wavelength 0.8 cm does not correspond to the 
collinear three-wave interaction theory in an infinite homogeneous plasma 


/ 
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background. Here, we shall examine the theoretical behavior of £, 
corresponding to the second observed idler ion-acoustic wave with 
wavelength roughly equal to 0.l6 cm (see figure 5*8)* Typical wave- 
numbers of the idler ion-acoustic and signal whistler waves were found 
to be 

k" 15 , e" =* 48° , lk , e'^it+ 52 0 . ( 5 . 32 ) 

Since obliquely propagating waves are involved in this wave interaction 
process, the linear characteristics of these waves may be very different 
from their corresponding waves with parallel propagation. Therefore, we 
shall proceed with a discussion of these obliquely propagating waves. 

5.5-1 Obliquely Propagating Ion-Acoustic Waves (l") , k" = {k^O, k" | 
_ ~ _±__ y 

In the hydromagnetic. region (ft" « 0 3-3 X 10 -5 ), the linear 

ion-acoustic wave dispersion properties may be significantly affected by 
the presence of a sufficiently strong dc magnetic field. For the experi- 
mental plasma of interest here, the Alfven speed, V. =“ (n n =» 

A ' i c' 

2.8 X 10 , is much larger than the ion-acoustic speed, V =- ^.7X 1 0~^ . 

The ion-acoustic wave dispersion relation^ and the polarization vector in 
the hydromagnetic region, are then (Ginzburg, 197 O) 

Q,/2 V a 2K#2 * e" - [0,1,0} . (5.33) 

Here e is linearly polarized perpendicular to k" and B. 

The ion-acoustic wave observed in the experiment, however, had a 
frequency Q ^ 4.9 X 10 , and therefore does not fall into the hydro- 

magnetic region. It is necessary to examine the properties of the wave 
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2 a 2 2 

in more detail. Under the condition (5.6), and with Cl^ » Q » 
we have from (4,4o) 
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2 2 

Although a , b « 1 in (5-34), they are included because the zero order 
terms are found to cancel each other in the derivation of the dispersion 
relation and the polarization vector. Substitution of (5.3*0 in (4 .42) -( 4 .44) , 
we then have 
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(l/V a 2 HK" 2 /V A 2 ) 


v 2 k " 2 


e 


revealing a nearly linearly polarized longitudinal wave, since fi ~ 15 . 

z /2 2 „2 2 2 2 

The approximation Q" V & K* , is appropriate because 1/V a ~ 30 /V . 
When Cl' « we note that je w reduces to (5-33) of the hydromagnetic 


case. 


We may now substitute (5.3*0 and (5-35) in (*K50) to obtain 


YX 2 (lY 2 
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(5-36) 


e e 
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When B" = 0, (5.36) reduces to the corresponding expressions in (5.9). 
Also, we have Y. = 0 since U. is parallel to K , 

5.5.2 Obliquely Propagating Whistlers (r'), K' = [k% 0, K '} 
^ X z 


Reduced forms of the plasma polarization constant, M S , can be 

obtained from (5.5) by use of ( 5 . 6 ), and the condition, q ' 2 » Y V 2 k ' 2 

e e 7 
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. (5.37') 


The dispersion relation and the polarization vector for the oblique 
whistler wave may be obtained from (4.42)-(4.44) by use of ( 5 . 6 ), ( 5 . 37 ), 
and the conditions □ ' <■ Q c and Q^O^sinS ' « 4 (q' 2 -i) 2 cos 2 0 ' (see 
for example, Helliwell, 1965 ). We have 
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. (5.38) 


When 0'= 0, e' reduces to e** of (5-22). Since K ' 2 Q c » Q ”, we 
have je =* {sin 0 , 0 , cos 0 '}, which represents a linearly polarized 
longitudinal wave, corresponding to a whistler wave that is near the 
oblique cyclotron resonance fi' - -Q cos 9 ' ( 3n/2 > 9 ' > ji/2) . 

Use of (5.37) and (5.38) in (4.50) now gives 
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(5.39) 


Here, the second equation for 3 g is obtained by replacing Q by 

2 

-Q^cosG '( 1 - 1/K' ), from (5.38), before eliminating terms involving 

p 

1/K' as small compared with unity. From (5.39 )j we see "that the direc- 
tion of U is perpendicular to k' and oriented so that f = 

5*5.3 Synchronism Conditions and Coupling Coefficient 

Synchronism conditions: Since the signal and idler waves propagate 
obliquely, the synchronism conditions become less obvious than those 
cases shown in figures 5*^, 5.6* and 5*9* The synchronism conditions for 
this case are shown schematically in figure 5*10* It contains the dis- 
persion curve R from (5.22), for whistler waves with parallel propaga- 
tion; the dispersion curve R / from (5*38), for whistler waves with 
oblique propagation near 0 ' ~ rt + 52 °, and the dispersion curve l" 
from (5.35) (plotted with respect to a displaced origin O'), for ion- 
acoustic waves with oblique propagation near 0 7/ ==“ 48°. Since the whistler 
signal on the R * curve is nearly at the oblique cyclotron resonance, we 
have 


Q ^ Q' =“ 0 cose' » Cl" 

C 


(5.^o) 


RR 

Coupling Coefficient; To obtain the coupling coefficient, * ir /t , 

■ KK. K 

LR ^1 77 R I ^ 

(4.53) shows that we have to evaluate h , u, since e - e . This 
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Figure 5 . 10 The synchronism conditions for interaction among a whistler 
(R) with parallel propagation as the pump (K, Q), an 
Obliquely propagating whistler (R ') as the signal (k',Q')> 
and an obliquely propagating ion-acoustic wave (i") as 
the idler (K",n"). The i" curve is plotted with respect 
to the displaced origin, 0% to show that the conditions, 

~ = iS + 2$' and n= ft' + Ci"f are satisfied. 
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requires M 1 and M & for a left-handed circularly polarized wave, which 
can be obtained from ( 5 . 36 ) as 


M = M 




I Q 0 


M = M 


* 




(5.4i) 


Then use of (5.22), (5.34), (5-35), (5-37), (5-38), and (5*4l) reduces the 

r R ' T » 

expression for Ay™, * v u from (4.46) to a simpler form. With the paiameters 
KK K 

given in (5.20), ( 5 . 32 ), and (5.40), the significant terms are, in order 
of ‘decreasing importance 
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The corresponding coupling coefficient, is 'th® 11 
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(5.43) 


where we have made the additional simplifications to A _ KR * R // , appropriate 

to the experimental conditions of interest, K ~ K » K^, 9^0 ^ 0 - it, 

2 2/2 

fi =“ fi" » fi", and fi » Y e V g K . The first term in the square brackets 
of (5.43) is the combined contribution, of the first two terms of (5.42), 

Note that (5.43) cannot be reduced to (5.24), for the coupling coefficient 
among whistlers and ion-acoustic waves with parallel propagation, by 
making 0=0. This is because the polarization vector, je , in (5*38) for 
the obliquely propagating whistlers, is nearly a linearly polarized longi- 
tudinal wave, rather than a circularly polarized transverse wave of (5.22). 
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5.5.4 Evolution of and t 

i s 

For interaction of obliquely propagating waves, the coupled mode 
equations of (4.79) are required. We make X = 0 because only along the 
axis of the plasma column were measured ( Porkolab, et al., 1972). After sub- 
stituting the appropriate values from ( 5 . 20 ), ( 5 . 2 l), and ( 5 . 32 ), for 

(°“ 0 . 0 l 6 ) and 0 (^ 50 °); and using the experimental pump wave 

-6 

amplitude, £p(0) = ‘7.2 X 10 , we obtain from (5-36), (5-39), (5.43), 

and (4.69), 


V*s« “ “3.7 


3 i U iz “■ 2k > c = 4.3 X 10- , 


M. 


‘O' 


c & (0) 

p 


%< 3 S l U sJ 3 i U l z > 


1/2 


21 


(5.44) 


The expression for |Mq| 1® modified from (4.69) by replacing U and 

v s 

\ by U sz and respectively, following the definitions of ( 4 . 77 ). 

To obtain p iz , we must use (5*27) and = 0" =“ 50° in (4.77). 

To obtain p gz , an expression for the damping of an obliquely propagating 
whistler is required. It is (Ginzburg, 1970 ) 
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Then use of (5.45), (5-39), and (4.78) yields the value for u . 
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(5.45) 


-Q^cosS 0 . 
We have 


4 iz ^0.34 , ja gz =- - 6.8 x 10 


(5.46) 


By use of (5.46), the value of (i (“* 0.0l6), and (4.69), we find 
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(5.^7) 


v = 




{1 + ) 

SZ XZ y 


12 . 


According to figure 5.8, we are interested in the region where z ^ 2 cm, 

which corresponds to Z < 5, and hence \i Z « 1* Therefore, from 

P 

( 5 . 44 ) and ( 5 . 47 ), we have X ^ |Mq| > v » 1 , so (4.7*0 can be used to 
approximate and Substitution of ( 5 . 44 ), ( 5 * 46 ), and ( 5 . 47 ) 

in ( 4 . 74 ) then gives 


£ i (z) - ^(0)003(0.34 Z) + 0.39 £ s (0)sin(0.34 Z)jexp(-0.l7 Z) , 

£ g (z) = [e- s (0)cos(0.34 z) - O .39 6^(0) sin (0. 34 Z)]exp(-0.17 Z) , (5-48) 


which have been plotted in figure 5.11* Note that good agreement of £ 

with the experimental observations has been obtained without changing the 

value of \i . given in (5*27). This indicates that this value of jj. p is 
i i 

acceptable as the collisional damping for short wavelength ion-acoustic 
waves ( K." =“ 15). For K " ** 3> however, the results of Section 5.4 
suggest that the long wavelength ion-acoustic wave may have been strongly 
modified by those factors in the experiment not included in the wave 
theory of an infinite homogeneous plasma background. Finally, to make the 
agreement more satisfactory, & g (z) and £ should be measured in 
the experiment to see whether the theoretical result in (5*48) is also 
acceptable. 
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Figure 5.11 Comparison between the observed evolution of the short 

wavelength idler ion-acoustic wave, deduced from a result 
obtained by Porkolab, et al. ( 1972) (figure 5.8), and the 
theoretical evolution of according to ( 5 . 48 ), with 

e i (0):e s (0) 

:£ in given approximately as 5 : 4 : 1 . The 
theoretical evolution of the signal whistler wave amplitude, 
^ s ( z )^ is subject to future experimental verif ication. 
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5 .6 Discussion 

The verification of linear plasma wave theory imposes only the 
necessity for verifying the dispersion relation, i.e. measuring 
K and ft. In the case of nonlinear wave-wave interaction, it is 

y~> 

necessary to establish not only that the linear waves satisfy the dis- 
persion relation, but that also that the synchronism conditions are satisfied 
and coupling coefficients follow theory. This requires more extended 
series of measurements, providing details of the amplitude variations 
of the pump, signal, and idler waves with distance. Few satisfying 
comparisons between theory and experiment have been made so far. 

We have tried to remedy this difficiency to some extent by predicting, 
in Section 4, the wave amplitude behavior in space, taking into account 
such practical factors as damping of the pump, signal, and idler waves. 

In this section we have tested this theory against available experimental 
results (Ohnuma and Hatta, 1970; Porkolab, e^t a]_*, 1972) . Unfortunately, 
these experiments do not provide all the data required for a complete 
quantitative comparison with theory, since they are confined to a qualita- 
tive presentation of the idler wave behavior without specifying the 
observed values of pump threshold, noise intensity, and, in one case 
(Porkolab, et al, 1972), the signal wave behavior. 

The four examples treated have all involved the excitation of ion- 
acoustic waves, either by other ion-acoustic waves, electron plasma waves, 
or whistlers. In Section 5.2 we have shown satisfactorily that the 
observed independence of the location of the idler wave peak from the ion- 
acoustic pump wave amplitude is due to the increased pump wave damping 
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with amplitude. In Section 5- 3, we have examined the difficulty of exciting 
an ion-acoustic wave to observable amplitude by injecting two counter- 
streaming electron plasma waves into a weakly ionized plasma when 
appreciable collisional damping occurs ■ In Section 5.4, the predicted idler 
wave behavior has disagreed with the experimental results. We have shown 
that the discrepancy can be explained if the effective ion-acoustic wave 
damping was substantially lower than that due to electron-neutral collisions. 
Further theoretical and experimental work is necessary to determine 
whether this discrepancy is caused by inhomogeneous plasma background, 
oblique wave propagation, or by other nonlinear wave processes 0 Good agree- 
ment has been obtained in Section where the excitation, of an obliquely 

propagating ion-acoustic wave of short wavelength was considered. The 
results have served to demonstrate the need for extensive and detailed 
measurements to be taken, in which wave and noise amplitudes are thoroughly 
documented, so that detailed comparisons can be made against theories such 
as those of Section 4. 

It will be realized by the reader that the averaged Lagrangian method 
is more efficient and versatile than the conventional iterative method 
for the description of nonlinear wave-wave interaction processes. The 
results of Sections 4 and 5 can be readily applied to wave interaction 
processes other than those studied in this section. 
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6. CONCLUSIONS 


The main aim of this thesis has been the establishment and applica- 
tion of Lagrangian methods to magnetoplasmas described by macroscopic 
equations* This involves derivation of a Lagrangian of which the 
given macroscopic equations are the Euler-Lagrange equations obtained 
by applying Hamilton's principle* We have taken two approaches to this 
problem. In Section 2 , we considered the general mathematical inverse 
problem of the calculus of variations, i.e. the derivation of Lagrangian 
densities from an arbitrary set of equations. We were able to establish 
sufficient conditions for systems of first and second order, quasilinear, 
differential equations, and used these conditions to transform equations 
apparently not in Euler-Lagrange form to Euler-Lagrange form. As examples, 
appropriate Lagrangians were obtained for a linear resistive transmission 
line, and for a linear, collisional, one-dimensional^ warm plasma. 

For reasons summarized in Section 2.5> it was found that substantial 
further development is necessary before this general mathematical approach 
can be used to obtain Lagrangians for plasmas. Consequently, in Chapter 
a suitable Lagrangian density was obtained for the macroscopic plasma 
description through energy considerations, taking account of the dual 
roles played by the dynamic variables, i.e. the macroscopic plasma cell 
position and electromagnetic potentials. The effects of viscosity, heat 
conduction, and elastic collisions were included by energy balance 
arguments. To obtain the corresponding Hamiltonian, the canonical 
momentum conjugate to the plasma cell position was defined in Lagrangian 
coordinates, while the canonical momentum conjugate to the electromagnetic 
vector potential was defined in Eulerian coordinates. The resulting 
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Hamiltonian was shown to equal the appropriate macroscopic plasma energy 
at t =s 0* For completeness, the Hamiltonian corresponding to the Low 
Lagrangian (Low, 1958), appropriate to the microscopic description of 
plasmas, was obtained, and shown to equal the appropriate microscopic 
plasma energy at time t. It was found that, for a collisionless (Vlasov) 
plasma, the macroscopic Lagrangian could not be obtained by a simple 
velocity integration of the Low Lagrangian, This fact was traced to 
the loss of particle discreteness in the macroscopic approximation. 

In Section 4, the macroscopic Lagrangian was expanded in terms of 
small perturbations for the case of scalar pressure and adiabatic compres- 
sion, and used to study nonlinear three- wave interactions in a homogeneous 
magnetoplasma. The averaged Lagrangian method was applied to obtain the 
coupled mode equations. These were extended to include phenomenologically 
the effects of wave damping. Solutions of the coupled mode equations were 
given describing the spatial variations of the signal and idler waves 
under the assumption of a strong but damped pump wave. 

In Section 5, these solutions were specialized to experimental 
conditions involving parametric excitation of ion-acoustic waves by 
other ion-acoustic waves, by electron plasma waves, and by whistlers. 
Quantitative comparison with available experimental data demonstrated good 
agreement for situations in which the propagation characteristics of the 
uncoupled waves could be measured accurately. Suggestions were given in 
Section 5.6 for improved theory and experiments to reduce the most 
significant discrepancies between them. 

Although Lagrangian techniques can be used to obtain complete solutions 
to plasma wave problems, such solutions may not always be necessary. For 
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example, in determining the oscillation properties of a bounded plasma, 
the variational properties of the action integral can be used to obtain 
the resonance frequencies, using only approximate trial functions for 
the electromagnetic fields and plasma dynamics. We have treated one such 
problem in Appendix B. The quadratic Lagrangian obtained in Section 4 
was specialized to the problem of electrostatic resonances in an inhomo- 
geneous plasma column. The Rayleigh-Ritz procedure was applied directly 
to the Lagrangian. For a low pressure positive column, it was found 
that accurate frequencies and eigenf unctions could be obtained efficiently 
for the first few resonances,, provided that appropriate coordinate functions 
were defined, in contrast to numerical solutions to this problem obtained 
by Parker, Nickel, and Gould (196^), the variational approach was found 
to be applicable for the entire range of the ratio (column radius/electron 
Debye length) without incurring serious numerical instability in the 
calculation of the first few resonance frequencies. 

Some extensions to the work described in this thesis have already 
been discussed in the individual sections. Among, and in addition to> 
these are the following: 

First, the general mathematical approach of the inverse problem of 
the calculus of variations was found to be less effective in practice than 
the more intuitive approach by energy considerations used in Section 3. 

One of the reasons for this is that the mathematical approach lacks an 
independent definition of the dynamic variables that represent the degrees 
of freedom in the corresponding physical problem. It would be of interest 
to consider a problem in which the dependent (generalized) variables are 
defined by the physical degrees of freedom before employing the general 
results of the inverse problem. We have not treated the question of 
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subsidiary constraints. When these are imposed, the results of Section 2 
must be modified to include them. 


Second, one of the long-standing difficulties associated with varia- 
tional principles is that of including energy dissipation effects in a 
physical system (see for example, Goldstein, 1959). We have only 
partially succeeded in eliminating this drawback: in Section 2, because 
of the difficulties associated with solving first order, nonlinear partial 
differential equations, the examples on dissipative systems were limited 
to linear equations; in Section 3, by closing the system of energy transfer 
in the macroscopic plasma, the effects of viscosity, heat conduction, and 
elastic collisions, were included in the Lagrangian. These effects had 
to be discarded in Section 4.2, however, because of difficulties involved 
m obtaining the corresponding perturbation approximations. The phenome- 
nological approach used in Section 4.4.2 includes the wave damping effects 
in the coupled mode equations, rather than the Lagrangian leading to them. 

It would be of interest to see if dissipation effects could be included 
in the Rayleigh-Ritz procedure used in Appendix B. A possible procedure 
has been suggested by Mikhlin (1964). 

Third, the results of Section 4.3 were written completely in terms 
of quantities representing the linear wave properties: the polarization 
vector e, the plasma polarization constant M S , the frequency Q, and 
the wavenumber K.. in the case of nonlinear wave-wave interactions, the 
averaged Lagrangian method is equivalent to expressing the approximate 
Lagrangians in terms of the linear eigenfunctions, e exp[i(QT-K • x) ] , 
of the homogeneous plasma. The Manley-Rowe relations, and the synchronism 
conditions, can then be viewed as results of the orthogonality property 
of these eigenfunctions. Extensions of the averaged Lagrangian method 
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to inhomogeneous plasmas should then he straightforward: we would start 
with the general energy conservation theorem derivable from the approximate 
Lagrangians obtained in Section 4 * A formal substitution of the eigen- 
functions, in place of the dynamic variables, could then be made* The 
orthogonality property of the eigenfunctions should lead to selection 
rules for the appropriate eigenvectors, and a power balance relationship 
between the interacting linear inodes* In those cases where the exact 
eigenfunctions are difficult to obtain, the approximate eigenfunctions 
could be obtained by the Rayleigh-Ritz method used in Appendix B. 
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APPENDIX A t SOLUTION OF ULTRA HYPERBOLIC EQUATIONS 


To obtain the solution of the homogeneous part of (2.26), 




( A . 1 ) 


when cc / and i 4 j, we first convert it to ultrahyperbolic form 
(Koshlyakov, Smirnov, and Gliner, 196^), 


2 2 p p 

U + O. £2. 0 

2 2 2 p “ u 

aS l ^ S 2 at l Bt 2 

by the substitutions, 


(A. 2) 



s +t 
11 





S 2 +t 2 



V*2 




*"?< «J, d“ II?) . Jt. 1( s 2 , t 2 ) . 


(*-3) 


Equation (A. 2) is similar to the wave equation in a homogeneous two- 
dimensional medium, but with two-dimensional time coordinates. Its 
general solution can readily be obtained by separation of variables. 
Assuming 


9 = v i( S i) V 2^ S 2 ) W 

leads to the system of ordinary differential equations, 


(A. 4) 


l 8 l 



d 2 v 

1 


,2 

d V 2 


.2 

d w 


a v 
11 

a 2 
d w' 

* 2 
ds 2 

“ a 2 V 2 * 

dtf 2 


2 _ 

b w , 

a„ + a . — 

b, + b,. 


dt 2 2_ 

2 2 ’ 

1 2 

1 2 


where a ^ and b^ are constants. The explicit solutions of 

etc. are sinusoidal functions. Because no boundary conditions or 
initial values restrict the solution for the general solution of 

(A. 2) takes the form. 


7 - j* da i da 2 db l C (VV * 


where the weighting function, C(a^^ag } b^), will be restricted by the 
sufficient conditions of ( 2 . 31 ), ( 2 . 32 ), and (2.34). 
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VARIATIONAL CALCULATIONS FOR RESONANCE OSCILLATIONS 
OF INHOMOGENEOUS PLASMAS 

*>y 

Y-K „ M. Peng and F, W. Crawford 
Institute for Plasma Research 
Stanford University 
Stanford, California 94305 

ABSTRACT 

In this paper, the electrostatic resonance properties of an inhomo- 
geneous plasma column are treated by application of the Rayleigh-Ritz 
method. In contrast to Parker, Nickel, and Gould (1964), who carried 
out an exact computation, we have used a description of the rf equation 
of motion and pressure term that allows us to express the system of 
equations in Euler-La grange form. The Rayleigh-Ritz procedure is then 
applied to the corresponding Lagrangian to obtain approximate resonance 
frequencies and eigenfunctions. An appropriate set of trial coordinate 
functions is defined, which leads to frequency and eigenfunction esti- 
mates in excellent agreement with the work of Parker, et al. (1964), 



1. INTRODUCTION 


This paper is concerned with the use of the Rayleigh-Ri tz procedure 
to estimate the electron resonance frequencies of a warm inhomogeneous 
plasma column. This procedure has been extensively applied to single self- 
adjoint equations with great success (Mikhlin, 1964), For a system of 
equations, however, theoretical extensions have been noted only for the 
case of elliptic equations (Mikhlin, 1965), For the equations to be 
used here, which are not elliptic, it will be shown that accurate resonance 
frequencies can be predicted, provided that a certain set of coordinate 
functions is defined. 

Previous theoretical treatments of the electron resonance problem 

give predictions which agree well with experimental data. However, these 

approaches have encountered difficulties stemming from the inhomogeneous 

electron density profile. For example, Parker, Nickel, and Gould (1964) 

solved numerically an appropriate fourth-order differential equation for 

the rf potential in a cylindrical positive column. Because of the exponential 

nature of the solutions in the cutoff region, their calculations were 

2 2 

limited by the condition, a /\^ < 4500, where a is the column 
radius, (- e 0 KT / n 0 e ) is the mean-squared Debye length, and 

is the mean electron density in the column. Baldwin (1969) used a 
kinetic model in the low-temperature approximation to obtain the external 
admittance of a cylindrical plasma capacitor. The appropriate differential 
equation was solved after using inner-outer expansions connected through 
the resonance region. A WKB description was used for the inner region, 
where resonant waves are essentially evanescent, while a travelling wave 
description was used for the outer region, where Landau damping is important. 
Because the wave nature of the solutions was assumed in the outer region, 
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this theory is appropriate only for higher order resonances. Similar 
difficulties have been shown, to occur in the simpler one-dimensional model 
(Harker, Kino, and Eitelbach, 1968; Miura and Barston, 1971; Peratt and 
Kuehl, 1972). 

Variational methods offer an attractive alternative to these treatments. 
They have been used previously to estimate plasma resonance frequencies 
with simplified trial functions. Resonances of a cold inhomogeneous plasma 
were treated by Crawford and Kino (1963). Using the variational principle 
established by Sturrock (1958), Barston (1963) approximated the dispersion 
relations for wave propagation along an infinite cold plasma slab, and 
along the interface between two semi-infinite, counter-streaming cold 
plasmas. Some general features of the guided waves on a cold, transversely 
inhomogeneous plasma column in an axial magnetic field were studied by 
Briggs and Paik (1968). These papers (Crawford and Kino, 1963; Barston, 

1963; Briggs and paik, 1968) show that, with appropriate variational 
principles and judicious choices of trial functions, useful results can 
be obtained with relative ease by the variational approach. 

A theoretical variational formulation for the electrostatic resonance 
oscillations of a warm, inhomogeneous plasma column in a dc electric or 
magnetic field of arbitrary direction was presented by Barston (1965) 
with the adiabatic index, y , taken as unity. The variational prin- 
ciple to be presented here, however, is not restricted in the values of 
7 . One important feature in Barston 's (1965) analysis is that the 
rf electric potential was treated as the solution of the rf Poisson 
equation, with the rf electron density considered given. It will be 
seen that the coordinate functions to be used here are defined in a 
similar fashion. A variational method of the Rayleigh-Ritz type has 
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been applied successfully by Dorman (1969) to a one-dimensional, warm, 
and field-free plasma with arbitrary dc density profile. A single 
second order differential equation for the electric field was obtained, 
and shown to have hermitian operators. The variational principle to 
be used here differs from Dorman’s (1969) in that we are dealing 
directly with a system of Euler-Lagrange equations. In so doing, we 
can keep down the order of the equations, and are able to consider 
warm inhomogeneous plasmas in more than one dimension. 

In this paper, we shall show that by appropriate definitions of the 

rf equation of motion and pressure term, the equations of the hydrodynamic 

model used by Parker, et al . (1964) become Euler-Lagrange equations. This 

will enable us to demonstrate the effectiveness of the Rayleigh-Ritz 

procedure in estimating the resonance frequencies of an inhomogeneous plasma. 

The associated numerical method mainly involves evaluations of definite 

integrals and solutions of finite algebraic eigenvalue equations, and is 

2 2 

applicable over the entire range of a » 1 for estimating the first 
few resonance frequencies. 

In most of the papers that deal with the electrostatic resonance 
problem (Crawford and Kino, 1963; Parker et al . , 1964; Harker et al. , 1968; 
Baldwin, 1969; Dorman, 1969; Miura and Barston, 1971; Peratt and Kuehl, 

1972), it is assumed that the rf plasma current normal to the glass wall 
is zero. However, in the low temperature limit, a -♦ co , the main 
resonance frequency seems to agree with that of cold plasma theory, in which 
the normal rf plasma current is retained. We consider this problem and 
show that the resulting difference in predicted resonance frequencies is 

/ 

negligibly small for low pressure positive columns. 
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In §2, we present the basic equations, the corresponding Lagrangian, 
and the procedure to be applied in the variational approach. In 53, the 
numerical methods are explained before comparing computations with those 
of Parker et al. (1964). The paper concludes with a brief discussion 
in §4. 
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2 . THEORY 


For a low pressure positive column, moment equations with scalar 

pressure and negligible heat conduction are appropriate when the wave 

phase velocity is much larger than the thermal speed. For the first few 

electrostatic resonances, the wave phase velocity may be scaled to yy a , 

where u> p [ = (e n 0 (°)/me 0 ) ' ] is the axial plasma frequency. Thus we require 
2 ~2 

a /^d >:> 1 4 ^ stationary ion background will be assumed, since we are 
interested only in electron resonances. Dissipation due to collisions, 
and Landau damping, will be neglected. Our analysis will consequently be 
valid only for the first few resonances. Also, the analysis will be quasi- 
static. Apart from some differences in definition of the rf equation 
of motion and pressure term, the equations we shall use are essentially 
those used by Parker, ^et a_l 4 (1964). The equations are generalized here 
to include dc magnetic field, , and electron drift velocity, 

We have, 


V ’ <nv) = 0 , 


€ 0 9 'E + e ( n_n I > - 0 » 


mn (|t + v-Vv| + <?P + en ( E + v x B) = 0 ( at ly + i) • CD 

Specialized to small perturbations, these reduce to the dc equations, 


9-UoV = 0 , 


e o 7 ‘£o + e <n 0 _n i ) = 0 ’ 


ran oV^O + * P 0 + en 0 ( V^0 XP 0> = 0 (at r 0 ) 


( 2 ) 


and rf equations, 


^ n i 

jr + '’•‘"oil + “iV = 0 • 




mn 0 ? + (V-1)VP 0 + VP X + l-nVP Q ) 

+ en o ( S/ 7 £o + h + i X S 0 + Zq * V = 0 (3) 


6 



where j* = d£/dt = BJ /dt + . In these equations, m and -e are 

the electron mass and charge; n and P are the electron density and 
pressure; E is the electric field; is the ion density; is the 

vacuum permittivity; and ^ is the perturbation displacement for the 
electrons (Figure 1). 

The magnitude of is relatively small in the plasma region, but 

increases in the sheath region from the ion-acoustic speed to roughly the 
electron thermal speed at the glass wall (Self, 1963; Parker, 1963). We 

consequently neglect it in our analysis. However, due to the presence 
of non-zero v^, jn^, and an rf electric field at the wall, a non-zero rf 
normal current term arises, and hence an rf surface charge term. The 
inclusion of this surface charge term, in the cases where the electron 
rf excursion exceeds the Debye length, is equivalent to the use of the 
dielectric model for a cold plasma column (Crawford, 1965). Fui’ther 
discussion of this surface charge term will be given in §3.4. With v 
neglected, the following relations become appropriate 

Vi) = n 0 (0)f( i: ) ’ f Cr) = ex P[-e<P 0 <i:)AT e ] , 

(4) 

n i = ’ A = a£/ at > 

where k is the Boltzmann constant, ^(r) is the dc column potential, 

V£> = “*P 0 <£> ’ < 5 > 


and the first equation of (3) has been used to obtain . The rf 
electron pressure, , is determined by the adiabatic equation of state, 


P( VP/W - ■ P 0 - V\ ' (6) 

The form of (6) can be understood by reference to Figure 1, and 


follows from the fact that when a cell is displaced from r to r + F 

-0 -0 ’ 
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P is defined for the given cell, rather than for the given position, r 

— 0 

The adiabatic equation of state must consequently be applied to the same 
cell, before and after the displacement. Using the usual definition of 
perturbations , 

P <£> *= P o ( £ ) + > n <jO = h Q (r) + n^r) , (7) 

and (4), we obtain 

P 1 = 1*^0 • (8) 
The rf force law in (3) is obtained by comparing the force laws in (1) 
and (2) in the same fashion (Newcomb, 1962). 

It is now straightforward to use (4), (5), and (8) to rewrite (3) 
in terms of only J? and the rf potential, ) , 

“"ol ” <7 -D?p 0 V'£ - v|-vp 0 - yP Q v (?■$_) 

- en c|( ttp 1 + i-vvcp 0 - IxBq) = 0 , 

e 0 v \ + * 0 . (9) 

Equations in (9) can be Fourier-transformed, normalized, and expressed 
in a cylindrical coordinate system, (r,6), for a column of cylindrical 
symmetry. 


oV r - n o 0fe 9 * a* [r(C> I C e > ♦ 5<« 8 -C r >] £' 


Q 


* A D 7'K * 1 C r - 5 * I<} oCr +f 'l> = 0 • 

X - ♦ A It' |ac r -c 9 ) * 4 7f i(c; + Ur-Ue) 

* f (1 *o + 5 *i) ' 0 • 


( 10 ) 


( 11 ) 
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( 12 ) 


T *1 + < R£ Cr>' ' C e = 0 . 

where the derivative with respect to R is denoted by (') and 

V 0) = 0 ’ 4 =lr c r + ie (iC e ) ’ 

00 f 

$ i (R,e,T) ( R)exp i(OT+£0) , 


(13) 


with i denoting unit vectors. The normalized quantities are defined as 
R = r/a , T = u) t , O = w/<» n , 

r r 

Q c = eB^mujp , C - S/a , 

A D = X d/ & 2 = € o KT e/ n o (0)e2a2 ' $ = cpe 0 /n 0 (°)ea 2 , (14) 

and a static axial magnetic field, , has been included* 

2 . 1 Lagrangian Density 

The forms of the force law in (3), and the rf pressure in (7) , 
represent the important differences from the paper by Parker, et al . (1964) 
in that they make (9),as well as (10)- (12), systems of Euler-Lagrange 
equations without having to restrict the values of y (Barston, 1965)* 

In one of the models used by Dorman (1969), an appropriate pressure term 
similar to (7) was used without the benefit of the rf force law in (3). 

As a result, he was able to establish the variational principle only for 
the one-dimensional case. 

The Lagrangian corresponding to (10)- (12) can be shown to be the 
following, by straightforward application of Hamilton's variational principle, 


10 





1 ) = n A - ft QB + H , 
* c f 


{,“00 0 


X 

A J f RdR (|e r | 2 ♦ | Ce | 2 ) , 

J r\ 


B = / fRdR k C* + CiC 


r b 9 ' 


H = -A D (?T+T'> - V Q S - I+F+F / < 0 


X 

T ,J f »dB| C ; + | c r - I g 2 * f [fr ? <c; A r r - 1Ce ; 
1 

T ' -/ tRdR [i * 5 Ce«Q - C*e> * <=•■=•] 


+ 2 t"-f«C e -C r > ♦ £ CS«-C p -s; 9 > + ■ 


1 

.y fRdR («'ig 2 + i(»'| Cf 




y RdR (|^| 2 + i! | §i | 


F ' =/ £ g MR (l*l ! 2 + ^ Kl 2 ) A / BdR (U'l 2 + £ lij 

n R „ ' R 


R, = ~ R _ c 

b a 1 


11 



where c,c, denotes complex conjugate. Equation (15) is appropriate 

to the configuration shown in Figure 2 of a concentric metal cylinder 

surrounding a glass tube, of relative permittivity e , that contains 

g 

the plasma column. In the expressions for T and T / , the boundary 
terms are included to modify the natural boundary conditions on r 
and (see for example, Courant and Hilbert, 1953) that would other- 

with be unphysical. In the expressions for H and S , 0^(r) and 
are defined as 


*0 <r) = - <f 0 (r)/V 0 , V Q . - 4 0 (») . 

Substitution of the exact solutions of (10)- (12) would make £ (Q,<L) 
zero, We see from (15) that, for negligible Q , A and V , the 

C u u 

resonance frequencies are determined essentially by the values of I, F, 

and F / . Since V Q is approximately proportional to (Self, 1963; 

Parker, 1963), the effect of higher electron temperature is to raise each 

resonance frequency. When Q ^ 0 , the roots of £, (Q,£) = 0 are 

c 2 

0 12 = n c B/2A ± [(Q c B/2A) 2 - H/A] 1 / 2 . (16) 

Since fi and - Q are indistinguishable in experimental observations, we 
.see that all the resonance frequencies are predicted by (16) to split 
in two, in agreement with the theoretical results of Barston (1965), and 
Vandenplas and Messiaen (1965). For sufficiently small 0^ , when the 
values of B and A are not greatly affected by the presence of an 

axial static magnetic field, the amount of- the split, | Q^B/A ] , will 

be proportional to This is in agreement with the observed splitting 

character of the main dipole resonance frequencies (Crawford, Kino, and 
Cannara, 1963; Messiaen and Vandenplas, 1962), Furthermore, since 
As | B | , where the equal sign applies when = Q e , this split is predicted 
to be always less than or equal to Q 
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Rayleigh-Ritz Procedure and Coordinate Functions 


For a single linear Euler-Lagrangian equation, the Rayleigh-Ritz 
procedure is efficient in obtaining an approximate solution, by use of a 
weighted summation of a set of judiciously chosen coordinate functions, 
These coordinate functions must be linearly independent and complete, 
and satisfy the boundary conditions specified by the problem. The 
weighting coefficients that appear in the approximated Lagrangian are 
varied independently. This results in a system of algebraic equations 
that take the place of the original equation. Theoretically, better 
approximations can be obtained by using more coordinate functions. When 
eigenvalues are involved, the approximate eigenvalues always converge to 
the exact values from above (see for example, Mikhlin, 1964). 

Much less attention has been given to the analogous problem for a 
system of linear Euler-Lagrangian equations, e.g. (10)- (12), though the 
theoretical extension of the variational method to a system of second order 
elliptic equations has been mentioned by Mikhlin (1965): the way to set 

up the corresponding coordinate functions is similar to that for the single 
equation case, i.e. the coefficients of each dependent variable are varied 
independently . 

For our problem, in which (10)- (12) are not in elliptic form, the 
coordinate functions and the coefficients must be more restricted. In the 
Appendix to the paper it is shown that acceptable estimates of resonance 
frequencies and eigenfunctions can be obtained for our problem provided 
that the coordinate functions chosen for each dependent variable are 
related by (10)-(12). 

By expanding , and $ in power series of R, and substituting 

in (10)- (12), we see that for small R , r °c £ a R^ -1 , and $ oc 

^ * Q 1 
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Thus, for £ £ 1 t the solutions are well behaved at R = 0 , If we 
also choose even, functions for f(R) and $^(R) , then and are 

even in R , and is odd in R , when l is odd, and vice versa , 

Since there are no other singular! t ies in (10)- (12), polynomials in R 
constitute appropriate coordinate functions for our problem. For convenience, 
the coordinate functions chosen for will be 

r . = R*'" 1 - / +2J_1 (j = 1,2,...) , (17) 

^rj 

which conform to the usual assumption of zero normal rf current, since 


w 1 ’ - ° • 

Rather than choosing r . and § . independently, we must determine 

"6j 1 , J 

them via the original differential equations and (17). After eliminating 
in (10) and (11), it follows by using the second equation of (4) that 


+ c 6 ■ = Lr Q c Q + Cj + C r -^C e ) • (18) 

With a given expression for $ Q (r) , and an assigned value of ft , e.g. 

0=1, Cqj can then be easily determined for any given £ - 

An immediate question arises concerning the dependence of the resulting 

variational estimate of resonance frequencies on the size of Q chosen 

arbitrarily here. We have found that the first few resonance frequencies 

-4 

do not change by more than 10 m when ft in (18) changes from 0.4 to 1 

for all the values of 1 /Aq used in this paper. If this were not the 
case, an iterative procedure would have to be used, i.e, the resulting 
variational estimate of ft would have to be used in (18) to obtain a 

new set of . These would be used in turn to obtain improved frequency 

estimates, and so on. 

The corresponding coordinate function, $> , , is obtained by solving 

I , J 

(12) (Barston, 1965), with conditions of continuity of potential and normal 
displacement across the boundaries defined in Figure 2. According to the 
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discussion given in the Appendix, identical coefficients are assigned to 
each set of coordinate functions , 


a j C rj ’ ^0 ' *1 " 2 a jVd ’ <1£ 

j=l j-1 j=l 

before substitution in the Lagrangian, £ r (0,-t) , of (15). The resulting 

2 

Lagrangian then gives the algebraic Euler-Lagrange equation below, 


(fl A.. - «i B.. + H )a 4 = 0 
J1 cjl jl j 


(i = 1,2 , ... , N> , 


I 

S. j = / fRdR (C ri G r j + ^'91^0 j) * 
J 0 


where A. . , B. . , and H. , are the matrix elements of the integrals, A, 
l J 1 3 ^ 3 

B, and H, given in (15), respectively, and are obtained by substituting 
the coordinate functions in a fashion given by the above A. . expression, 

3. j 

Equation (20) can be transformed into a generalized eigenvalue problem 
(Dorman, 1969) 


» >1 


( ^- w ° 


(i = 1,2, ... , 2N) , 


b . = (a. ,Qa. ) 
2 i l 


(j = 1,2, . . ., 2N, i = 1,2, . .. , N) , 


/-H 

0 \ 



-H V 


U 

A + ) 

5 

-u 

a b + ) 

C ' 

(21) 


with superscript + signifying the transposition of a matrix. Equation (21) 
is now solvable by standard computer codes . 
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3. NUMERICAL METHODS AND RESULTS 

The computing procedure is straightforward: 

(i) read in physical and computational parameters > 


(ii) compute coordinate functions , and all other functions 

appearing in the Lagrangian, * at intervals Ar , 

(iii) compute (^>&) by Simpson’s rule to obtain A„ , , and 


H . . , and 
ij 


(iv) solve (21) for ft and a^ , and compute relevant eigenfunctions for 
N, N-l, and N-2 coordinate functions. 


3.1. Approximate DC Density Profile 

A density profile which approximates Parker's results (Parker, 

1963), and is convenient for both analytical and numerical manipulation, 
is given by 


f(R) = exp [- T) w <j> 0 (R) ] , <t> 0 (R) = pR 2 + (l-p)R h , (22) 

where h(> 2) is an even integer, p < 1 , and T| w [= V^/A^] is the wall 
value of the potential function, T)(r) , used by Self (1963) and Parker 
(1963). The particular form of (22) is used to justify the choice of 


coordinate functions of (17). Furthermore, with the use of (22), Q 

’ 6j 

3nd h,j Can n ° W be solved analytically in terms of power series in R, 

in addition to the numerical solutions of (12) and (18). Comparison of 

the solutions by the two methods will provide an estimate of the degree 

of accuracy achieved in obtaining Q q . and $ . 

The values of p and h are varied until f(r) best approximates, 

by least-square deviation, the profile given by Parker for specified values 

° f and ’ The resultin e profiles are shown in Figure 3, and the 

corresponding values of p and h are given in Table 1. It will be 
2 

seen that as 1 /A d increases in value, (22) decreases in accuracy because 
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densit 



Table 1 

Parameters used in (22), and Figure 3, for a mercury- vapor 

2 2 

positive column, The conversion between 1 /a and 1A 

1 D D 

is obtained by the calculations of Parker (1963). 
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of the increasing degree of steepness displayed by the density profile 
in the sheath region. At zero temperature, the sheath is omitted and 
an accurate approximation can again be obtained. The approximation of 
(22), however, was found to be sufficient to give accurate frequency 
and eigenfunction estimates for all of the values of l/A^ listed in 
Table 1. 
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3.2, Solutions for r and $ 

By the use of (14), (15), and (19), Cqj can be P ut in the 
form of a rapidly converging series. 


r . = R 




-1) (h-2)-l 


Q+(£+2j 


j 


i-2) - 


2 6 (7-1) V f 


f- (S,+2j-2) + ttti-2 6 (7-1) AV 0 


<7-l)(l-p)hV [£a -(^+2j-2)] 

-n J ; 

Mj >2 ^2 

^2(-t+2j+h-4)+£-) C2-2p( 7 -lHV n 
l co 


^ r( 7 -i>(i-s)hv 0 

J>1 iQ 2 (h-2) 


(i > 3). 


Since the maximum value of of interest is roughly 0.02 and 

^ 1 , one needs at most five or six terms in (23) to attain a precision 
of 10 ® for C Q . . 

To solve for . , the predictor-corrector method of Adams-Bashforth 

> 3 

(Fox, 1962) has been used on (12), which can be reduced to the form, 

y[ = y 2 / R > y 2 = g(R) + & 2 y^/R , 


y i = Vi ’ 


v « R4> 
y 2 K *l,j ’ 


g(» = f(-R\ * 0 'c rj + RC r 'j + C rj - *c ej > - (2« 
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Since the complementary solution of (24) is = R^ , and its 

P £+2 j 

particular solution, y , is proportional to R near r =: 0 , 

p' p 7 

y (0) and y (0) are both equal to zero. So the starting values of 
1 2 

P P 

y and y are well-behaved, and easily obtained by Taylor series 

X A 

expansion near R s= 0 . The total solution of $ can then be written 

i f J 

as 


■ °/ + *i,j • 


(25) 


where c . is determined by imposing the boundary conditions of $ . 

By making the interval A r =r 0,01, and using double precision, $ 

x , j 

— 8 

can be calculated to within 10 , This is arrived at, first, by com- 

paring results that use different values of A r > and secondly, checking 
against solutions of (12) obtained by power-series expansions in R . 

3 . 3 Numerical Instability 

In the process of solving the algebraic equation, (21), the size 

of N is limited by the inaccuracy involved in obtaining A_ , etc. 

This inaccuracy introduces a numerical instability whenever the coordinate 

functions are not orthogonal functions with respect to the differential 

operators of (10)-(12) (Mikhlin, 1971, Chap. 2). The situation is best 

illustrated by an example in which the dipole resonance frequencies 

corresponding to Figure 3(b) are calculated for different values of N 

while holding the size of A r constant at 0.05, 

As shown in Table 2, as N is increased from 2 , the first few 

resonance frequencies are approached from above with rapidly stabilized 

estimates. When N is increased beyond 8, undesirable fluctuations 
-4 

larger than 10 , and clearly erratic changes in the values of Q , start 

to appear. In the case N = 9, for example, one would obtain an erroneous 
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Table 2 


Dipole resonance frequency estimates obtained with £r = 0.05 
for the case of Figure 3(b). Significant numerical insta- 
bility sets in when N £ 8 . The ’best estimates* are 
obtained with £r = 0.01 and N = 9 . 


N 

Main 

First 

Second 

Third 

Fourth 

Best estimates: 

0.4447 

0.7099 

0.8746 

0.9924 

1.080 

2 

0.4703 

0.8813 




3 

0.4486 

0.7446 

0.9250 



4 

0.4450 

0.7155 

0.8818 

1.022 


5 

0.4448 

0.7119 

0.8757 

0.9958 

1.084 

6 

0.4448 

0.7117 

0.8753 

0.9956 

1.084 

7 

0.4448 

0.7104 

0.8751 

0.9947 

1.082 

8 

0.4448 

0.7098 

0.8748 

0.9921 

1.081 
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fundamental resonance frequency. Characteristic of the variational 
nature of the Lagrangian, ^(0,-t), more serious errors are found in 
the approximate eigenfunctions, than in the resonance frequencies. The 
optimal combination of N and £r , that produces acceptable results 
in the shortest computation time, can be obtained by trial and error. 
Repeated solution of (21) for a few adjacent values of N thus becomes 
an economical technique; this requires computation of the matrices A 

ij ’ 

etc. only once, and offers safeguards against obtaining erroneous results 
due to numerical instabilities. 

3.4 Computer Results 

As a practical example, the approximate density profiles given 
in Table 1 have been used to predict dipole resonances for Tube No, 1 
used by Parker, ^t jal. (1964) (£ = 1; a = 0,5 cm; effective relative 
permittivity at the surface of the column K = $ '/£$ = 2.1). The 

61 I 11 

o 

computation time varies roughly as N . With N = 10, a typical calcu- 
lation takes about 40 seconds, and requires a core space of less than 
100K bytes in an IBM 370/67 machine. 

The resulting approximate solutions for | r >i| e » n^ , cp^, and 
are plotted in Figures 4-6. The density and the radial electric field 
solutions, n^(R) and Cp' resemble very closely those given by Parker, et al . 
(1964), and Parbhakar and Gregory (1971), respectively. The relative 

amplitudes of and shown in these figures are retained, revealing 

~ 2 

that as 1 /A d increases, g progressively dominates over g . For 
- — 0 

V^d > 4500 ' it wil1 be seen that the Perturbations should be progressively 
compressed toward the sheath region as l/A^ increases in value. These 
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solutions are not reproduced here because they also exhibit undesirable 
oscillations with wavenumber equal to N , an expected characteristic 
when we try to approximate rapidly varying functions with truncated 
polynomials. Since only a moderate computer storage is used for N £ 10 , 
there is room to increase N, and decrease Jr, to obtain better approxi- 
mate solutions. However, this is considered unimportant for our purpose, 
since we ai'e able to obtain good frequency estimates for this region with 
N £ 10, as Figure 7 reveals. 

The corresponding estimates of resonance frequencies are shown in 

Figure 7. Since the electron temperature corresponding to the experimental 

resonance data of Parker et al. (1964) was adjusted to fit their theoretical 

spectrum, it would be reasonable for us to make a similar adjustment* As 

is evident from Figure 7, however, no such adjustment is necessary. Indeed, 

our result seems to be in slightly better agreement with the t - 3 eV 

e 

data. The minor differences between the two theoretical results probably 
come from the differences in the rf equation of motion and pressure 
term used in the two treatments 

Similar to other papers (Crawford, 1964; Parker et al * 1964; Harker 

et al, 1968; Baldwin, 1969; Dorman, 1969; Miura and Barston, 1971; Peratt 

and Kuehl, 1972), we have assumed zero normal rf plasma current density 

at the glass wall, through the form of Q . in (17). This is appropriate 

t* j 

when the plasma is sufficiently warm that the electron excursion velocity 
is much smaller than the thermal speed, and f(l) « 1 , This assumption, 
however, is inconsistent with the dielectric model for a cold plasma 
column, where normal rf plasma current must be included. It is of 
interest to ask why the main resonance frequency of a wa^m plasma column, 
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Approximate solutions for 1/X^ = 510, l = 1 


27 






+ 



FIG, 7. Dipole resonance spect 


estimates by the variational 





in the limit of low temperature, approaches the principal resonance of 
a cold plasma column. To answer this question, we need only change (17) 
to 


C rj = R * / ~ 1 + 


(j — 1,2,...)., 


(26) 


and impose the requirement of continuity of normal displacement in the 
form 


+ f(DC r (D = 6 §'(1 + ) . C 27 > 

8 t 

The resulting solutions of are found to be only slightly different 

from the previous case near R = 1 [Figures 4-6, where the dashed lines 
correspond to the use of (17)]. Furthermore, the main resonance is 
lowered by less than 1 per cent for all of the values of l/A^ used 
here, including the case 1/A^J . This is well within the experi- 

mental errors . 
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4. DISCUSSION 


In this paper, we have applied the Rayleigh-Ritz procedure to a system 

of three Euler-Lagrangian equations that describe the electron resonances 

of a nonuniform warm plasma column. It is shown that accurate frequencies 

for the first few resonances can be obtained for the entire range of 
2 

1/Ajj >:> 1 • Results which agree closely with those of Parker, et al . 

(1964) have been obtained 

Contrary to the case of a system of elliptic equations, where the 
coefficients are assigned independently to each dependent variable 
(Mikhlin, 1971), we have found that for (10)-(12), the same coefficient 
must be assigned to each set of coordinate functions, e.g. (19). In 
addition to the usual requirements , that the coordinate functions must be 
linearly independent and complete, we have chosen that they be set up in 
accordance with (10)-(12), 

The present method can be easily modified to include the effects 
of electron dc drift, dc magnetic field, and ion motion. With the axial 
dimensions and rf magnetic field included, this procedure would be efficient 
in solving travelling wave problems in a nonuniform plasma waveguide. 
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APPENDIX 


Here we shall show that the coordinate functions, r . , L. , and 

b rj b 0J 

$ . , must be assigned the same coefficients, a, /, for the Rayleigh- 

>,J J 

Ritz procedure applied to £ (Q,£) of (15) to be successful: it will 

be shown that, by restricting these coordinate functions according to 
(10)-(12), the appropriate Rayleigh-Ritz procedure can be established. 
A system of second order elliptic equations can be written as 


'VP ' + X V = 0 • V ■ % ’ 


E,.. - E. . , 
jk kj 


(A.l) 


th 

where U is the k dependent variable, the matrices D and E 

k jk j K 

are functions of the independent variable x , the eigenvalue is 0), 

and the summation convention has been used. Ellipticity demands that 


D jk“A > 0 ’ B jk°A 


^ 0 , 


(A. 2) 


for any real non-zero vector cr . The solutions of (A.l) then admit 
of variational estimates, as outlined by Mikhlin (1965). 

When we apply the Rayleigh-Ritz procedure to the Lagrangian for (A.l), 




dx (D., U.U, - \ E., U.U, ) 
jk ] k jk j k 


(A. 3) 


the coefficients preceding the coordinate function for each U can be 

k 

varied independently. Suppose a set of legitimate trial functions, 

(summation convention not used here) with coefficient , are 

used in (A.3). Because of (A. 2), we can obtain crude estimates of \ 
larger than the lowest eigenvalue, even if only one of the C is non-zero, 

.K 

If the same procedure is applied to the Lagrangian in (15), we will 

obtain erroneous estimates of Q . Consider the case of a cold, uniform 

plasma column with B_ = 0 , so that H=-I+F+F / . Suppose C , 

0 r 
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C . and C are the coefficients for the coordinate functions £ . * 

0 9 0 

r ^ , and $. . , respectively* Then if C = 0 , the value of the Q 

*r0 ’ i,J 0 

estimate will be zero, and fail the requirement of the Rayleigh-Ri tz 


sequence . 

Consequently, to make (-H/A)^* from (15) at least non-zero, we 
must use a single coefficient for each set of the coordinate functions 


>V2 


r . , r . , and 6 . . Merely making C = C = C = a . is insuffi- 

cient to produce a legitimate Rayleigh-Ritz sequence, because the value 

of , for example, can be arbitrarily small in comparison with 

1 , J 

and £ , making the Cl estimates also arbitrarily small. We 

see that restricting these coordinate functions according to (10)- (12), 
is sufficient to reduce the resulting approximated Lagrangian to a single 
Euler-Lagrange equation. The appropriateness of the resulting Rayleigh- 
Ritz procedure can then be guaranteed* It is conceivable, of course, 
that there may be less restrictive choices of appropriate coordinate 
functions corresponding to the Lagrangian in (15), but we have not chosen 
to pursue this point* 


36 



